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ABSTRACT 

In this paper, we describe the design and data analysis of the DEEP2 Galaxy Redshift Survey, the densest 
and largest high-precision redshift survey of galaxies at z ^ 1 completed to date. The survey was designed to 
conduct a comprehensive census of massive galaxies, their properties, environments, and large-scale structure 
down to absolute magnitude Mb = -20 at z ^ 1 via ^ 90 nights of observation on the Keck telescope. The 
survey covers an area of 2.8 deg^ divided into four separate fields observed to a limiting apparent magnitude 
of Rab = 24.1. Objects with z ^ 0.7 are readily identifiable using BRI photometry and rejected in three of 
the four DEEP2 fields, allowing galaxies with z > 0.7 to be targeted ~ 2.5 times more efficiently than in a 
purely magnitude-limited sample. Approximately sixty percent of eligible targets are chosen for spectroscopy, 
yielding nearly 53,000 spectra and more than 38,000 reliable redshift measurements. Most of the targets which 
fail to yield secure redshifts are blue objects that lie beyond z ^ 1.45, where the [O ii] 3727 A doublet lies 
in the infrared. The DEIMOS 1200-line/mm grating used for the survey delivers high spectral resolution 
(R ^ 6000), accurate and secure redshifts, and unique internal kinematic information. Extensive ancillary data 
are available in the DEEP2 fields, particularly in the Extended Groth Strip, which has evolved into one of the 
richest multiwavelength regions on the sky. 

This paper is intended as a handbook for users of the DEEP2 Data Release 4, which includes all DEEP2 
spectra and redshifts, as well as for the DEEP2 DEIMOS data reduction pipelines. Extensive details are 
provided on object selection, mask design, biases in target selection and redshift measurements, the spec2d 
two-dimensional data-reduction pipeline, the specld automated redshift pipeline, and the zspec visual redshift 
verification process, along with examples of instrumental signatures or other artifacts that in some cases remain 
after data reduction. Redshift errors and catastrophic failure rates are assessed through more than 2000 objects 
with duplicate observations. Sky subtraction is essentially photon-limited even under bright OH sky lines; we 
describe the strategies that permitted this, based on high image stability, accurate wavelength solutions, and 
powerful b-spline modeling methods. Summary data are given that demonstrate the superiority of DEEP2 over 
other deep redshift surveys at z ^ 1 in terms of galaxy numbers, redshift accuracy, sample number density, and 
amount of spectral information. We also provide an overview of the scientific highlights of the DEEP2 survey 
thus far. 

Subject headings: Surveys - galaxies: distances and redshifts - galaxies: fundamental parameters - galaxies: 
evolution - galaxies: high-redshift 
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1. INTRODUCTION 

Spectroscopic redshift surveys have been a major gateway 
to understanding galaxy evolution. Locating each galaxy in 
space and time, they provide the raw material for compiling a 
census of galaxy properties as a function of cosmic epoch and 
position in the cosmic web. Clustering measures derived from 
redshift surveys of sufficiently large volume - such as the cor- 
relation function, ^(r^,!!) diagram, and measures of the local 
overdensity of galaxies - probe the gravitational growth of 
structure. Pair counts contain information on merger rates, 
while satellite motions probe dark-halo masses and dynam- 
ics. When the galaxy sampling is sufficiently dense, groups 
of galaxies become well defined and allow the estimation of 
the number density of groups (constraining cosmology), as 
well as the determination of the environments in which in- 
dividual galaxies are found. The spectra themselves can re- 
veal a wealth of detail on emission-line strengths, stellar pop- 
ulations, star-formation rates, AGN activity, gas and stellar 
metallicities, internal motions, and dynamical masses. Fi- 
nally, spectroscopic redshifts provide the fundamental data 
needed for calibrating photometric redshifts, which are the 
only way to estimate distances for the hundreds of thousands 
of fainter galaxies that are beyond the reach of current redshift 
surveys. 

Pioneering low-redshift surveys out to z ^ 0.2, such as the 
Sloan Digital Sky Survey (SDSS, York et al. 2000; Abajazian 
et al. 2003) and the Two Degree Field Survey (2dF, Colless 
et al. 2001) have demonstrated the value of mammoth sam- 
ples of galaxies containing hundreds of thousands of objects. 
Huge samples enable finer slicing in redshift and galaxy- 
parameter space, reveal rare phases of galactic evolution, and 
provide sufficient weight to measure statistical scatter about 
the various galactic scaling laws. 

Redshift surveys of more distant galaxies at z ^ 1 pro- 
gressed more slowly, primarily as they have targeted galax- 
ies more than 100 times fainter than SDSS or 2dF. The first 
substantial surveys were conducted in the 1990's, measuring 
redshifts of several hundred objects at intermediate redshift; 
e.g., the LDSS survey (Colless et al. 1990), the ESO-Sculptor 
Survey (Bellanger et al. 1995, Arnouts et al. 1997), the CNOC 
and CNOC2 surveys (Yee et al. 1996, 2000), and the Hawaii 
Deep Fields Survey (Cowie et al. 1996). 

Deeper surveys also began taking shape at roughly the same 
time. The pioneering Canada-France Redshift Survey (CFRS; 
Lilly et al. 1995) garnered ^ 600 redshifts and was the first to 
provide a dense, statistical sample of galaxies out to z ^ 1 . 
Though conducted with only a 3.6-m diameter telescope, it 
had a fairly deep magnitude limit of /ab = 22.5 and yielded a 
median redshift of z = 0.56. Augmented by Hubble imaging, 
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it established landmark norms for the luminosities and sizes 
of galaxies and the star-formation history of the Universe to 
z-1. 

The Keck telescopes were also used for deeper surveys. 
Cowie and collaborators (Cowie et al. 1996) used roughly 
400 redshifts down to / ^ 23 in the Hawaii Deep Fields to 
elucidate the phenomenon of "downsizing", that the peak of 
star formation occurred first in the most massive galaxies and 
swept down later to smaller galaxies. This concept has since 
become central to understanding galaxy evolution. Later sur- 
veys by the same group have probed the nature of distant X- 
ray sources (e.g., Barger et al. 2005, 2008). 

The Caltech Faint Galaxy Redshift Survey (CFGRS, Co- 
hen et al. 2000) exploited the Keck 10 m telescope and went 
roughly one magnitude fainter than CFRS (R = 23-24), as- 
sembling several hundred redshifts with median z = 0.70 in 
two fields (one of them containing the Hubble Deep Field-N). 
Notable discoveries included the extreme amount of cluster- 
ing found even at z ^ 1, the relatively normal morphologies of 
most field galaxies observed (i.e., the rarity of true peculiars 
and mergers), the measurement of the shape of the luminos- 
ity function to deeper magnitudes, and the exploration of the 
nature of distant mid-IR and radio sources (van den Bergh et 
al. 2000; Hogg, Cohen, & Blanford 2000; Cohen 2002). 

The largest of this first generation of deep redshift surveys 
was DEEPl. This survey was conducted using the Keck I 
telescope to a magnitude limit of /? ^ 24 in and around the 
HDF-N and the Groth Survey Strip (GSS), a 127 arcmin^ re- 
gion observed with WFPC2 on HST in both / and V (Groth 
et al. 1994; Vogt et al. 2005). DEEPl was conceived and exe- 
cuted as a pilot survey for the DEEP2 Galaxy Redshift Survey 
described in this paper. Over 720 redshifts were measured in 
total (Phillips et al. 1997, Weiner et al. 2005) having a me- 
dian value of ^0.65, and bulge and disk parameters were fit- 
ted to the Groth Strip HST images (Simard et al. 2002). Sci- 
ence highlights include an early survey that helped to establish 
the nature of Lyman-break galaxies (Lowenthal et al. 1997), 
the detection of a strong color bimodality out to z ^ 1 (Im et 
al. 2002, Weiner et al. 2005), the measurement of fundamental 
plane evolution for distant field spheroidal galaxies (Gebhardt 
et al. 2003), the detection of "downsizing" in the star forma- 
tion rate of compact galaxies (Guzman et al. 1997), measure- 
ments of ISM metallicity evolution to z ^ 1 (Kobulnicky et 
al. 2003), an extensive analysis of selection effects on disk- 
galaxy radii and surface brightnesses (Simard et al. 2002), 
and the discovery of the prevelance of already-red bulges in 
typical disk galaxies at z ^ 1 (Koo et al. 2005; Weiner et 
al. 2005). A distinguishing feature of DEEPl was the use 
of relatively high spectral resolution (R ^ 3000, FWHM = 
3-4 A), which resolved redshift ambiguities by splitting the 
[O ii] A3727 doublet and also yielded linewidths and rotation 
curves. These were used to extend the Tully-Fisher relation 
to z -- 1 (Vogt et al. 1996, 1997), a strategy that was later 
used in the Team Keck Redshift Survey (GOODS-N; Wirth et 
al. 2004) to yield an even larger sample of linewidths (Weiner 
et al. 2006a,b). 

DEEPl was a valuable learning experience. In addition to 
proving the scientific value of high spectral resolution, it con- 
firmed that faint spectral features could be more easily seen 
in the dark stretches of spectrum between the OH lines. On 
the other hand, achieving photon-limited sky subtraction un- 
der the newly concentrated OH lines places even higher re- 
quirements on flat-fielding accuracy, wavelength calibration. 



and night-sky spectral modeling. Furthermore, pushing be- 
yond the de facto redshift limit that then existed near z ^ 1 
(the so-called "redshift desert") would require pursuing [O ii] 
A3727 to well beyond 8000 A. A new spectrograph on one of 
the largest optical telescopes would be necessary to achieve 
these goals. 

The present paper describes the design and execution of 
the DEEP2 Galaxy Redshift Survey (commonly referred to as 
DEEP2 in the remainder of this paper), which was conducted 
by a collaboration of observers primarily based at UC Berke- 
ley and UC Santa Cruz, together with astronomers from the 
University of Hawaii and other institutions. The survey was 
executed with the DEIMOS multi-object spectrograph (Faber 
et al. 2003) on Keck 2; aspects of the survey have previously 
been described in Davis et al. (2003); Davis, Gerke & New- 
man (2005); Coil et al. (2004b); and Willmer et al. (2006). 
The survey obtained roughly 53,000 spectra, coming close to 
the target of ^ 60,000 galaxies in the original survey design. 

DEEP2 was conceived as a legacy survey for the astronomi- 
cal community, and all data are now public. The present paper 
is intended as a handbook for users to understand and utilize 
the data products in Data Release 4 (hereafter DR4), which is 
available on the web at http://deep.berkeley.edu/DR4/. Brief 
descriptions of the data and methods are given there, but the 
present account is more complete and presents new details 
on the final data reduction techniques, redshift measurement 
procedures, the current generation of redshifts, and summary 
data tables. Additional measurements that were not part of 
previous DEEP2 data products are also included. All data in 
this paper are available in the form of FITS BINTABLE files 
at the DR4 website. 

The outline of this paper is as follows. Section 2 sets forth 
the overarching goals and tradeoffs that determined the de- 
sign of the DEEP2 survey. Section 3 summarizes some of 
the major science results so far, giving examples of how vari- 
ous aspects of the design have translated into specific science 
gains. Section 4 gives the basic DEEP2 survey parameters 
and field properties and describes the extensive existing an- 
cillary data that are available in the DEEP2 fields. DEEP2 is 
put into context in Section 5, which compares DEEP2 to other 
leading z- 1 surveys (VVDS-deep, VVDS-wide, zCOSMOS, 
and PRIMUS) using several different measures of survey size 
and information content. 

Section 6 describes the parent galaxy photometric catalogs 
and selection cuts that were used to determine potential tar- 
gets for DEEP2 spectroscopy. The slitmask design process 
is covered in Section 7, which describes the DEIMOS detec- 
tor and slitmask geometry, as well as the algorithms used to 
place objects on masks. Section 8 is an exhaustive discussion 
of all known selection effects resulting from the slitmask de- 
sign process, caused either by loss of objects from the initial 
target pool or by systematics in the rate at which we select 
pool objects of given properties for masks. 

Sections 9 and 10 describe the spectroscopic exposures and 
their reduction, including detailed information on the spec2d 
pipeline that produces the reduced 1-d and 2-d spectra. This 
pipeline has been employed by many DEIMOS users; this 
paper provides its definitive description. The features of 
the DEIMOS spectrograph and specid pipeline that enable 
photon-limited sky subtraction accuracy even under the OH 
lines are described in the Appendix. The range in quality of 
DEEP2 data is discussed in Section 10, where examples of 
typical data, bad data, and interesting spectra are given. 



Section 1 1 describes the redshift measurement process, in- 
cluding the automated specld redshift pipeline that produces 
a menu of redshift possibilities and the subsequent visual in- 
spection process to select among them via the interactive user 
GUI zspec. Final results on redshift success, quality codes, 
completeness, biases, and redshift accuracy are also given 
here. Section 12 investigates the impact of multiple galax- 
ies masquerading as single galaxies in the ground-based peat 
photometry. Section 13 is a guide to DR4 and the two main 
data tables in this paper, which summarize the mask data and 
final redshifts and quality codes. 

Section 14 is provided to enable correlation function mea- 
surements with DR4 data and to aid theorists who wish to 
make mock DEEP2 catalogs from model galaxy data. It de- 
scribes a set of files being distributed along with DR4, which 
describe the selection and redshift success probability for a 
DEEP2 target as a function of location on the sky. Finally, 
Section 15 summarizes the properties of the DEEP2 dataset 
via plots of numbers, colors, magnitudes, etc., as a function 
of redshift. 

Depending on their interests, a reader may wish to focus on 
certain sections of this paper For instance: 

• Users of DEEP2 Data Release 4 data: Key sections 
are §4 (survey parameters), §8 (assessment of biases 
in targeting), §10.3 (presenting sample spectra), §11.3 
and §11.4 (redshift results and completeness), §13 (de- 
scribing the released data tables), §14 (describing the 
released 2-d selection functions) and §15 (describing 
trends in the sample with z). 

• Individuals interested in the design of DEEP2 and 
comparison to other deep surveys: Key sections are 
§2 (design goals), §3 (science results), §4 (survey pa- 
rameters), §5 (comparison to other surveys), §8 (as- 
sessment of biases in targeting), §12 (evaluating how 
frequently spectroscopic targets consisted of multiple 
unresolved galaxies), and §15 (describing trends in the 
sample with z). 

• Users of the DEEP2/DEIMOS data reduction 
pipelines: Key sections are §10 (on the 2-d data re- 
duction pipeline and spectral extractions) and § 1 1 (on 
the redshift determination pipeline). 

• Individuals interested in details of DEEP2 target se- 
lection and techniques: Key sections are §6 (describ- 
ing the target selection criteria), §7 (slitmask design al- 
gorithms), and §9 (observational strategy). 

Throughout this paper, unless specified otherwise we utilize 
magnitudes on the AB system and assume a ACDM concor- 
dance cosmology with Vl,„ = 0.3, ilA=0.7, and Hq = 100 h km 
s"' Mpc~'. The DR4 data tables also use these quantities. 

2. CONSIDERATIONS GUIDING THE DESIGN OF THE DEEP2 
SURVEY 

DEEP2 was originally envisioned as a tool for simultane- 
ously studying galaxy evolution and the growth of large-scale 
structure since z ^ 1 . However, while the survey was being 
designed it was realized that such data could also be used to 
constrain the nature of dark energy (i.e., its equation of state) 
by counting the abundance of groups and clusters as a func- 
tion of velocity dispersion and redshift (Newman et al. 2002; 
Newman & Davis 2000, 2002; Gerke et al. 2005). This was 



one of the first methods proposed for studying dark energy by 
counting massive objects (see also Haiman, Mohr, & Holder 
2001). The velocity and redshift distribution of clusters de- 
pends on both the growth rate of large-scale structure and on 
the volume element (i.e., the amount of volume per redshift 
interval per angular area), which are each related to funda- 
mental cosmological parameters in predictable ways. Count- 
ing the abundance of intermediate-mass groups to high accu- 
racy therefore became a third major goal guiding the design 
ofDEEP2. 

It was clear from the start that one of the most important 
DEEP2 data products would be precision counts, to be com- 
pared with various types of model predictions. Obtaining 
high-precision counts requires: 

• Many galaxies, 

• An accurately known and simple selection function - 
DEEP2 is basically magnitude-limited with nearly uni- 
form sampling on the sky and in redshift space at z > 
0.75 - and, 

• Low sample (commonly referred to as "cosmic") vari- 
ance. 

These aims - large numbers, simple selection function, and 
low cosmic variance - were the central guiding principles for 
DEEP2. They all had to be achieved, however, while using as 
little total telescope time as possible. 

Optimizing the measurement of environment and clustering 
statistics played a major role in setting the sizes and shapes of 
the DEEP2 survey fields. Galaxy properties vary systemati- 
cally with environment, and the length-scale of this influence 
can shed important light on the physical mechanisms causing 
these effects: for example, differentiating between processes 
associated with the parent dark matter halo of a galaxy from 
those that act on larger scales. Since the typical comoving ra- 
dius of galaxy groups and clusters is of order ^ 1 /;"' Mpc, 
measuring the local overdensity of other galaxies around a 
given object (the most common measure of environment used 
today) requires counting neighbors to a separation at least 
this large, corresponding to 1 .5 at z ~ 1. As a result, even 
when corrections for lost area are applied, neighbor counts 
become inaccurate for galaxies too near to survey boundaries 
(cf. Cooper et al. 2005); minimizing the number of lost ob- 
jects in environment studies requires: 

• Fields that are at least ten times wider in each direction 
than the radius within which environments are mea- 
sured, corresponding to a minimum dimension of ^15' 
on the sky for a survey at z ^ 1 . 

This restriction also ensures that fields will be well-sized for 
characterizing groups and clusters, in support of our third sci- 
ence goal. 

The measurement of correlation functions drives the choice 
of a minimum field dimension in similar ways, as the domi- 
nant error terms on large scales increase rapidly for pair sep- 
arations larger than half the minimum field dimension (e.g., 
Bernstein 1994). It is also necessary to survey a large vol- 
ume to minimize the effects of cosmic variance on correla- 
tion function estimates. Typical values of the scale length of 
galaxy clustering, r^, for galaxy populations of interest are 
generally < 5 /z~' co-moving Mpc at z ^ 1 (Coil et al. 2004a, 

2008), corresponding to an angular separation of 7 .4. Mea- 
suring correlation functions well therefore requires: 



• Contiguous volumes spanning several times ro in one 
dimension, with the other two dimensions spanning 
many times ro. 

A field size of 0.5° x 2° was accordingly selected for three out 
of the four DEEP2 fields, corresponding to 4 ro by 16 ro at z = 
1 . Rectangular fields rather than square fields were chosen to 
increase the statistical power of the sample: the opposite ends 
of a strip will be more statistically independent than regions 
closer to each other, reducing cosmic variance (Newman & 
Davis 2002). 

Since environment measures involve the pairwise count- 
ing of objects, the statistical weight of such data increases 
as the square of the volume density of galaxies sampled, 
all other things being equal. ^** Essentially all environmental 
studies in the Poisson-dominated regime should benefit from 
this scaling, including those which employ a wide variety of 
galaxy environment estimators, the measurement of correla- 
tion functions and cross-correlations for dilute samples or at 
small scales, pair counts, comparisons of group and field pop- 
ulations, and studies incorporating the dynamical masses of 
groups. The DEEP2 targeting strategy is therefore designed 
to obtain: 

• A high number density of galaxies (or equivalently, 
number of objects per unit area and redshift) down to 
the faintest feasible magnitude limit. 

However, under the constraint of finite telescope time, the 
need for large, densely sampled fields to enable accurate 
correlation function and environment measurements conflicts 
with the desirability of covering many independent regions 
on the sky to minimize (and allow assessment of) cosmic vari- 
ance. Since its subregions will be correlated due to large-scale 
modes of the power spectrum, the sample variance in a sin- 
gle large field will always be greater than in a larger number 
of widely-separated fields on the sky covering the same total 
area. Our preferences for large field sizes but small cosmic 
variance are most efficiently reconciled by choosing: 

• Several widely separated areas on the sky, each one of 
which exceeding the minimum width and length needed 
for correlation and environment studies. 

The final choice to survey four fields was driven by cosmic 
variance calculations utilizing the QUICKCV code of New- 
man & Davis (2002)^'' and by the instability of standard devi- 
ation measurements with fewer samples. Surveying multiple 
fields spaced well apart in RA also makes telescope schedul- 
ing easier. 

Obtaining a dense sample of galaxies can be achieved by 
selecting down to a very faint magnitude limit; but that would 
require prohibitive observation time. Instead, it is desirable 
to go as faint as one must to obtain enough targets to fill slit- 
masks, but no fainter Since the DEIMOS spectrograph gen- 
erally has 10 slots available for user slitmasks, observation 

^* This is straightforward to understand when considering the use of counts 
within some aperture (cylinder, sphere, etc.) as a measure of local overden- 
sity. When the number density, n, is increased, the fractional uncertainty in 
each object's neighbor count (i.e., its environment estimate) goes down as 
-J^, but the number of objects in the sample also increases as n. Hence, the 
standard deviation of the mean overdensity determined for a given sample 
scales as 1/n, corresponding to an n^ increase in inverse variance/statistical 
weight. 

2' The QUICKCV code is publicly available at 
http://www.phyast.pitt.edu/~janewman/research.html. 



times of order 1 hour per mask are optimal. In order to obtain 
redshifts for a majority of objects in an hour's observation 
time, we desire: 

• A target population brighter than Rab ^ 24. 

In order to include a large enough sample to optimize the fill- 
ing of slitmasks, the final DEEP2 sample extends down to 

i?AB=24.1. 

A single DEIMOS slitmask covers an area on the sky of 
roughly 5 .3 x 16 .7 (cf. §7), of which more than 30% is un- 
suitable for galaxy target slitlets due to vignetting, chip gaps, 
or variation in wavelength coverage. DEIMOS masks can 
generally include ^ 130- 150 objects, corresponding to a sur- 
face density of 2.2-2.6 galaxies per square arcminute, compa- 
rable to the number counts of galaxies down to /?ab ^21.4. 
However, targeting just those objects would result in a rela- 
tively bright and quite dilute sample of galaxies, weighted to- 
wards more luminous galaxies and lower redshifts than were 
of interest for DEEP2. Furthermore, many objects of interest 
would be lost as we are unable to observe multiple objects at 
the same displacement in the spatial direction on a slitmask, 
lest spectra overlap, reducing our sampling of overdense re- 
gions. Given our desire to maximize the number density of 
high-redshift galaxies, we therefore require: 

• Rejection of lower-redshift (z < 0.7) galaxies, and cov- 
erage of each point on the sky by at least two slitmasks. 

In the Extended Groth Strip (EGS, DEEP2 Field 1), however, 
we do not reject low-z galaxies, both to test our selection 
methods and to take advantage of the rich multiwavelength 
data in that field; as a consequence, in that field we require 
coverage of each point on the sky by at least four slitmasks 
to compensate. To keep the number of slitmasks per field 
roughly constant, the width of the area covered is half as large 
in the EGS as for the remaining fields (i.e., the survey area 
goal in this field was 0.25° x 2°, rather than 0.5° x 2° as in the 
other DEEP2 fields). 

The field area of 0.5° x 2° is roughly 60 times the area 
covered by a single DEIMOS slitmask; double-coverage thus 
requires 120 slitmasks per field. Assuming an average of 135 
targets per slitmask, this led to a survey goal of ^65,000 spec- 
tra. 

With the field geometry and sample sizes determined, the 
main parameters remaining were spectral dispersion and spec- 
tral coverage. These also involve tradeoffs. For a fixed detec- 
tor size, low dispersion captures more spectral features and 
thereby enhances spectral information and (potentially) red- 
shift completeness. However, as noted above, high disper- 
sion yields more information per angstrom by resolving inter- 
nal galaxy motions and yielding more accurate redshifts. It 
splits the [O ii] A3727 doublet, converting that feature into a 
unique, robust redshift indicator even when it is the only fea- 
ture visible. At high resolution, OH sky lines are concentrated 
into a few pixels, leaving the remainder of the pixels dark 
and thereby shortening required exposure times (see §10.2). 
High dispersion also reduces flat-fielding errors due to "fring- 
ing" (see Appendix A. 2). In conventional grating spectro- 
graphs, it increases anamorphic demagnification, which nar- 
rows slitwidth images still further and reduces the displace- 
ment of spectra along the spectral direction due to variations 
in slitlet position, maximizing the spectral region common to 
all objects. Finally, high dispersion allows the use of wider 
slits for the same net spectral resolution, thereby capturing 
more galaxy light. 



Evidently, we would like to have both high dispersion and 
broad spectral coverage. The DEIMOS CCD detector/camera 
system was designed to do this by providing a full 8192 pix- 
els along the spectral direction, with minimal chip gap. Our 
ideal choice balancing resolution and coverage would have 
used an 830-line/mm grating, which covers 3900 A of spec- 
trum with DEIMOS. However, this grating suffers signifi- 
cantly from ghosting, which we feared would interfere with 
accurate sky subtraction. We therefore conservatively chose 
the 1200-line/mm grating, which captures 2600 A of spec- 
trum over 2000 separate resolution elements. When centered 
at 7800 A, this spectral range yields at least one strong spec- 
tral feature for essentially all galaxies from z = to z = 1 .4. 
Redshift ambiguities are rare and, if present, can largely be 
resolved in the future by using photometric redshifts (Kirby 
et al. 2007). We therefore selected: 

• DEIMOS' highest-resolution grating, 1200 line mm"', 
typically covering 6500 A to 9100 A and yielding a 
spectral resolution R = A/AA ^ 6000 with a l"-wide 
sHt.3n 

Practical considerations influenced many of the final design 
decisions for DEEP2. The detailed layout of slitmasks on the 
sky was governed by the DEIMOS slitmask geometry (cf. §7). 
The final fields were chosen to have low Galactic reddening 
utilizing the Schlegel, Finkbeiner, & Davis (1998) dust map 
and to be well spaced in RA to permit good-weather observing 
at Keck. Field 1 (the Extended Groth Strip) was specially 
selected for its excellent multiwavelength supplemental data 
(as described in §4). 

Finally, detailed attention was paid to ensuring good sky- 
subtraction accuracy even under the bright OH lines. These 
lines would occupy only a small portion of the spectrum (due 
to the use of high dispersion), but any sky-subtraction errors 
would spread to other wavelengths whenever spectra were 
smoothed. Requirements for excellent sky subtraction include 
(1) a highly accurate and stable calibration of the wavelength 
corresponding to every pixel on the detector, (2) a constant 
PSF for OH lines along each slitlet, and (3) highly repro- 
ducible CCD flat-fielding, which was to be attained by main- 
taining exactly the same wavelength on each pixel between the 
afternoon flat-field calibration and the night-time observation. 
These requirements in turn require a very stable spectrograph, 
uniform image quality over the field of view, careful calibra- 
tion and observing techniques, and high-precision data reduc- 
tion methods. Further details are provided in the Appendix, 
which discusses CCD fringing, sky subtraction, the DEIMOS 
flexure compensation system, image stability, and wavelength 
calibration methods. 

3. SCIENCE HIGHLIGHTS OF DEEP2 

The design resulting from the above considerations has en- 
abled a wide variety of scientific investigations; more than 75 
refereed papers have appeared as of this writing. An early 
summary appeared in the May 1, 2007, special issue of Astro- 
physical Journal Letters (devoted to the AEGIS collaboration, 
described in §4); Table 1 of this paper presents an updated list 
of selected highlights. Amongst our major findings are: 

• DEEP2 confirmed and buttressed the conclusion of Bell 
et al. (2004) that the number of red-and-dead, quenched 

3" Here AA is the FWHM of an OH sky line. 



galaxies has at least doubled since z ^ 1, implying 
that some galaxies arrived on the red sequence rela- 
tively recently (Willmer et al. 2006; Faber et al. 2007). 
Stellar populations imply late quenching (Schiavon et 
al. (2007), and the mean restframe U-V color of galax- 
ies evolves slowly due to these late arivals (Marker et 
al. 2006). 

• Star-formation rates have been measured as a function 
of stellar mass using both 24-/im fluxes and GALEX 
photometry (Noeske et al. 2007a; Salim et al. 2009). 
In star-forming galaxies, the star formation rate follows 
stellar mass closely at all epochs out to z ^ 1, displaying 
a "star-forming main sequence" that can be modeled by 
assuming that larger galaxies start to form stars earlier 
and shut down sooner (Noeske et al. 2007b). The RMS 
scatter about this relation is of order a factor of two. 
Most star-formation at z ^ 1 is taking place in normal 
galaxies, not in merger-triggered starbursts (Noeske et 
al. 2007b, Marker 2008). Red spheroidal galaxies are 
truly quenched, and the strong [O ii] emission in some 
is due to AGN and/or LINER activity (Konidaris 2008, 
Yan et al. 2006). 

• The fraction of quenched galaxies is larger in high- 
density environments (Cooper et al. 2006) and is higher 
in groups versus the field (Gerke et al. 2007). How- 
ever, the decline in global star formation rate since z ^ 1 
is strong in all environments, and cannot be driven by 
environmental effects (Cooper et al. 2008). The frac- 
tion of quenched galaxies in high-density environments 
decreases with lookback time and essentially vanishes 
at z > 1.3, implying that massive galaxies started shut- 
ting down in bulk near z > 2. Environment correlates 
much more strongly with galaxy color than luminos- 
ity; a color-density relation persists to z > 1 (Cooper 
et al. 2006, 2007, 2010). Many massive star-forming 
galaxies were still found in high-density regions at z ^ 1 
but have no analogs today; they must have quenched 
in the interim time (Cooper et al. 2006, 2008). Post- 
starburst galaxies are found in similar environments as 
quenched galaxies at z 0.8, but similar to star- forming 
galaxies at z 0.1; given the overall movement of galax- 
ies towards denser environments at lower z, the absolute 
(as opposed to relative) overdensity around quenching 
galaxies may have not have evolved with redshift (Yan 
et al. 2009). Finally, massive early-type galaxies in 
high-density regions are >25% larger in size relative 
to their counterparts in low-density environs at z ^ 1, 
strongly suggesting that minor mergers may be impor- 
tant in the size evolution of massive ellipticals (Cooper 
et al. 2012a). 

• The roots of the Tully-Fisher relation are visible out 
to z ^ 1 (Kassin et al. 2007). Disturbed and merg- 
ing galaxies have low rotation velocities but high line- 
widths that will place them on the TF relation after set- 
tling (Covington et al. 2009). The zero point of the 
stellar-mass TF relation has evolved very little since 
z ^ 1, in agreement with LCDM based models (Dut- 
ton et al. 2011), and the Faber-Jackson relation for 
spheroids likely arose from the TF relation via merg- 
ers and quenching. 

• The number of disky galaxies has declined and the 
number of bulge-dominated galaxies has increased 



since z ^ 1 (Lotz et al. 2008). However, X-ray de- 
tected AGN are not preferentially found in major merg- 
ers, and the fraction of merging galaxies is not appre- 
ciably higher at z ^ 1 than now (Pierce et al. 2007; Lin 
et al. 2004, 2008; Lotz et al. 2008). Many distant galax- 
ies previously classed as peculiar or merging are more 
likely normal disk galaxies in the process of settling. 

• The clustering properties of a wide variety of galaxy 
and AGN samples have been measured. The auto- 
correlation function of typical DEEP2 galaxies implies 
halo masses near 10'^ M0. The clustering amplitude 
is a stronger function of color (Coil et al. 2008) than 
luminosity (Coil et al. 2006b), matching the results 
from environment studies. Red galaxies cluster more 
strongly than blue ones, with intermediate-color "green 
valley" galaxies preferentially found on the outskirts of 
the same massive halos that host red galaxies (Coil et 
al. 2008). The halo mass associated with star formation 
quenching has been quantified both through the corre- 
lation function (Coil et al. 2008) and through the mo- 
tions of satellite galaxies (Conroy et al. 2007). Coil 
et al. (2006a) use the clustering of galaxies with group 
centers to separate out the "one-halo" and "two-halo" 
contributions to the correlation function. Using the 
cross-correlation of AGN and DEEP2 galaxies. Coil 
et al. (2007) find that quasars cluster like blue, star- 
forming galaxies at z = 1, while lower-accretion X-ray 
detected AGN cluster similarly to red, quiescent galax- 
ies and are more likely to reside in galaxy groups (Coil 
et al. 2009). Comparison between X-ray and optically- 
selected AGNs indicates that the fractions of obscured 
AGN and Compton-thick AGNs at z ^ 0.6 are at least 
as large as those fractions in the local universe (Yan et 
al. 2011). 

• Ubiquitous outflowing winds have been detected in nor- 
mal star-forming galaxies for the first time, at z ^ 1.3 
using [Mg ii] absorption (Weiner et al. 2009). Out- 
flow speeds are higher in more massive galaxies and 
are comparable to the escape velocity at all masses. At 
z=0. 1-0.5, outflows are detected in both blue and red 
galaxies but are more frequent in galaxies with high IR 
luminosity or recently truncated star formation (Sato 
et al. 2009). Overall, these results suggest that stel- 
lar feedback is an important process in regulating the 
gas content and star formation rates of star-forming 
galaxies and that galactic-scale outflows play an impor- 
tant role in the quenching and migration of blue-cloud 
galaxies to the red sequence. 

• The cause of AGN activity and the role of AGN in 
galaxy quenching remain major puzzles. X-ray AGN 
are found in massive galaxies at the top of the blue 
cloud, on the red sequence, and in the green valley 
(Nandra et al. 2007), as expected if AGN activity is 
associated with black-hole growth and bulge-building. 
However, X-ray AGNs are not especially associated 
with mergers (Pierce et al. 2007) or with high star for- 
mation rates (Georgakakis et al. 2008). A surprise is 
the large number of luminous, hard, obscured X-ray 
sources in morphologically normal red spheroidals and 
post-starburst galaxes (Georgakakis et al. 2008), well 
after star formation should nominally have stopped. If 



AGN feedback causes galaxy quenching, these persis- 
tent hard AGN imply that considerable nuclear gas and 
dust must remain in the central regions, where it con- 
tinues to feed the black hole. 

To summarize these findings: the roots of most present-day 
galaxy properties are visible far back in time to the edge of 
DEEP2 (z = 1 .4, more than 8 billion years in the past); galaxy 
evolution over this time has on balance been a fairly regular 
and predictable process; and a good predictor of a galaxy's 
properties at any era is its stellar mass. An overarching hy- 
pothesis that may unite all of these findings is that the most 
important physical driver controlling a galaxy's evolution at 
z < 1.4 is the growth of its dark halo mass versus time. On 
the other hand, several of the above findings point to close 
couplings between AGNs, starbursts, quenching, and galactic 
winds, and the causal connections among these seems more 
complex than the simple halo-driven model would predict. 
These issues are the great frontier for future work. 

The present data in DR4 do not exhaust the information 
from DEEP2. Many additional quantities have been derived 
from survey data, and have either been published separately 
(see Table 1, or may be obtained by contacting the indi- 
vidual scientists involved. These include measurements of 
emission-line equivalent widths (from B. Weiner, R. Yan, 
N. Konidaris, and J. Harker); improved velocity widths for 
emission and absorption lines (from B. Weiner and S. Kassin); 
local environmental densities (from M. Cooper); group iden- 
tifications and group memberships (from B. Gerke); stellar 
masses (from K. Bundy and C. N. A. Willmer); ISM metal- 
licities (from A. Phillips); D4000 and Balmer-line absorption 
strengths (from J. Harker and R. Yan), optical AGN IDs (from 
R. Yan), morphologies (from J. Lotz), pair catalogs (from 
L. Lin), satellite kinematics (from C. Conroy), star-formation 
rates (from K. Noeske and S. Salim), and indicators for post- 
starbust galaxies (from R. Yan). 

4. SURVEY OVERVIEW AND FIELD SELECTION 

In Section 2, we presented the goals which were used to fix 
the basic properties of the DEEP2 survey. In the remainder of 
this paper, we describe the design and execution of the project 
and the resulting datasets in more detail. This first section 
serves as a basic introduction to the DEEP2 survey. 

To first order, DEEP2 is a magnitude-limited spectroscopic 
galaxy redshift survey with limiting magnitude of Rab = 
24.1. Instrument and exposure parameters are summarized 
in Table 2. The survey was designed to be executed in four 
separate 0.5° (or 0.25° in the case of the Extended Groth 
Strip) X 2° rectangular fields widely spaced in Right Ascen- 
sion, with an average of 120 DEIMOS slitmasks per field. 
The spectral setup used the 1200-line/mm high-resolution 
DEIMOS grating with a spectral resolution of R^ 6000 and a 
central wavelength of 7800 A; the typical exposure time was 
1 hour per mask. The survey is primarily sensitive to galax- 
ies below a redshift of z ^ 1 .45, past which the [O ii] 3727 
A doublet moves beyond the red limit of our typical spectral 
coverage. The total number of spectra obtained is 52,989, and 
the total number of objects with secure (classes with > 95% 

" In some cases, spectra of fainter objects which serendipitously fell 
on DEIMOS slitlets were obtained and extracted by our data reduction 
pipelines; such objects are not presented here, as their selection is less well- 
characterized than the main DEEP2 sample, though they may be included in 
a later data release. 



repeatability) redshifts is 38,348. Due to the total allocation of 
telescope time being somewhat smaller than originally envi- 
sioned, the total area of sky covered by the final survey sample 
is roughly 3.0 square degrees (rather than 3.5), and 41 1 (rather 
than 480) slitmasks were observed. 

Objects are pre-selected in DEEP2 Fields 2, 3, and 4 us- 
ing broad-band CFHT 12K BRI photometry to remove fore- 
ground galaxies below z ^ 0.7 (cf Section 6.3); this selection 
is essentially complete for galaxies at z > 0.75. This preselec- 
tion removes 60% of galaxies brighter than Rab = 24.1, thus 
multiplying the efficiency of the survey for studying galaxies 
at z ^ 1 by nearly 2.5 x. Field 1 is the Extended Groth Strip 
(EGS, Groth et al. 1994) and is treated differently owing to the 
wealth of ancillary data there. In that field, we do not reject 
low-redshift galaxies on the basis of their colors. 

To select objects, we first define a target pool of candi- 
date galaxies in each field, consisting of all objects between 
^AB = 18.5 and R^b = 24.1 except those with exceptionally 
low surface brightness or large (> 80%) Bayesian probability 
of being a star. The photometric catalogs ("pcaf's) used for 
this selection are described in §6.1. We then apply a color 
pre-selection to this pool to remove z < 0.7 objects (in Fields 

2, 3, and 4). The remaining galaxies are used to design the 
masks; in a typical field, roughly 60% of potential targets re- 
ceive slitlets. 

The chosen depth of Rab = 24. 1 is suitable for several rea- 
sons. At this depth, DEEP2 probes down to L* in the lu- 
minosity function at z ^ 1 .2. The surface density of 20,000 
objects/n° (after high-z preselection) is large enough to pack 
slitmasks efficiently while leaving a modest excess for flexi- 
bility. Finally, the majority of galaxies at this magnitude limit 
still yield secure redshifts in one hour of observation time. 

Information on the four DEEP2 fields is given in Table 3, in- 
cluding coordinates, average reddening, field sizes, number of 
slitmasks planned and observed, and number of target galax- 
ies and redshifts. By ranging in RA from 14'' 20" to 2*' 30™, 
with spacings of 3-7 hours between fields, the four DEEP2 
fields are well suited for observing in what is historically a 
good-weather period at Keck, ranging from March through 
October Masks were designed over 0.5° wide by 2° long re- 
gions in Fields 2, 3, and 4, but not all of these masks were ul- 
timately observed. As noted. Field 1 (EGS) is special: it hosts 
a wide variety of deep multiwavelength data from X-rays to 
radio. It is also one of the darkest and most dust-free regions 
of the sky. Since no z > 0.7 preselection was done in Field 
1, the surface density of targets on the sky there is roughly 
twice that in the other DEEP2 fields; in order to sample high- 
redshift objects at the same number density as in other fields, 
only half as much area could be covered using the same num- 
ber of slitmasks. As a result, the spectroscopy in EGS covers 
an area that is only ^16' wide (but still 2° long). All masks in 
EGS were observed. 

The field layouts and final slitmask coverage in Fields 2, 

3, and 4 are shown in Figure 1. The boundary of each in- 
dividual CFHT 12k BRI peat photometric pointing is shown 
by dashed lines; individual fields are labelled by their field 
and pointing number (e.g.. Pointing 12 is the second CFHT 
pointing in Field 1). The chevron pattern permits slitlets to 
be aligned with atmospheric dispersion both east and west of 
the meridian (the odd pattern in Field 2 was a result of this 

'^ Three sets of weak deweighting functions are applied to suppress very 
faint galaxies, very nearby galaxies, and very blue galaxies, as explained in 
§6.2. 
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requirement, combined with the high priority assigned to ob- 
servations in Field 1); see Section 7 for details. The grey 
scale at each point represents the probability that a galaxy in 
that mask meeting the DEEP2 sample selection criteria was 
targeted for spectroscopy and that a secure (Q = 3 or Q = 4, cf. 
§11) redshift was measured. This spatial selection function 
can be downloaded by users who wish to perform cluster- 
ing measurements or make theoretical models of the survey 
(see §14). Certain areas were redesigned following the first 
semester of survey operations, which caused a few regions to 
be observed twice with different masks (darker stripes). All 
DEEP2 masks, including these regions of double coverage, 
are included in DR4. The density of masks is especially high 
where CFHT 12K pointings overlap in Field 4. This dupli- 
cation of coverage occurred because slitmasks were designed 
separately pointing by pointing in Fields 2^ (but not EGS). 
Although photometry was obtained for three CFHT pointings 
in each of fields 2-4, due to limited time for spectroscopic ob- 
servations field 23 (which had inferior photometry) and most 
of field 43 were not observed with DEIMOS. The photomet- 
ric (peat) catalogs for these fields are still included in the data 
release accompanying this paper. 

A similar plot of combined selection and redshift success 
probability is given for Field 1 (the Extended Groth Strip) in 
Figure 2. The layout of this field is long and narrow both 
to allow dense coverage and in order to follow the geome- 
try of preexisting space-based data. The strip is divided into 
eight blocks, each containing 15 masks. Within each block, 
eight masks run perpendicular to the strip and seven run par- 
allel to it. Unlike Fields 2, 3, and 4, a single peat was cre- 
ated by merging the photometry from all four pointings be- 
fore the masks were defined. There are thus no regions where 
separately-designed masks overlap, unlike in Field 4. 

A summary of the existing data in all DEEP2 fields is given 
in Table 4, and skymaps of the major multiwavelength EGS 
surveys are shown in Figures 3 and 3. The data in EGS are 
particularly deep and rich. EGS is one of only two deep-wide 
Spitzer MIPS and IRAC regions (Fazio et al. 2004, Papovich 
et al. 2004, Dickinson 2007 [FIDEL survey]) and is receiving 
ultra-deep imaging at 3.6 fim and 4.8 /im in the SEDS survey, 
a Warm Spitzer mission. It has one of only two large two- 
color HST ACS mosaics (Davis et al. 2007, Lotz et al. 2008), 
the other being GEMS in ECDFS (Rix et al. 2004). It has 
the deepest GALEX imaging on the sky and the deepest wide- 
area Chandra mosaic (800 ksec. Laird et al. 2008). It is a 
CFHT Legacy Deep Survey field, is slated for deep Herschel 
imaging, and many other deep optical, IR, sub-mm, and radio 
data have been taken or planned. This includes the follow-on 
DEEP3 survey (Cooper et al. 2011, 2012a), which is acquir- 
ing ~8,000 new spectra down to /?ab = 25.5 and will dou- 
ble the redshift sampling density over the area covered by the 
HST ACS mosaic (see Table 4). Most recently, EGS has be- 
come one of five fields targeted by the CANDELS (Cosmic 
Assembly Near-infrared Deep Extragalactic Legacy Survey) 
Multi-Cycle Treasury program on HST (Grogin et al. 2011; 
Koekemoer et al. 2011); it is being imaged deeply with both 
the ACS and WFC3 instruments as part of this project. 

To exploit the great richness of data in EGS, a broader re- 
search collaboration has been formed called AEGIS (the All- 
wavelength Extended Groth Strip International Survey). This 
collaboration is combining efforts from more than a dozen 
teams who have obtained data in this field ranging from X- 
ray to radio wavelengths. More information on AEGIS may 
be found at the AEGIS website (http://aegis.ucolick.org) and 
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Figure 1. 2-d redshift completeness maps, w(a, 5), in DEEP2 Fields 2, 3, 
and 4, updated from Cooper et al. (2006). The boundaries for each point- 
ing of the DEEP2 CFHT 12K BRI photometry are indicated by the dashed 
lines, labelled by the pointing number (e.g., pointing 21 is the first pointing 
in DEEP2 Field 2). The grey scale at each point represents the probability 
that a galaxy in that mask meeting the DEEP2 sample selection criteria was 
targeted for spectroscopy and that a2 = 3or2 = 4 redshift was measured. 
Most of the slitmasks for each pointing are in two rows of approximately 20 
masks each. The small white lines correspond to gaps between the DEIMOS 
CCDs. The darker regions show areas where masks overlap and objects are 
therefore observed with higher probability. At the intersection of the top and 
bottom rows, objects may be observed twice, allowing verification of redshift 
repeatability, etc.; however, the vertical masks filling in the "fishtails" at the 
east end of each pointing do not contain duplicate objects. 



in the May 1, 2007, special issue of Astrophysieal Journal 
Letters. 

With Field 1 designated to be the Extended Groth Strip, 
Fields 2, 3, and 4 were selected by finding regions free of 
bright stars, with low reddenings based on the IRAS dust map 
(Schlegel, Finkbeiner & Davis 1998), spaced roughly 4 hours 
apart in RA, and weighted toward higher RA where Keck 
weather is better. Their distribution also avoids the prime 
Galactic observing season and balances the number of nights 
required in each semester, in order to ease telescope schedul- 
ing. The chosen fields all have declinations such that they 
are observable for more than 6 hours with an airmass less 
than 1.5, to ensure that differential refraction across a slit- 
mask is below .2 (DEIMOS does not have an atmospheric 
dispersion compensator). Fields 3 and 4 are also in the multi- 
ply observed Equatorial Strip (Stripe 82) of the Sloan Digital 
Sky Survey (SDSS, York et al. 2000, Adelman-McCarthy et 
al. 2007, Ivezic et al. 2007), and are visible from both northern 
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Figure 2. As Figure 1, but for tlie Extended Grotli Snip (DEEP2 Field 1). 
Pointing boundaries for the CFHT 12K BRI photometry catalog are again 
shown as dashed lines. Masks were designed in eight blocks along the long 
direction of the Strip. Each block has 15 masks, eight perpendicular to the 
strip and 7 parallel. Masks overlap extensively but unlike in Fields 2, 3, and 
4 there are very few duplicate observations built into the design. DEEP2 
pointing 14 is omitted here, as mask design in that region followed different 
algorithms to account for its poorer photometry. Pointing 14 is therefore not 
included in any DEEP2 large-scale structure analyses. 
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Figure 3. Initial data in the Extended Groth Strip (DEEP2 Field 1). De- 
tails of the individual data sets are given in Table 4. The CFHT 12K imag- 
ing provides the BRI peat photometry used for DEEP2 target selection. The 
HST/ ACS mosaic is one of the largest two-color V +1 mosaics on the sky. 
The Guaranteed Time Observation Spitzer/IRAC data and the Chandra 200 
ksec data are the deepest/widest of their kind, and the CFHT Legacy Survey 
provides valuable synoptic variability data. EGS has also been imaged deeply 
with the VLA at 20 cm and with the GMRT (not shown) at 50 cm. 
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Figure 4. Data in the Extended Groth Strip are continually expanding. De- 
tails of the individual data sets are given in Table 4. The FIDEL SpitzerlMlVS, 
AEGIS-X Chandra/ACIS, and GALfiX/FUV-l-NUV imaging are the deepest 
exposures of their size on the sky. The ongoing DEEP3 survey (Cooper et 
al. 2012, in prep) will increase the number of DEEP2 spectra by 50% and 
triple the weight of environmental data in the upper three fields of EGS. The 
Warm Spitzer Deep Extragalactic Survey (SEDS) will provide a total of 10 
hours of integration time per pointing in IRAC Channels 1 and 2 (centered 
at 3.6 and 4.5 ^m). The NEWFIRM survey will provide JK photometry to 
~24 AB mag, while the NEWFIRM Intermediate Band Survey (PI: Pieter 
van Dokkum; cf. van Dokkum et al. 2009) will measure photoz's to roughly 
Xab ~ 23.4. In addition (not shown), EGS is a deep SCUBA-2 Legacy field 
and a deep field for the Herschel HERMES survey with the PACS and SPIRE 
instruments. 

and southern hemispheres. The recently-released deep SDSS 
photometry in this stripe should yield useful photometric red- 
shifts down to DEEP2-like depths. 

The original design of the DEEP2 survey consisted of 480 
slitmasks. The expected number of 135 slitlets per mask 
would then target ^65,000 galaxies, of which 52,000 were ex- 
pected to yield secure redshifts assuming an 80% success rate. 
One can therefore crudely think of DEEP2 as having roughly 
twice the number of objects and three times the volume as 
the Las Campanas Redshift survey (Shectman et al. 1996), 
but at z ~ 1 rather than at z ^ 0.1. The design volume of 
9 X lOf'hr^ Mpc^ would be expected to contain only a hand- 
ful of rich clusters but is large enough to count both galaxies 
and groups of galaxies to good accuracy. At one hour per 
mask and eight masks per night, such a survey would take 60 
clear nights of Keck time. Ninety nights were allocated via 
a Keck LMAP (Large Multi-Year Approved Project) proposal 
with Marc Davis as PI, which left a cushion of 50% for bad 
weather and equipment malfunction. 

The above plan was followed closely, but the yield of re- 
liable redshifts per galaxy targeted was only 71%, not 80%, 
owing to the fact that more galaxies than expected were be- 
yond the effective redshift limit of z ^ 1.4 (§11). Weather 
also did not fully cooperate, with the result that only 411 out 
of the 480 original slitmasks were observed, yielding a total 
of ^53,000 spectra. The final area and numbers of slitmasks 
and targets observed in each field are given in Table 3, and 
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sky maps of the regions covered are shown in Figure 1 and 
Figure 2. The final survey observed 86% percent of the pro- 
posed number of slitmasks and 88% of the proposed number 
of galaxies. 

5 . DEEP2 IN COMPARISON TO OTHER Z ~ 1 SURVEYS 

In this section and Table 6, we compare the properties of 
DEEP2 (considering both Fields 2-4 and EGS separately) 
to other large z ^ 1 surveys to date in a variety of ways. 
The other projects considered include the Team Keck Red- 
shift Survey in GOODS-North (TKRS, Wirth et al. 2004), the 
VVDS-deep survey in CDFS and in one other field (Le Fevre 
et al. 2005), the VVDS-wide survey in four fields (Garilli et 
al. 2008), zCOSMOS-bright in the COSMOS field (Lilly et 
al. 2007), and PRIMUS, which covers 7 fields including sev- 
eral of the above (Coil et al. 2011).^^ The first five are conven- 
tional spectroscopic surveys, but PRIMUS is a low-resolution 
prism survey yielding less precise redshifts but targeting many 
more galaxies. DEEP2 and TKRS have released all data, but 
the properties of the others have been deduced from stated 
plans and/or partial published data. Details of the data and 
sources used for these comparisons are given in the footnotes 
to Table 6. 

Given the significant redshift failure rate in all distant 
galaxy surveys, it is important to have a consistent definition 
of "redshift success." "Reliable", "secure", or "robust" red- 
shifts are defined in this paper to be those with >95% proba- 
bility of being correct, according to their authors (for all sur- 
veys described, this means quality codes 3 and 4 [or equiv- 
alent, e.g. including codes 13, 14, 23, and 24 for VVDS], 
but not 1 or 2.) In comparing surveys we use the sizes of 
the presently published datasets for all surveys except zCOS- 
MOS, for which we use the design values. 

We now proceed to compare the power of the various z ^ 
1 surveys using a variety of metrics. Numerical results are 
summarize in Table 6. 

• Redshift efficency. A useful concept is "redshift effi- 
ciency," for which we offer two definitions. The first is the 
overall redshift efficiency, Zetfici- which is defined as is the 
fraction of all slitlets (including stars) that eventually yield 
reliable galaxy (or QSO) redshifts. Values of Zeffici are shown 
for the various surveys in Table 6. Apart from the two VVDS 
surveys, there is considerable degree of homogeneity in this 
quantity, with efficiencies hovering between 56-72%. The ef- 
ficiencies of VVDS-deep and VVDS-wide ai-e only 28-40%, 
however, partly because these surveys do not exclude stars, 
which hence take up a larger fraction of the slitlets than in 
other surveys, but mostly because of our high standard for 
redshift reliability, >95% repeatability, which excludes many 
VVDS redshifts. 

The second definition of efficiency is the galaxy-only effi- 
ciency, Zeffici, which is defined as the fraction of galaxy-only 
slitlets that yield reliable galaxy redshifts.^'* The percentages 
rise by 3-4%, to 59-76%, for the non-VVDS surveys and are 
now 42-43% for the two VVDS surveys. The efficiency of 
VVDS-wide rises the most due to the high fraction of stars 
amongst its targets. 

^^ Although they include some objects at comparable redshifts, the AGES 
(Hickox et al. 2009; Kochanek et al. 201 1), BOSS (Aihara et al. 201 1), and 
WiggleZ (Drinkwater et al. 2010) redshift surveys are not included since their 
magnitude limits are 3-4 magnitudes brighter than these other surveys . 

^* To calculate this number, we assume that the fraction of objects with 
uncertain redshifts that are stars is the same as their fraction amongst objects 
with secure redshifts/identifications. 
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Figure 5. Number of redshifts vs. volume sampled for major spectroscopic 
redshift surveys out to z ~ 1 . For all surveys except zCOSMOS-bright, we 
use the number of redshifts and areas covered to date (as given in Table 6 for 
the distant surveys). For zCOSMOS-bright, which is still in progress, design 
data are used. Numbers in all cases use reliable redshifts only (those with 
probabilities of being correct >95%). Volumes covered are computed from 
the field areas and magnitude limits of each survey as described in §4. The 
TKRS survey is omitted from this figure (although included in Table 6) due 
to the relatively small sample size and volume. PRIMUS has low-resolution 
spectra, and provides coarser redshift information for brighter samples than 
the other higher-z surveys shown. If it is set aside, DEEP2 leads amongst 
distant surveys in both volume surveyed and number of reliable redshifts. 

• Volume sampled vs. number of objects. Figure 5 com- 
pares the number of reliable redshifts versus volume sampled 
for both high- and low-redshift surveys. The volume used 
is computed from the areas in Table 6, bounded by the red- 
shifts corresponding to the 2.5 percentile and 97.5 percentile 
point within a given survey (or by z = 0.2 and z= 1.2, the red- 
shift range the team uses for science, in the case of PRIMUS). 
For surveys other than DEEP2 and TKRS (where the actual 
percentiles are used), these redshift limits are computed from 
analytic redshift distributions of the form n{z) oc z^e~^°, us- 
ing fits to zo as a function of limiting magnitude determined 
as described in Coil et al. (2004b), but using the full DEEP2 
dataset for calibration. This results in values of zq of 0.26, for 
DEEP2 and TKRS, 0.224 for VVDS-wide and zCOSMOS- 
bright, and 0.283 for VVDS-deep. Among the contemporane- 
ous distant spectroscopic surveys (i.e., excluding PRIMUS), 
DEEP2 has the best combination of total number and volume 
sampled, with 2.5 times as many redshifts as the next com- 
petitor, zCOSMOS-bright, and roughly twice its effective vol- 
ume. PRIMUS surpasses DEEP2 and all other high redshift 
surveys in this space, targeting a larger number of galaxies 
over a larger volume (but with lower spectral resolution). 

• Co-moving number density. Figure 6 uses the same 
model to compute cosmology-independent sample densities, 
expressed as the number of galaxies with secure redshifts per 
square degree per unit redshift interval; we also include the 
corresponding curve for the Sloan Digital Sky Survey (for 
which we take n{z) oc ^^gC-z/o.ovs) .5^ ^^^ comparison. Figure 7 
converts this to the co-moving number density of objects, as- 
suming our standard cosmology. The diamonds mark the me- 
dian redshifts of the various surveys. The sample density is 
important since the statistical weight of a survey for environ- 
mental and small-scale clustering measurements increases as 
the square of this density if sky area is held constant (see be- 
low). Figure 7 shows that DEEP2/EGS is ~' 4 times denser 
than zCOSMOS-bright and VVDS-deep in the redshift range 
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Figure 6. Sample number densities of various spectroscopic surveys, ex- 
pressed as the number of galaxies per square degree per unit redshift inter- 
val (eliminating any dependence on cosmological assumptions). Curves are 
based on the stated magnitude limits, the number of reliable galaxy redshifts 
per square degree (>95%, see Table 6), and the redshift distribution model 
described in §4, applied to each survey. Note that the statistical weight of 
a survey for environmental and clustering purposes increases as the square 
of the co-moving number density if area and magnitude limit are held fixed. 
Because of its comparably high number density at z ^ 1 , DEEP2 is signif- 
icantly better suited for environmental studies at intermediate redshifts than 
other deep surveys. The PRIMUS survey is not shown here, as in that sample 
spectral resolution, rather than sample number density, limits the ability to 
measure Mpc-scale environments for individual objects. Its number density 
at peak (z ~ 0.4-0.6) is roughly 40% that of DEEP2, or approximately 50% 
larger than that of zCOSMOS-bright or WDS-deep. 

z = 0.6-1.0 and -- 20 times denser than WDS-wide. Both 
the higher secure-redshift rate and denser targeting of DEEP2 
contribute to this difference. We omit PRIMUS from these 
figures and the environment figure of merit calculation (be- 
low) as the ability to determine local environment for individ- 
ual objects in the PRIMUS sample is limited by redshift er- 
rors ((T;/(l +z) ^ 0.005 ^ 15 /i~' Mpc comoving) rather than 
sample number density; its number density is intermediate be- 
tween that of zCOSMOS-bright or WDS-deep and that of 
DEEP2. 

• Statistical weight for environment measures. The figure 
of merit that describes a survey's total statistical weight for 
environmental purposes scales as N^/A, where A^ is the num- 
ber of reliable redshifts and A is the area covered (assuming 
matching redshift distributions).^'' This metric is given for the 
various surveys in Table 6 ("Env metric"). The combined en- 
vironmental weight of DEEP2H-EGS will be 4 times that of 
zCOSMOS-bright when that survey is fully completed, and 
is 14-35 times the weight of WDS-deep and WDS-wide in 
their present states. 

• Spectral power Other measures of survey power include 
spectral resolution and the total number of independent spec- 
tral resolution elements sampled. With a spectral resolution 
that is 10 and 26 times higher than zCOSMOS-bright and 
WDS respectively, DEEP2 is the only distant survey that can 
measure the internal kinematics of galaxies, satellite motions, 

^^ On scales where Poisson variance, rather than cosmic variance, dom- 
inates, the statistical weight of a survey for environment measures scales is 
proportional to the number of pairs of galaxies it contains on the relevant 
scale. For a fixed redshift distribution, the number of pairs will be propor- 
tional to the number of galaxies with secure redshifts, N, times the number 
of companions per galaxy. This latter is proportional to the number density 
of galaxies, and hence (for constant z. distribution) the surface density N/A. 
Hence, the total number of pairs scales as N^/A. 



Figure 7. Number densities of various spectroscopic surveys, expressed as 
number of galaxies per comoving Mpc"^^ (assuming our standard ACDM cos- 
mology). Curves are based on the same models and data used in Figure 6. The 
tick marks indicate the estimated median redshifts of each survey; the open 
tick mark is for EGS. Since the weight of a survey for environmental and 
clustering data scales as the square of the number density, the square of the 
relative heights of the tick marks ranks the surveys according to the strength 
of their clustering information per unit volume at their peak redshifts. DEEP2 
and its EGS subset are the most densely sampled surveys shown, followed by 
SDSS, zCOSMOS-bright, WDS-deep. and WDS-wide. The number den- 
sities of all distant surveys considered decline roughly exponentially with 
redshift. As in Figure 6, the PRIMUS sample is not shown here as its number 
density is not the limiting factor for small-scale environment measures; it is 
intermediate between the DEEP2 and zCOSMOS-bright surveys. 

and small group velocity dispersions. Each DEEP2/DE1MOS 
spectrum has ^ 2000 independent spectral resolution ele- 
ments, versus - 320 for zCOSMOS-bright, -120 for WDS- 
wide and WDS-deep, and - 20 for PRIMUS. Multiplying 
by the number of design targets (excluding stars) gives the 
total number of spectral elements present in galaxy spectra. 
By this metric, DEEP2H-EGS has 16 times more information 
than zCOSMOS-bright, 10 and 33 times more information 
than WDS-wide and WDS-deep, and 185 times more in- 
formation than PRIMUS. These figures are generous, as we 
have here assumed design numbers for total survey size and 
that all spectra in all surveys are equally likely to yield reliable 
redshifts once stars are excluded. 

• Galaxy counts. A final metric of survey power is the ac- 
curacy of galaxy counts. Two sources of noise contribute, 
Poisson statistics in the number of galaxies counted and sam- 
ple (or "cosmic" variance), which is determined by the num- 
ber of fields, their areas, and geometries. Using the publicly- 
available QUICKCV code of Newman & Davis (2002), we 
find that the cosmic variance in the count of a population 
of galaxies with ro = 4 /i"' Mpc (comoving) and 7 = 1.8 at 
z = 0.75-0.95 will be 7.6%, as opposed to 30% for TKRS, 
14% for the completed area of WDS-deep, 6.6% for WDS- 
wide, 12% for the full zCOSMOS-bright area, and 5.1% for 
PRIMUS (here we assume square field geometry and equal 
area per field for all surveys but DEEP2). As a result, the noise 
in DEEP2 galaxy counts (e.g. measurements of the abundance 
of any particular population) will be smaller than in other dis- 
tant z — 1 grating spectroscopic surveys (though not than in 
PRIMUS); the larger sample size yields smaller Poisson er- 
rors, while the relatively long, narrow field geometry utilized 
gives smaller cosmic variance than nearly-square fields of the 
same area do (Newman & Davis 2002). 

To conclude, by all of these measures of survey power, 
DEEP2H-EGS is more powerful than any contemporaneous. 
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high-resolution spectroscopic survey at z ^ 1; in many cases, 
by a factor of 4-30 x. This power enables a wide variety of 
unique science. 

We have excluded PRIMUS in many of these comparisons; 
it has largely followed after DEEP2 in time, and operates in 
a very different domain. Due to the very low resolution spec- 
troscopy employed, redshift uncertainties are relatively large, 
^ 1800-3300 km s"' (depending on redshift, due to the wave- 
length dependence of resolution in prism spectroscopy). The 
large redshift errors make it unsuitable for measuring envi- 
ronment on ^Mpc scales for individual objects (see Cooper 
et al. 2005), which is why we have not computed an environ- 
mental figure of merit for it in Table 6. However, with nearly 
80,000 redshifts over 7 different fields, it has, by a significant 
factor, the lowest cosmic variance of any of these surveys and 
will be able to count rare objects and measure mass and lumi- 
nosity functions for bright objects with unparalleled accuracy. 
It will also be effective at measuring clustering extending to 
relatively large scales; this allows measurement of the average 
overdensity of a population, an important measure of environ- 
ment complementary to the measurement of overdensities for 
individual objects. However, information on redshift-space 
distortions will be lost (compared to these other surveys) due 
to the larger redshift errors. 

6. OBJECT SELECTION 

6.1. Photometric Catalog and Pgai Probabilities 

The photometric catalogs used for selecting DEEP2 targets 
(the peat catalogs, described in detail in Coil et al. 2004b) 
are derived from Canada-France-Hawaii Telescope (CFHT) 
images taken with the 12K x 8K mosaic camera (Cuillandre 
et al. 2001) in the B, R, and / bands. These images were taken 
as part of a major weak-lensing survey by Nick Kaiser and 
Gerry Luppino, and we were very fortunate to piggy-back on 
their efforts. The field of view of the camera is 28' x 42', as 
illustrated by the pointing boundaries in Figures 1 and 2. 

The /?-band images have the highest signal-to-noise and 
are used to define the object catalog, which is complete to 
Rab > 24.15 in all pointings. Objects were identified us- 
ing the imcat software package written by N. Kaiser and de- 
scribed in Kaiser, Squires & Broadhurst (1995). This package 
was also used to calculate other image parameters used in ob- 
ject selection, such as object sizes and photometry. In Fields 
2, 3, and 4, the calibrated peat catalogs for separate point- 
ings are only considerated separately, and slitmasks in these 
fields were designed independently (though pointings overlap 
in Field 4, which means that one object can appear in two 
peat's). In Field 1 (EGS), the peat's were merged to make a 
single catalog for the whole area before the masks were de- 
signed (retaining the higher-quality photometry in overlap re- 
gions), such that each galaxy can appear only once. Each ob- 
ject is assigned an 8-digit object number; the first digit of the 
object number indicates the field the object was found in, the 
second digit indicates the CFHT pointing number, and the re- 
maining digits provide a unique identifier, counting upwards 
from zero. As an example, object 32001226 is the 1,227th 
object in the imeat catalog for Field 3, Pointing 2. 

Objects with a high probability of being a star (i.e., a low 
probability of being a galaxy, p^„i, as described in Coil et 
al. 2004b) were excluded when defining the pool of objects 
from which DEEP2 targets are selected. Figure 8 plots ob- 
ject radius rg versus apparent /?ab magnitude; where rg is the 
1-CT radius in pixels of a circular Gaussian fitted to the CFHT 
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Figure 8. Size-magnitude diagram resulting from peat photometry in CFHT 

pointing 32, wliich lias median seeing for the CFHT data set (0 . 84). Here 
Tg is the RMS radius in pixels of a circular Gaussian fit to an object's image 

(1 CFHT pixel = .207). The stellar locus is visible as the narrow horizon- 
tal band below and to the left of galaxies. The left vertical red line is where 
stars begin to saturate, and the right vertical line is the magnitude limit of 
the DEEP2 survey, Rab = 24.1. The blue horizontal lines indicate the 95% 
size range for for bright stellar sources in this pointing; it is determined in- 
dividually for each CFHT pointing. The orange slanting line marks a line of 
constant surface brightness, 26.5 mag/D", above which objects are excluded 
from DEEP2 targeting. Most exluded objects prove to be multiple objects 
incorrectly identified as single in the peat photometry (see § 12). Stars begin 
to blend with galaxies near R ~ 22.5. A Bayesian probability for star-galaxy 
separation, < pgai < 1, is calculated for all objects with size below the 
maximum of the stellar band by combining information on the magnitude 
and color distributions of stars and galaxies in DEEP2 (cf. Figure 9); objects 
with Pgai < 0.2 are colored purple in this plot. §6.1 provides more details on 
the algorithms used. 



photometry. The plot shown is for DEEP2 Field 3, pointing 
2 (also referred to as pointing 32), which has median seeing 
for the CFHT/DEEP2 data (0'.'84 FWHM); however, appai-ent 
object sizes are analyzed separately for each CFHT pointing. 
Stars are identifiable as the relatively tight, horizontal locus 
of points below and to the left of galaxies. Stellar radii in- 
crease at bright magnitudes due to saturation; the vertical line 
marks the point at which saturation becomes detectable in this 
particular pointing. Saturation sets in at Rab ^ 16.7- 18.0 
depending upon pointing, establishing the bright magnitude 
Umit for DEEP2 targets, Rab > 18.5. 

Information from the size-magnitude diagram in Figure 8 
is used in conjunction with the B-R vs. R-I color-color di- 
agram to compute pgai- First, objects in any part of Figure 8 
that have colors unlike any other stars or bright galaxies are 
given Pgai = 2. Second, the upper horizontal dashed line de- 
notes the 95% upper radius range for stars (determined sep- 
arately for each CFHT pointing), and all sources above this 
line are considered extended and assigned pg„i = 3. However, 
stars have sizes indistinguishable from galaxies at dim magni- 
tudes, below Rab ^23. In this regime, we can use the fact that 
stars tend to be brighter than and occupy a different locus in 
color-color space from compact galaxies to differentiate the 
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Figure 9. BRI Color-color diagram for stars and galaxies. Grey points are 
candidate target galaxies in the Extended Groth Strip; i.e., they are objects in 
the pcaf with 18.5 <i?AB <24. 1 andpga/ > 0.2. Stars, defined here as objects 
with p^cil < 0-2, are indicated by the red points. Straight line segments show 
the boundary used for color pre-selection to screen out low-redshift galaxies 
in Fields 2, 3. and 4; candidate objects at high redshift lie below and to the 
right of these lines. Specifically, the lines show the locus in color-color space 
where objects are given 10% weight for selection in Fields 2-4, the minimum 
possible for inclusion in the DEEP2 sample. The weight falls off from a 
maximum of 1 (a little below and to the right of this locus) as a 2-d Gaussian 
having a = 0.05 mag in both coordinates. This "pre-whitening" of the bound- 
ary is done to reduce the impact of any systematic color errors within a CFHT 
12K pointing, which are of order 0.02-0.04 mag (Coil et al. 2004b). The divi- 
sion between objects treated as low- and high-redshift in EGS is similar to the 
lines shown (cf. §7.2). Stars are distinguished from galaxies by computing 
a Bayesian probability, p^„i, based on the distiibutions of stars and compact 
galaxies in both magnitude and color-color space. §6.1 gives more details. 

two classes. The differences in colors are illustrated in Fig- 
ure 9, which shows the BRI color-color diagrams for stars and 
extended galaxies in a DEEP2 peat. For all objects having a 
size consistent with the stars in a given pointing, the value of 
Pgai indicates the net Bayesian probability that an object is a 
galaxy (as opposed to a star) based on both color information 
and brightness; it is therefore a value between and 1. The 
pink points in Figure 8 are the objects identified as likely stars, 
defined as those having pgai < 0.2. 

To summarize: pgai = 3 indicates an extended source, all 
of which are included for target selection. Objects with pe- 
culiar colors have pgai = 2; these are also all taken as candi- 
dates (they vast majority have turned out to be objects with 
large photometric errors). Finally, values of pg^i from to 
1 are used for compact objects with a size consistent with 
stars. More details are provided in Coil et al. (2004b). All 
objects with pg^i > 0.2 are included as candidate targets for 
DEEP2, with selection weight proportional to pgai (to a max- 
imum weight of 1). Further details on sample selection are 
given in Section 6.2, and possible biases are described in Sec- 
tion 8. 

Note that the trajectory of stellar colors in Figure 9 is ap- 
proximately described by two straight lines with a "knee" sep- 
arating the two. Demanding that this knee be identical in all 
fields allows us to place all CFHT 12k pointings on a com- 
mon color system, with an overall zero point determined via 
comparison to stars in the Sloan Digital Sky Survey. From a 



variety of tests, we find that pointing-to-pointing variations in 
the average photometric zero points are at the < 0.01 magni- 
tude level, while variations within pointings and absolute zero 
point uncertainties are ^ 0.02-0.04 mag (Coil et al. 2004b). 

6.2. Target Pool Selection Procedure 

With Pgai determined, a pool of potentially acceptable 
DEEP2 target galaxies can be defined. The final slitmasks are 
designed by selecting objects from this pool and placing them 
on slitlets. All galaxies selected as targets must fulfill three 
separate selection criteria that apply in all 4 DEEP2 fields. 

First, the object must be classified as a potential galaxy 
by our probabilistic star-galaxy separation procedure. All 
extended objects and objects with peculiar colors are auto- 
matically included (pg^i = 2 or 3). Objects whose sizes are 
consistent with stars are taken as eligible targets if they have 
Pgai > 0.2, i.e., the probability of being a galaxy is greater 
than 20%. The distribution of Pgai is extremely bimodal, with 
values piling up around or 1 (Coil et al. 2004b); as a result, 
varying the cut level between 0.2 and 0.8 makes very little 
difference to the sample. 

The next cut is on apparent /?-band magnitude: to be con- 
sidered for DEEP2 targeting, an object must lie in the range 
18.5 < /?AB < lA.l?^ For science purposes, it would have 
been more ideal to use a redder band like / to set the survey 
magnitude limit, but the CFHT / photometry obtained was 
less uniform in depth, necessitating the adoption of R?^ 

Finally, the object must have surface brightness above a 
standard limit (indicated by the slanting line in Figure 8). This 
surface brightness (SB) requirement is defined by the equa- 
tion: 



5B = /?AB + 2.51ogio[7r(3r,)2] <26.5, 



(1) 



where rg here is the Gaussian profile radius for that object, 

given in arc seconds. The minimum rg is set to .33, so that 
objects with 3rg < 1" are presumed to have an effective 3(t 
radius of 1" minimizing the effect of noise in measuring the 
sizes of compact objects on this cut. Visual inspection of the 
low-surface-brightness objects which are excluded by this cut 
indicates that nearly all are double or multiple sources in the 
BRI images that were incorrectly identified as single objects 
in the peat due to blending. Their surface brightnesses are 
low because their radii are falsely inflated. Since the data for 
multiple objects are ambiguous, they are not good targets for 
DEEP2, and this cut properly excludes many of them. Further 
discussion of multiple interlopers is given in §12. 

6.3. Color Pre-selection in Fields 2, 3, and 4 

In Fields 2, 3, and 4, the target pool is further refined to 
remove galaxies with redshifts likely to be below z ^ 0.75. 
All but the bluest galaxies exhibit a significant break in their 
spectra at A ~' 4000 A. This break causes a strong separation 

'* The DEEP2 BRI photometric system is an approximately AB magnitude 
system in which the observed magnitude in each filter has been zeropointed 
by requiring that it matches SDSS photometry in the closest available band 
for stars with zero AB color (Coil et al. 2004b): e.g. we choose a zero point 
such that DEEP2 R equals SDSS r plus a color term proportional to SDSS 
(r-i). To the degree to which the SDSS system is on AB, DEEP2 magnitudes 
will be AB magnitudes within the native filter-l-telescope system; however, 
this assumption is imperfect at the ~ 0.02 mag level (D. Eisenstein, private 
communication). System response curves including the filters and detector 
may be found at http://deep.berkeley.edu/DRl/photo.primer.html. 

'' Ironically, the originally-planned fields with worst / ended up being 
discarded due to the descope from the originally-planned coverage. 
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Figure 10. Assessment of the accuracy of the DEEP2 BRI color-color pre- 
selection for rejecting low-redshift (z < 0.7) galaxies. Grey points are the 
same set of DEEP2 candidate galaxies shown in Figure 9; red points are those 
EGS galaxies whose spectroscopic redshifts are below z = 0.7 (only secure, 
Q = 3 and Q = 4 redshifts are used). We here take advantage of the fact that 
low-redshift galaxies are not excluded from DEEP2 targeting in the Extended 
Groth Strip. The boundary is the 10%-weight locus repeated from Figure 9. 
Nearly all red points lie above the boundary, indicating that they would be 
successfully excluded from DEEP2 by the color-color pre-selection cut. 



in color-color space between objects where the break falls in 
the B band or on the blue side of R, versus those where it oc- 
curs at redder wavelengths. We therefore use galaxies' B-R 
and R-I colors to identify the low-z interlopers. Since the 
bluest, flattest-spectrum objects exhibit minimal break (hav- 
ing colors near at all redshifts), we always include those 
objects within our selection region, in order to ensure that our 
sample is complete for galaxies at z > 0.75. 

Figure 10, Figure 11, and Figure 12 show color-color plots 
of galaxies in different redshift regimes in the Extended Groth 
Strip, where we obtained spectroscopy of galaxies at all red- 
shifts (i.e., galaxies are not excluded based on this color cut). 
The lines show the nominal color boundary used to screen out 
low-redshift galaxies - essentially, only objects to the right 
and below this line are chosen in Fields 2, 3, and 4, while 
objects of all colors are observed in EGS. 

The boundary used is not perfectly sharp but rather the 
weight given an object for selection falls off (from 1) as a 
two-dimensional Gaussian with a = 0.05 mag in both coordi- 
nates; we call the resulting weight the "color-weighting fac- 
tor" for Fields 2, 3, and 4 in the discussion below. This "pre- 
whitening" is done to reduce the impact of possible systematic 
color errors within a single CFHT 12K pointing, which are 
of order 0.02-0.04 mag (Coil et al. 2004b). The Unes shown 
are not the nominal (i.e., sharp) boundary but rather show (for 
R-I > 0.2) or approximate (for R-I < 0.2) the locus of points 
where the color- weighting factor equals 0.1 = 10%, the mini- 
mum for objects selected for DEEP2 spectroscopy. Although 
this weight is used in target selection, selection probability is 
a slow function of the weight^'^. The boundary in these figures 
corresponds to selection criteria for DEEP2 targets of: 

^^ The exact position of these lines changed slightly after 2002, taking 
account of lessons learned from the first semester's data. 



Figure 11. A test of the accuracy of the DEEP2 BRI color-color pre-selection 
for galaxies near the transition redshift. This figure is similar to Figure 10, 
but colored points indicate EGS galaxies with spectroscopic redshifts in the 
range 0.70-0.75. Most of these lie below the boundary, indicating that they 
would be accepted by the color-color pre-selection cut despite their redshifts 
below 0.75. This is conservative, ensuring that objects with redshifts above 
z = 0.75 are not lost. 



Pi 
I 
m 2 





.'-'■' 1 '' - 




z>0.75; .. 


r 


Low Z i . . : 


L 










■■ . ■-■•■-j]gS^*''S-5'\i,i- / ■'..,'■.■ .f^ ■''■•' • 










— 


■" "''-^iSS&^Mvij^f'" '' > ■ - r""" ■" ' ~ 




\i^^^^M-' \ 




'^" ^^"^^ffl^P- "''''^'''•'^^^^W^^'^'''''- 




-, ' ■•:'.• • ijUh-t^Hlfl^k ■ y , i^^aSJ^Kg^pffiKTif ■ r-' : 








^^j^fei^^^^^'^-'-- Highz ; 




. ■■ '".(t^'X.-^i^^^W'ff;:' '■■' 




1 1 1 1 1 1 1 '1 1 1 1 1 1 1 1 1 1 1 1 1 ' 



0.0 



0.5 1.0 

R-I 



1.5 



Figure 12. The accuracy of color-color pre-selection for known high-redshift 
galaxies. This figure is similar to Figure 10, but colored points indicate EGS 
galaxies with spectroscopic redshifts above 0.75. Nearly all of these lie be- 
low the boundary, indicating that nearly all would be accepted by the color- 
color cut. The DEEP2 color selection produces a sample which is highly 
(approaching 100%) complete for high-redshift objects. 

1. (B-/?)< 0.5; and/or: 



2. (R-I)> 1.1; and/or: 

3. (B-R) < 2.45 X {R-I)-0.2916 . 

The selection weight is 100% for objects with (B-R) < 0.389; 
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Figure 13. Model color-color distributions for redshifted spectral energy dis- 
tributions (SEDs) superimposed on DEEP2 targets. The SEDs used are based 
on spectroscopy of normal nearby galaxies from Coleman. Wu, and Weed- 
man (1980) and selected starburst galaxies from Kinney et al. (1996). Each 
SED executes a loop as redshift increases as the Balmer/4000-A break moves 
first through B-R and then through R-I at higher redshifts. Each model 
loop runs from z = to z = 2.0 with diamonds shown at increments of 0.2; 
the line type switches from dashed to solid at z = 0.7. The lines are colored 
to improve visibility, with eai'ly-type spectra in red, later-types in puiple, and 
starbursts in blue. The dashed lines are the boundaries used to pre-select tar- 
get galaxies with z. > 0.75 in Fields 2. 3, and 4. Arrows indicate the impact of 
/i\/ = 0.2 of extinction at a variety of redshifts. Dust reddening tends to move 
galaxies parallel to, rather than across, the DEEP2 color cut. These models 
are illustrative only - the final positions of the boundaries were tuned using 
EGS data, which contain galaxies at all redshifts. 

(R-I) > 1.211; and/oi (B-R) < 2.45 x (R-I)-0.311. 

Figure 13 demonstrates why we would expect such a color 
pre-selection to be effective. This figure shows the tracks 
of redshifted galaxy SEDs from Coleman, Wu & Weed- 
man (1980) (augmented by starburst galaxies from Kinney 
et al. 1996) in color-color space. Galaxies move across the 
color-color boundary rapidly as their Balmer and 4000-A 
breaks move first through B-R and then through R-I near 
z = 0.7. The resulting valley in the colors can be used to sort 
galaxies into two groups below and above z ~ 0.75. 

The CWW models are illustrative only - the actual bound- 
ary was tuned using DEEPl data and early redshifts in Field 
1 (EGS), where galaxies are targeted regardless of their color. 
Figure 14 shows the final redshift histograms of objects pass- 
ing and failing the color cut in EGS. The vertical dashed line is 
the desired redshift cut at z = 0.75, and the 50-50% crossover 
point is at z = 0.714. The crossover has been placed signifi- 
cantly below the target cut in order to enhance sample com- 
pleteness at z > 0.75 (at the expense of somewhat decreased 
efficiency at targeted only high-z objects). Figure 15 shows 
the targeted fraction at all redshifts, which indicates the com- 
pleteness of our color cut at a particular z. For the bin cen- 
tered at z = 0.75, the pre-selection color cut correctly captures 
87% of all objects and incorrectly rejects 13%. For bins above 
z = 0.8, these fractions are constant at 97.5% and 2.5% respec- 
tively. Figures 10 and 12 show that the loss of high-z galaxies 
and the creeping in of low-z galaxies are primarily due to blue 
objects in the lower-left comer of the color-color diagrams, 
whose colors tend to blur together at low and high redshift, as 



Figure 14. Redshift histograms for EGS galaxies with secure (2 = 3 or 4) 
redshifts that either pass or fail the pre-selection color cut. We correct here 
for the selection probability of each object, as otherwise the S-magnitude 
weighting scheme used in EGS would influence the redshift distributions (see 
§7). The 50-50% crossover point has been deliberately placed at z ~ 0.71 
to insure high sample completeness above the nominal target redshift cut at 
z = 0.75 (vertical dashed line). 
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Figure 15. The fraction of EGS galaxies passing the pre-selection color cut 
relative to the total number of galaxies with measured redshifts, as a function 
of z. Numbers are taken from Figure 14. This figure illustrates the high 
completeness of the color cut, which is 91% at z = 0.75 and averages 97.5% 
for0.8<z< 1.4. 

well as to objects with catastrophic photometric errors. 

The net effect of the color pre-selection in Fields 2, 3, and 4 
is to remove ^ 55% of R < 24.1 objects from the target list - 
almost all unwanted foreground galaxies - with very minimal 
loss of the desired distant galaxies, and thus to increase the 
survey's efficiency for studying objects at the target redshift 
of z ~ 1 by a factor of ^ 2.2 at the cost of a ^ 3% loss of 
desirable targets. 
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Figure 16. Schematic of tiie DEIMOS focal plane, showing the geometry of 
the 2x4 CCD mosaic. The dispersion directions of the individual spectra run 
vertically. Superposed is a slitmask image as projected onto the detector in 
direct imaging mode. The TV guider camera stares directly at a fixed of the 
focal plane. Light is reflected into the camera from the shiny slitmask (upper, 
light grey region) and more efficiently from a reflective mirror (lower, dark 
grey region). The geometry is precisely known, so that placing a given star at 
a predicted TV guider pixel at the proper spectrograph angle should succeed 
in placing the alignnment stars for that mask within their respective 4"x4" 
alignment boxes. Alignment proceeds by taking successive direct images 
through the mask followed by small corrections to the telescope RA/Dec and 
spectrograph position angle. Flanking the main detector array are the flexure 
compensation CCDs. These are fed by ThAi' emission spectra, which are 
piped into the slitmask focal plane by optical fibers (there are actually two 
fiber spectra at each end, but only one is shown). The ThAr spectra are so 
rich that at least one emission feature falls on the FCS CCDs at all grating 
tilts. 



7. SLITMASK DESIGN 

7.1. DEIMOS Detector and Slitmask Geometry 

Having described the factors which determine whether an 
object is targeted by DEEP2, as well as the factors which con- 
tribute to their selection weight, we now turn to the methods 
used to design individual slitmasks from this target pool. 

The design and selection of objects is intimately inter- 
twined with the geometry of the DEIMOS detector array and 
slitmasks, which are illustrated in Figure 16. Related in- 
strument parameters are given in Table 2. The telescope fo- 
cal plane and slitmask geometry are shown as they are im- 
aged onto the CCD detector array. For simplicity, we adopt 
the same coordinate system convention used for DEIMOS 
pipeline software outputs: "columns" run horizontally in the 
figure, i.e. parallel to the slitlet direction, while "rows" run 
vertically in the figure, parallel to the spectral direction. A 
typical extracted spectrum will therefore have 8192 columns, 
but considerably fewer rows.^^ 

The CCD detector is composed of eight 2K x 4K MIT Lin- 
coln CCDs with 15-/X pixels, for a total size of 8192 x 8192 
pixels and total physical dimensions 126 mm x 126 mm. The 
CCD thickness is greater than normal (40 /i rather than 15 /i), 
which renders fringing negligible at wavelengths below 8000 
A and still yields low fringe amplitudes at longer wavelength 
(e.g. only ±2% at 9000 A). The thickness of these CCDs also 
boosts total system quantum efficiency at far red wavelengths, 

^' Note that the reduced data keyword hcol as explained in § 1 3 refers to 
readout columns on the individual DEIMOS CCDs. which run parallel to the 
dispersion direction, opposite from our standard convention. 



e.g. to 16% at 9000 A.'*" The between CCDs in the spectral di- 
rection is only seven pixels wide, while the gap between the 
long sides of the chips is roughly 1 mm, or 8". A single spec- 
trum in general falls on only two CCDs. 

The slitmasks are flexible flat aluminum sheets that are 23.6 
cm wide, 75.7 cm long, and 0.50 mm thick. The total length 
of the mask on the sky is 16 .7, while the maximum usable 
slit length (subtracting the CCD gaps) is 16.3. The masks 
are milled flat on a computer-controlled milling machine and 
spring-loaded onto a cylindrically curved mandrel in the spec- 
trograph that approximates the spherically curved telescope 
focal surface. Eleven slitmasks can be held in a jukebox-like 
cassette at one time. 

7.2. Mask Design Algorithm 

The widest part of a DEIMOS slitmask projects to 5 .2 
across on the sky, but we place slitlets only within the in- 
ner 4' because the outline on one side is irregular. The wave- 
length range of the portion of the spectrum from a given slitlet 
which actually falls on the DEIMOS CCDs will vary accord- 
ing to slitlet placement. To maximize the wavelength range 
that is common to all objects (given anamporphic shifts in 
central wavelength with position on the sky), we select targets 
for spectroscopy only within a strip which is curved slightly 
"backwards" on the sky, but whose edges correspond to lines 
of constant wavelength on the detector. To ensure good sky 
subtraction, we require that slitlets generally be at least ^ 3" 
long, implemented by requiring that successive targets be at 
least 3" apart in the long direction of a DEIMOS slitmasks 
(i.e. along a 'row'). Objects may not be placed within 0. .3 of 
the end of a slitlet, in order to ensure uniform slit width along 
each object for sky subtraction. We also require .6 of dead 
space in the slitmask's long direction between slitlet ends to 
limit cross-contamination between light from adjacent spectra 
(this is part of, rather than adding to, the minimum spacing be- 
tween objects). The average sHtlet length resulting from this 
procedure is 7". 

In Fields 2, 3, and 4, the pattern of slitmasks on the sky 
is chevron-shaped, as shown in Figure 1. The two arms of 
the chevron are observed either east or west of the meridian 
to better align the slitlets with the atmospheric dispersion di- 
rection. Masks are designed separately for each of these two 
sets of masks; where they overlap at the point of the "V", 
some objects will be multiply-observed, allowing tests of the 
repeatability of DEEP2 measurements. The slanting pattern 
also keeps the masks within the boundaries of the 28'-wide 
CFHT 12k pointings, despite the masks being > 14' long. 
Mask centers are separated along by 2' on average, which 
gives each object roughly two chances to be included on a 
slitmask. The actual step size is modulated by an adaptive 
tiling scheme (see below). 

The algorithm to place objects on masks in Fields 2, 3, and 
4 assigns each object a weight W between and 1 that is the 
product of four factors. The first factor, Wjc, comes from the 
star-galaxy separation probability pgai; it is unity if Pgai = 2 or 
3 but equals pgai if Pgai is in the range to 1 (i.e., for small 
galaxies or stars); only objects with pgai > 0.2 are considered 
for targeting. The second factor. We, which we will refer to as 
the "color-weight" is a measure of each object's consistency 



■*" See throughput plots at 
http://www.ucolick.org/~ripisc/Gol200/through_gol200_tilts.gif 
details. 



for 
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with the DEEP2 color-color cut. As explained in §6.3, this 
factor runs between and 1 and is the result of smoothing a 
sharp-edged step function in the color-color diagram by a 2- 
d Gaussian of width a = 0.05 mag in each coordinate. Only 
objects with color-weight > 0.1 are considered for targeting. 

The third factor, Wr, is a function only of R magnitude and 
is designed to flatten the distribution of objects in redshift 
and magnitude. In Fields 2, 3, and 4, this weight is given by 
the smaller of the two numbers (0.75 x io-0-4*(24i-fi) and 1, 
which yields a weight that rolls off smoothly from 1.0 to 0.75 
in the last 0.31 magnitude, de-weighting the faintest galaxies 
moderately. The final factor, applied only for DEEP2 objects 
observed in 2003 or later, is designed to reduce the number 
of high-redshift galaxies beyond z = 1 .4 targeted, as they of- 
ten fail to yield redshifts in DEEP2 spectroscopy. The tar- 
get selection weights for these "blue corner" objects (those 
with B-R<Q.5 and R-I < 0.45) are multipHed by a factor 
Whc that rises linearly from 0.1 for objects with R-I < 0.045 
to 1 for objects with R-I > 0.45. Hence, in fields 2-4, 
W = Wscy<WcXWRxWbc. 

In the Extended Groth Strip, pgai is still used for weight- 
ing, but no color pre-selection or blue-corner deweighting are 
applied, eliminating the second and fourth weighting factor. 
However, ignoring color information entirely yields a sam- 
ple which is dominated by faint foreground galaxies. There- 
fore, the /? -magnitude-based weight, Wr is modified in EGS 
to counteract this. 

The magnitude-based weight factor is calculated differently 
accordingly to which of two possible different scenarios ap- 
plies: 

• If /?< 21.5 or the galaxy passes the color cut to be at high 
redshift,"*' Wr takes the same value as it would in Fields 2, 3, 
and 4, fafling from 1 at /? < 23.8 to 0.75 at /? = 24.1. 

• If R > 21.5 and the galaxy fails the color cut (i.e., it is both 
faint and at low redshift), the weight is the lesser of the num- 
bers 0.1 X io-0 4(24.i-R) and 1. This yields a weight that rofls 
off smoothly from 1 .0 at 7? = 2 1 .6 to 0. 1 at /? = 24. 1 , decreas- 
ing the sampling rate for the faintest low-z galaxies. 

Thus, distant EGS galaxies are weighted very like those in 
Fields 2, 3, and 4 (except for small differences in the color 
boundary), and rare, intrinsically bright low-redshift galaxies 
are favored in target selection; so in EGS W = WSG x Wr, but 
unlike in Fields 2-4 Wr depends indirectly on BRI color. We 
summarize the effects of the two weighting schemes on the 
resulting sample at the end of this section after describing the 
actual mask design process."*^ 

In Fields 2, 3, and 4, objects are allocated among the slit- 
masks in a two-stage process"'-'. The first pass places slitlets 
only in a central strip that is 2' wide on average (though this 

*" This cut uses the nominal color-color boundary rather than a smoothed 
version; so objects with (B-R) < 0.5, (R-I) > 1.1, and/or (B-R) < 2.45 X 
(R — /) — 0.5 are treated as high-redshift in EGS, while those failing all these 
criteria are considered to be at low z- 

"'^ The weighting scheme and target selection in EGS is complicated some- 
what due to the fact the photometry in pointing CFHT 12K pointing 14 
was comparatively poor and less well-calibrated than in the remainder of the 
DEEP2 spectroscopic fields, which leads to less well-characterized B-R and 
i?-/ colors there. In this pointing, Pg„i is still retained as a factor, but the R- 
magnitude weighting is the same as in Fields 2, 3, and 4 regardless of object 
color, and there is no surface-brightness cut or blue-corner de-weighting. Im- 
proved photometry for this field will be separately released in an upcoming 
paper, Matthews et al. (2012). 

^^ This process begins after the selection of four to six bright stars to be 
used for slitmask alignment. Each star is then used as the center of a 4" X 
4" square alignment box. To the degree possible, these ai'e placed evenly on 
both ends of a mask for maximum leverage. 
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Figure 17. Blowup of mask designs in CFHT pointing 32 showing a region 
where many masks overlap, near the "V" of the chevron pattern. Dotted lines 
show the prime inner ~2' strip within each mask (actual widths vary due to 
the adaptive tiling algorithm). Target objects are indicated by black ellipses 
with size 3rg and axis ratio (1 -e^). Slitlets on a given mask have the same 
color. Alignment star's are indicated by the red asterisks. 



width can vary due to adaptive tiling, q.v. below). This cen- 
tral, curved strip is favored in order to maximize the wave- 
length region common to all spectra (with our chosen grating, 
the wavelength shift across 2' is ^200 A). 

Cosmic variance in the number of targets within a 16' x 2' 
region is high, so that some masks will have fewer slitlets than 
others if masks are evenly spaced. To minimize this variation, 
we iteratively adjust the positions (and correspondingly the 
widths of the central regions) of masks in Fields 2, 3, and 4 
such that the number of targets selected in the first pass is held 
roughly constant from mask to mask (see Figure 17, which 
reflects these variations). This "adaptive tiling" strategy helps 
to ensure that a uniform fraction of target galaxies is sampled 
over the survey region. 

At the beginning of the first pass, the list of targets is tem- 
porarily trimmed by generating a random number between 
and 1 and retaining only objects for which this random num- 
ber is less than their weight. Thus, low-W objects (e.g., likely 
stars and galaxies just outside the nominal color cut) are only 
rarely placed on slits in this first pass. The trimmed list is then 
searched for cases where we can place two objects (separated 
by less than 3" in the long direction) on a single slitlet that 
has a PA relative to the long axis of the mask of less than 30°. 
Slitlets are allocated first to all such cases, boosting our ability 
to study close pairs of galaxies. 

A random priority uniformly distributed between and 1, 
Pi, is then generated for each of the remaining objects in the 
trimmed list. Slit allocation then begins at one end of the mask 
(along its long axis) and proceeds to the opposite limit. Each 
object that can be observed without precluding the observa- 
tion of any other is selected for observation, while in cases of 
conflicts the object with greatest Pi is taken. A new object 
will be allocated a slitlet only if it is more than a minimum 
distance along the long axis of the mask away from any ob- 
ject already targeted. This minimum buffer distance is the 
greater of 3"or (0 .3 -i- twice the radius of the new object pro- 
jected along the slit, r^proj), where the latter is computed from 
the new object's Gaussian radius r^, axis ratio, and PA. New 
objects that are large must therefore be farther from existing 
objects and will be allocated longer slitlets as a result. 

This procedure (initial selection according to W , and nearly 
random selection amongst the surviving targets) yields a sam- 
ple where the probability of selection is almost perfectly pro- 
portional to an object's weight, making the best possible use 



18 



of the central portion of each mask. In Fields 2, 3, and 4, 
77-80% of slitlets are allocated in the first pass. 

After the first pass through all masks is completed, the re- 
maining galaxies that overlap a given mask (over the full 16' 
X 4' area) are given new weight values W2 according to their 
initial selection weights, with j^5^ - This new weight ac- 
counts for the fact that the probaliility a given object is se- 
lected in pass 1 is very close to 0.55W, so the probability it is 
not selected is 1 -0.55W. A new random number uniformly 
distributed between and 1, P2, is then generated for each 
object; all potential second-pass targets are then ranked in de- 
scending order of ^, providing a priority to each one. If all 
objects in pass 2 were subject to conflicts (i.e., have to com- 
pete with other objects for slitlets), this priority would ensure 
that the final probability an object is selected after passes 1 
and 2 is proportional to its initial weight, W. 

Objects are then allocated slitlets in descending order of 
their priority, subject to the requirement that so long as they 
do not conflict with any slitlet that has already been added to 
the mask (e.g. as an alignment star, in Pass 1, or as an object 
with higher priority in Pass 2). In the end, any target that can 
be placed on the mask without causing a conflict is included, 
so in the second pass, low-weight objects will frequently be 
selected if they are cost-free. This causes the final selection 
probability of an object not to be simply proportional to W; 
in fields 2-4, the probability is instead closely approximated 
hy p,,uc,ion = 0.27976 + 0A47nW -0.09137W^. Hence, an 
object with weight of 0.1 will be selected for observation 32% 
of the time, while one with W = 1 will be selected 64% of the 
time. Since most eligible targets have a weight near 1, the 
overall selection rate is close to the latter value. The second 
pass proceeds from west to east, assigning targets to all the 
masks in a row (i.e., all the masks with similar PA) in order 

The choice of targets in the Extended Groth Strip differs in 
two respects from Fields 2, 3, and 4. First, as noted above, 
there is no pre-selection based on color, and the /^-magnitude 
weighting factor used is different in order to suppress the 
number of faint foreground galaxies. The tile pattern in EGS 
is also different owing to the narrower width of the region, 16' 
rather than 30'. Masks are divided into 8 groups with 15 masks 
per group, each group approximately square in size (see Fig- 
ure 2). In each square, eight masks have their long dimensions 
perpendicular to the long axis of the strip, and seven masks 
have their long axes parallel to the strip; masks parallel to the 
strip (within a given pointing) are given larger mask numbers. 

In order to maximize the number of galaxies having mea- 
sured rotation curves, we first divide all extended, elliptical 
objects (those with eccentricity e > 0.081 and size 3rg > 1 
arcsec) into two sets, one consisting of objects with their ma- 
jor axis within 45° of the long axis of the strip and the other 
including the remainder. The two-pass maskmaking algo- 
rithm described previously is then run on the set of objects 
in the second category, combined with with those which are 
too small or round to have a well-defined position angle, to 
assign objects to masks with their long direction perpendic- 
ular to the strip. Then, any galaxies not assigned to a mask, 
as well as the extended objects with long axis along the strip, 
are used to design the masks with long dimension parallel to 
the strip, again using our standard slit assignment algorithm. 
The crossed orientations again permit observations both east 
and west of the meridian, as well as enabling rotation curve 
measurements. However, in EGS the step size between masks 
is only 1 ' owing to the greater density of targets on the sky. 



and here it is not varied adaptively. 

In EGS, 74% of objects on slitmasks are selected in the first 
pass, which again only places targets in the central 2' width 
of a mask. The second pass is similar to the second pass in 
Fields 2, 3, and 4; again galaxies are selected over the full 
4' active width. By limiting the first pass selection to a field 
of only 2', the spectral coverage of the dispersed light for the 
galaxies is made more uniform. In the EGS, the probability 
an object is selected for observation depends on the weight, 
W (which here is the product of Wsc and the magnitude- 
dependent weight Wr) as prelection = 0.33398 + 0.42687W; so 
an object with W = 0.1 is selected 38% of the time, while 76% 
of objects with W = 1 are targeted. 

Once all objects are chosen, the final step is to set the 
lengths of the slitlets. Since no account was taken of the sizes 
of previously selected objects when new ones were added, it 
may happen that two projected buffer distances, each of size 

2rg proj+0 .3, overlap for two neighboring objects. In this case 
the available space is evenly allocated between the objects 
regardless of which one is bigger. Much more commonly, 
when the separation is larger than the two buffer distances, 
the available slit lengths are allocated in proportion to their 

values of rg^proj- In both cases, half the dead space of .6 
between slitlets is subtracted from the extent of each slitlet in 
the direction of overlap. 

Once slitlets have been placed on targets, the design of all 
DEEP2 masks are adjusted to include ^8 sky-only slitlets 
placed in open spaces on the masks. Their positions along 
the long axis of the mask are chosen to fall in between the 
objects with largest separations, while their positions along 
the short axis are chosen by choosing empty regions from 
the CFHT 12K imaging. The sky slitlets are used to perform 
fallback, "non-local" sky subtraction useful in the presence 
of certain rare instrumental anomalies and for the occasional 
very short slitlet. These sky slitlets are only .7 wide; to pro- 
duce a nonlocal sky model, the spectrum from a sky slitlet is 
convolved with a kernel that varies with wavelength, fit us- 
ing an individual slitlet's estimated sky spectrum. In practice, 
this method proved inferior to local sky subtraction using the 
object's own sky spectrum even for short slitlets, and it has 
not generally been used extensively. However, non-local-sky- 
subtracted versions of both 1-d and 2-d spectra are included 
in pipeline-produced data files for completeness. 

As a final step, in order to preserve rotation curve infor- 
mation, the long axis of each slitlet is oriented as closely as 
possible along the major axis of its galaxy (as measured from 
the CFHT photometry) provided the following conditions are 
met: 1) the axis is within ±60° of the long axis of the mask; 
2) the galaxy has rg > 1"; 3) the slitlet has not already been 
tilted to observe another object; and 4) the galaxy is mea- 
surably ellipsoidal with CFHT axis ratio less than 0.85. As 
slitlet tilts, 9, can be up to 30° away from the long axis of the 
mask but no more, the mismatch between slitlet and galaxy 
PA can be up to 30° (note also that PAs measured from high- 
resolution ACS images can differ considerably from the PAs 
from ground-based images, as found by Wirth et al. 2004). 

In the central portions of a slitmask, slitlets are still tilted 
by ±5° compared to the mask direction for objects which do 
not fulfill these requirements. This ensures better wavelength 
sampling of the night sky spectrum, improving sky subtrac- 
tion (optical distortions provide that sampling for untilted slits 
at the ends of a mask, so those slitlets are designed with PA=0 
relative to the mask). The crosswise widths of tilted slitlets 
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are reduced by cos(9) in order to maintain a fixed projected 
slitwidth along the dispersion direction. 

Figure 17 shows a blowup of a region covered by several 
slitmasks in Field 3, indicating the central first-pass regions, 
eligible targets, selected targets, and slitlets with their final 
lengths and orientations. The overlap of slitmasks at the "V" 
of the chevron pattern provides a sample of multiply observed 
objects for statistical and reliability assessment purposes. 

In the end, 61% of eligible targets are assigned slitlets in 
Fields 2, 3, and 4. More precisely, this is the fraction selected 
out of targets that have 18.5 < R < 24.1, pg^i > 0.2, pass the 
surface-brightness cut, and have probability > 0.1 of meeting 
the color-color cut. Slitmasks in these fields contain ^110- 
140 target slitlets (mean 121) that are typically 7" long but 
span the range 3-15H-". 

In EGS, the average number of slitlets per mask is 148. In 
pointings 11-13 of EGS, /?-band de- weighting of low-redshift 
galaxies boosts the relative number of high-z galaxies, with 
the result that roughly 73% of targets out of a Field 2-4 equiv- 
alent sample (as defined above) are selected in EGS^'* In con- 
trast, due to the deweighting of faint objects that fail the color 
cut, only 53% of low-z galaxies are selected. Without this 
deweighting, only 55% of the galaxies in a Field 2-4 equiv- 
alent sample would be selected in EGS, making it sparser in 
number density than other fields. 

8. SELECTION EFFECTS IN THE FINAL TARGET SAMPLE 

The final DEEP2 target sample suffers from a variety of 
known selection effects. These effects should in general be 
taken into account and/or corrected for when utilizing the 
DEEP2 dataset. A variety of methods for doing so have been 
implemented in DEEP2 science papers; see §3 and Table 1 for 
examples. 

• Galaxy color bias due to the R magnitude limit Selection 
according to a fixed apparent /?-band magnitude limit causes 
the resulting sample to reach different depths (in terms of ab- 
solute magnitude) for galaxies of different intrinsic colors at 
the same redshift. At z ^ 0.5, the /?-band filter corresponds 
closely to restframe B, so objects are selected down to the 
same B-band absolute magnitude independent of color Below 
that redshift, /?-band selection favors redder galaxies, while 
above this, it favors bluer galaxies. The effect was first illus- 
trated in Willmer et al. (2006) and is discussed in more detail 
in §15. Though fundamental, this bias is well understood and 
can readily be applied to mock catalogs and other statistical 
predictions. The same type of bias will inevitably affect any 
magnitude-limited sample which spans a range of redshifts, 
due to the broad span of intrinsic galaxy colors. 

• Loss of bright star-like objects The use of Pgai discrim- 
inates against point sources (i.e., objects with apparent size 
consistent with stars; cf. Figure 8), and thus potentially dis- 
criminates against bright Type 1 AGNs, i.e., QSOs. AGNs 
will be lost only if they are both point-like and lie on or near 
the stellar locus in the color-color diagram. As a result. Type 1 
AGN will be included in the DEEP2 sample at some redshifts, 
but excluded at others. 

• Misclassification of faint stars as galaxies Although the 
differential number counts of objects classified as stars (pgai < 

^'^ It is thus possible to simulate a high-z sample from these three EGS 
pointings by drawing a subset of the high-z galaxies at random, compensat- 
ing for the dependence of weight on color in the other fields. The resulting 
sample will be essentially identical to that in Fields 2, 3, and 4 except for 
slightly different statistics in a small corner of color space near the color-cut 
boundary. 



0.2) in the DEEP2/?cflfs rise approximately linearly at bright 
magnitudes, they flatten out at Rab ~ 22.5; it is likely not co- 
incidental that this is where large numbers of galaxies begin 
to populate the stellar regime in Figure 8. It appears that the 
Pgai values overestimate the probability that faint objects are 
galaxies rather than stars; we therefore tend to include more 
stars than expected within the DEEP2 target list. This is con- 
servative from the standpoint of galaxy selection, as the con- 
taminant objects are trivial to weed out once spectra are ob- 
tained.'*^ However, as a result of this misclassification, a list 
of photometric stars selected purely via their pgai values will 
be incomplete at the faint end and would be of limited use for 
stellar studies. 

• Loss of objects due to missing B and I photometry The 
definition of pgai and the DEEP2 color pre-selection both de- 
pend on having measured B-R and R-I colors.'*^ Though the 
/?-band photometry is complete to the /?ab = 24. 1 magnitude 
limit in all fields (Coil et al. 2004b), very red or blue objects 
may not be detectable in either the B or I imaging. However, 
in the three EGS peat pointings with accurate B and / pho- 
tometry (11-13), which contain 59,594 objects brighter than 
/?AB = 24.1, only 153 galaxies are missing B. They tend to be 
faint and red and constitute only 4% of distant (red-sequence) 
galaxies with /?ab > 23 and R-I > 1.0. Their loss is not 
significant. Similarly, 71 galaxies with /?ab < 24.1 are not 
detected in /; they are predominantly very blue. 

• Loss of small, distant, faint red galaxies Discriminating 
faint stars from galaxies depends on the contrast between stars 
and galaxies in the color-color plot (Figure 9) (as well as on 
the fitting functions used for the magnitude distributions of 
each sample). The blue stellar locus is quite narrow in this fig- 
ure, so the loss of blue galaxies misclassified as stars is frac- 
tionally small. The situation changes where B-R > 1.5 and 
R-I > 0.6, where the surface density of stars becomes more 
comparable to that of galaxies. As shown by Figure 13, this 
region of the color-color plot is populated by intrinsically red 
galaxies beyond z ^ 0.7. Loss of such galaxies was estimated 
by Willmer et al. (2006) using the higher-resolution HST ACS 
images in EGS and partial redshift data. We redo that calcu- 
lation using final EGS redshifts and a larger database. 

Figure 18 plots the color-color diagram for EGS but re- 
stricts to objects in the ACS mosaic fainter than R = 23.0, 
where misclassification is most likely (cf. Figure 8). DEEP2 
stellar candidates (pgai < 0.2) that are confirmed as stars in 
EGS ACS images are plotted in blue; those shown instead to 
be galaxies are red. The latter are missed by DEEP2 and are 
thus a concern. Galaxies are lost in a narrow band in B-R, 
which corresponds to certain combinations of intrinsic red- 
sequence color and redshift (cf. Figure 13). 

The central rectangle in Figure 18 delineates the region in 
which galaxies are lost. Setting true stars aside (the blue 
points) leaves a total of 99 galaxies in the rectangle, of which 
27 were misclassified as stars by pgai (red points). How- 
ever, only about half of all R > 23.1 red-sequence galaxies 
at z = 0.75-1.1 lie in the rectangle - the others lie in the rect- 
angles above or below - and this fraction is found to be vir- 
tually constant at all redshifts above z = 0.75. Since galax- 
ies outside the central rectangle are not lost, the net result is 
that ^ 0.27 X 0.50 ^0.13=13% of red-sequence galaxies in 

'^^ Overall, 1.8% of the objects targeted by DEEP2 turned out to be stars; 
69% of those stars have R > 22.5. 

^* Note: colors are used to select objects in Fields 2, 3, and 4 and (less 
directly) in Pointings 11, 12, and 13, but not in Field 14; see §6.3. 
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Figure 18. Color-color diagram for faint (R > 23) objects in the Extended 
Groth Strip with HST/ ACS imaging. This is the domain where we expect 
confusion between stars and galaxies to be most likely, based on Figure 8. 
Objects classed as stars based on the CFHT 12K photometry (p^ai < 0.2) 
are the colored points. Objects confirmed as stars in the ACS images are 
blue; those which proved to be galaxies are red. The red points are missed 
by DEEP2 and are thus a concern. The loss is worst (as a fraction of all 
objects) in the dashed rectangle, where high-redshift red sequence galaxies in 
the range z = 0.75-1.0 cross the stellar locus. The loss within that rectangle is 
27% . However, roughly half of galaxies in the affected redshift and parameter 
regime lie in the upper and lower rectangles, which are not lost, giving a final 
estimated loss fraction of ~13.5% (see §8). 



the dimmest magnitude bin are lost, and this loss is uniform 
beyond z = 0.75. This matches the estimate of Willmer et 
al. (2006). A loss of this size is not generally important for 
analyses of total counts along the red sequence but might be 
important (effects as large as 27%) for studies in individual 
restframe color bins within the red sequence that have appar- 
ent B-/? near 2.25. 

• Loss of objects at small separations Because spectra of 
galaxies on the same mask cannot be allowed to overlap, the 
number of slitlets that can be put down in a small area is lim- 
ited for objects that would have similar positions on a mask 
(along its long axis); this affects the sampling rate in dense 
regions. This effect is illustrated in Figure 19, expanded from 
Gerke et al. (2005), which shows that the sampling rate in the 
very densest regions (in terms of nearest neighbor distance on 
the sky) is reduced by almost a factor of two from nominal, 
but that the fraction of sampled galaxies that yield reliable 
redshifts is unaffected. The upturn of the sampling rate at 
small separations is due to the placement of close pairs on a 
single slitlet. 

It should be noted that this effect is much weaker as a func- 
tion of three dimensional density: unlike in nearby galaxy 
surveys like SDSS, most conflicting objects (and most sets 
of nearest-neighbors on the sky) are actually at very different 
redshifts from each other (a consequence both of L* brighten- 
ing with z and the comparatively small contrast in luminosity 
distance between z ^ 0.7 and z ^ 1.4). The effect is already 
substantially smaller in impact than fiber collisions in nearby- 
galaxy surveys, since in DEEP2 objects can conflict with each 
other along only one dimension, and at scales not much larger 
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Figure 19. Rate of spectroscopic observations and redshift success as a func- 
tion of local density of DEEP2 target galaxies on the sky, as measured by 
the third-nearest-neighbor distance, D3. Upper panel: dashed lines show the 
probability that a galaxy meeting the DEEP2 targeting criteria receives a slit- 
let. The heavy solid lines give the total probability that its redshift is reliably 
measured (Q > 3). The light solid lines repeat this with all galaxies in the 
blue corner excluded, as most of those do not yield reliable redshifts. The 
top axis shows the cumulative percentage of z ~ 1 galaxies with D3 up to a 
given value. Vertical lines show the apparent sizes of dense cluster cores at 
two redshifts (co-moving). The sharp increase in targeted fractions at low D3 
arises because close pairs may be placed on a single slitlet. The lower panel 
separately shows the redshift success rate once a galaxy is targeted. The prob- 
ability of observation is reduced by almost a factor of two in regions of high 
density, but redshift success rate is unaffected. 

than an individual galaxy (l"spans ~ 10/i~' kpc comoving at 
z = 1, while the minimum separation between objects along 
the mask direction - i.e. the separation within which they 
conflict with each other - is 3"). 

Nonetheless, the diminished sampling rate in dense regions 
is a large enough effect that it does need to be taken into 
account for certain statistics. For example. Coil, Davis & 
Szapudi (2001) investigated the impact of prototype DEEP2 
target selection algorithms on the £^(rp,Tr) diagram; they 
found that it introduces moderate distortion into the two-point 
redshift-space autocorrelation function at small projected sep- 
arations, reducing the elongation of contours along the "finger 
of God" by about 20%. The data shown in Figure 19 can be 
used for modeling purposes. In all DEEP2 clustering analy- 
ses (e.g. Coil et al. 2006a, 2008) this effect is quantified and 
corrected using mock catalogs that have been passed through 
the DEEP2 target selection code. 

• Multiple galaxies masquerading as single galaxies The fi- 
nal known bias in the DEEP2 target sample is due to the finite 
resolution of the CFHT BRI ground-based images, in which 
separate galaxies may appear blended. This blending can dis- 
tort the magnitudes and colors of DEEP2 galaxies, and it can 
also elevate galaxies that would otherwise be too faint to meet 
the 7?AB = 24.1 magnitude limit above that threshold. Investi- 
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gating these composite objects requires redshift information, 
and so this topic is deferred to §12 after DEEP2 spectroscopy, 
data processing, and redshift measurements have been dis- 
cussed. 

9. SPECTROSCOPIC PROCEDURES 

Instrument and exposure parameters for the DEEP2 survey 
are summarized in Table 2. Each mask is observed for at least 
three 20-minute exposures with the DEIMOS 1200 1 mm"' 
gold-coated grating centered at 7800 A. No dithering is per- 
formed between exposures because sky subtraction is natively 
almost photon-limited (see §10.2), and dithering would waste 
roughly half of the detector real estate. 

DEIMOS' guider camera is rigidly mounted in the slitmask 
coordinate frame with known orientation and offset (Fig- 
ure 16). Locating a guide star at a precomputed pixel loca- 
tion while the mask is at the proper PA therefore succeeds in 
placing all alignment stars within their 4' square slitlets. In a 
direct image of the night sky through the mask, the locations 
of alignment stars and of all slitlets, including the alignment 
boxes, are visible (the latter are illuminated by the sky and 
hence much brighter than the interslit background). 

We obtain the necessary direct alignment images efficiently 
by tilting the grating to zeroth-order, rather than swapping the 
grating for a mirror; the images obtained are then analyzed 
using custom-designed software using a graphical interface. 
The pointing of the telescope and the PA of DEIMOS are 
then adjusted to center all alignment stars within their boxes, 
which in general centers all galaxies in their slitlets to within 
0.1 "(modulo astrometric errors). Slitmask alignment gener- 
ally converges after two direct images and takes roughly 5 
minutes. It is generally checked once using the same proce- 
dure after the first of the three spectroscopic exposures (see 
below). 

The fraction of wavelengths common to all spectra is shown 
in Figure 20; 50% coverage limits are 6500 A to 9100 A. 
The 1 .0" slitwidth used for all targets yields OH skylines with 
a FWHM of ^ 1.3 A, which produces a spectral resolution 
R = AA/A = 5900 at 7800 A. A custom data verification pro- 
gram called quicklook^^ is run during observing. This pro- 
gram performs a coarse reduction of 10% of the slitlets on 
each exposure and computes the cumulative signal-to-noise 
(S/N) for an /? = 23.5 source obtained by combining all ex- 
posures so far. We continue obtaining exposures until a mini- 
mum S/N threshold is met, which enables productive observ- 
ing even under clouds or poor seeing. The quicklook software 
also checks for missing lamps in arcs, the occasional slitmask 
buckling, flexure compensation system misalignments, guid- 
ing errors along the slit, and bad or degrading seeing (based 
on the FWHM of the alignment star spectra). Run on a Sun 
workstation from circa 2002, it is able to evaluate a single 
DEIMOS frame in roughly 10 minutes (compared to our 20 
minute exposure time). 

All slitmasks are observed at airmass below 1.5. Under 
optimal conditions, eight masks can be observed per night, 
yielding 1000-1200 spectra. 

10. REDUCED SPECTROSCOPIC DATA 

10.1. The spec2d Data Reduction Pipeline 

The DEIMOS sf>ec2(i pipeline (Cooper et al. 2012b) is an 
IDL-based code that is modeled on and borrows significant 

"''' See http://deep.berkeley.edu/deep2manual/quicldook.htmI for details. 



C 
o 

o 
« 

u 




6000 7000 8000 9000 
Air Wavelength (A) 



10000 



Figure 20. The fraction of DEEP2 spectra having coverage at a given wave- 
length. The dip around 7800 A is due to the gap between blue and red chips 
in the DEIMOS detector; this gap falls at different wavelengths for different 
spectra due to variation in slitlet position along the wavelength direction. 

code from the SDSS spectral pipeline of David Schlegel and 
Scott Buries (Schlegel & Buries 2012). Most of it was writ- 
ten by Douglas Finkbeiner, Marc Davis, Jeffrey Newman, 
and Michael Cooper, with important contributions from Brian 
Gerke regarding non-local sky subtraction (see below). The 
pipeline operates in five separate stages, each of which pro- 
duces its own output files. The code operates without super- 
vision and usually does an excellent job of wavelength fitting, 
sky subtraction, and object fitting. "*** 

The first step in data reduction is creation of a "planfile" (a 
plain-text file with extension .plan) which may be generated 
automatically after a night's observing and defines the loca- 
tions of all files needed to reduce a given mask's data. The 
planfile may optionally control various parameters of the re- 
duction such as the chips to be reduced, the arc line list to be 
used, etc. The planfile is used by a controlling IDL program, 
DOMASK, which calls separate routines which perform a se- 
ries of reduction steps: 

Step I: The first stage of data reduction is to pro- 
duce calibSlit files (i.e., a set of FITS BINTABLE for- 
mat files whose filename begins with "calibSlit") from 
the flatfield and arc frames taken through a given mask. 
This processing is controlled by the IDL procedure 
DEIMOS_MASK_CALlBRATE. Each of the 8 DEIMOS 
CCDs is analyzed completely independently through all 
stages of the 2-d reductions, allowing trivial parallelization; in 
general, every spectrum will span 2 of these 8 CCDs, which 
we will generally refer to as 'blue' and 'red'; hence for each 
object's spectrum there are two separate calibSlit, spSlit, and 
slit files (q.v. below). The calibSlit files indicate where the 
slitlets fall on the CCD array for each mask and contain flat- 
field information and a 2-d wavelength solution for each slit- 
let. 

To produce the calibSlit files, first the multiple flats are 
first read in, corrected for pixel-to-pixel response variations, 
and processed to reject cosmic -rays. To identify every slit- 

^* The code is publicly available, and further information may be found 
about how it works and how to download it at http://deep.berkeley.edu/DR4 
and http://deep.berkeley.edu/spec2d/. 
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let, edges that are detected in the flats using an unsharp mask 
are compared to a table in the DEIMOS header data that de- 
scribes where each slitlet should appear; a smooth polyno- 
mial is fitted to the edges of eacqh curved spectrum to define 
its upper and lower limits. Each chip's combined flat frame 
is then mapped into rectangular arrays for each slitlet, which 
greatly simplifies further processing, by shifting its edges to 
be parallel, using interpolation (over the spatial direction) to 
accomodate subpixel shifts. Hence, we turn the entire flatfield 
into a set of individual, rectangular flatfields for each slitlet. 
Each slitlet's flatfield is first used to measure the "slit func- 
tion," i.e. the relative throughput (maximum 1) of each slitlet 
as a function of position (i.e. row number) along it.'*^ After 
normalizing by the slit function and any large-scale variation 
in light intensity with wavelength and position along a given 
slitlet, the flats are smoothed on small scales to generate a 2-d 
fringing-correction map with mean 1 . 

Step II: The next step, also controlled by 
DEIMOS_MASK_CALIBRATE, is to solve for the wave- 
length of each pixel in the 2-d data array. First, the entire 
arc (or arcs, if multiple ones are available^") is read in and 
corrected for pixel-to-pixel response variations. Next, a shift 
in the spatial direction between a given chip's arc data and 
the original flatfield, dx, is solved for by cross-correlating 
the intensity summed over the central (in the wavelength 
direction) region of the slit with the light profile predicted for 
the edges found in Step I. Then, individual slitlet arc spectra 
are extracted and rectified in the spatial direction using the 
same algorithms applied for the flatfield, but using slit edges 
shifted from the flat by dx; the individual spectra are then 
corrected for the slit function and fringing using the results of 
Step I. 

Starting with the DEIMOS optical model (produced by 
Drew Phillips), which gives a rough initial guess, we fit for 
the relationship between wavelength and pixel in a spectrum 
in multiple stages, obtaining a wavelength solution that is typ- 
ically accurate at the 0.007 A (--0.02 px) RMS level or bet- 
ter over the entire 2-d spectrum; 0.003 A RMS is commonly 
achieved. The fits used are linear regressions for the coeffi- 
cients of Legendre polynomials of up to fifth order (for wave- 
length as a function of pixel number along the central row of 
the slitlet) or second order (for specifying the local tilt of the 
constant- A locus as function of central pixel number). The 
mean residual about the fit in each row of the slitlet is also 
stored and incorporated into the wavelength solution (i.e., we 
apply a row-by-row shift from the smooth polynomial-based 
solution, e.g. to account for defects or dust along the slit 
which shift the mean wavelength of a row). The requirements 
that accurate OH-line subtraction places on sky subtraction 
are discussed in Section 10.2 and in the Appendix. We gen- 
erally require errors below ^ 0.01 A RMS for good sky sub- 
traction; to be useful in attaining wavelength solutions, indi- 
vidual arclines must contain ^ 10,000 photons in each row. 
The result of Steps I and II is a single calibSUt file containing 
all necessary calibration information per slitlet per chip (i.e., 
each spectrum is split into blue and red files). 

^' Throughout this discussion, we will use 'rows' and 'columns' to refer to 
the rows of a rectified image stored by the pipeline; hence rows run along the 
wavelength direction, with columns along the spatial direction, perpendicular 
to the definitions of rows and columns in raw DEIMOS data. 

'" In general, shifts between successive arc frames are greater than our 
wavelength-solution tolerance for sky subtraction; a single, KrArNeXe arc 
yielded best results. 



Step III: The third stage of the reductions, controlled by the 
procedures DEIMOS_2DREDUCE and DEIMOS_SPSLIT, 
uses the calibSUt files produced in Steps I & II to flatfield and 
rectify the on-sky ("science") data for each slitlet. Response 
correction, slitlet extraction/rectification, and flatfielding is all 
done using the same procedures as for the arc files. All CCDs 
are again treated separately. An inverse-variance image is also 
produced based on the photon and read noise (with compen- 
sations for all multiplicative corrections applied), as well as 
maps of bad pixels (which will have zero inverse variance; 
this map provides flag information on why a given pixel is 
bad) and questionable ones (e.g., those which may be affected 
by imperfect vignetting corrections or cosmic rays, which will 
not have their inverse variance set to zero; again, this map pro- 
vides flag information for each pixel). 

The most important operation in this step is computing the 
proper sky background level (including bright OH lines) at 
every location in a 2-d science spectrum as a function of both 
wavelength and position along the slitlet. First, a constant 
wavelength shift (relative to the arc wavelength solution) for 
each frame is determined using cross-correlation against a 
high-resolution night-sky spectrum. Then, a B-spline model 
for sky intensity as a function of wavelength is fitted to the sky 
regions of the slitlet (i.e., regions where, based on its location 
and size, both contamination from DEEP2 targets and slitlet 
edge effects are expected to be minimal). A key strategem is 
tilting all slitlets by at least 5° with respect to the detector ar- 
ray so that each night sky line is substantially oversampled, 
making the B-spline fit more accurate; see Appendix A. 4. 
The result of Step IV is one spSlit file per slitlet per chip, 
which contains the reduced two-dimensional spectra (with as- 
sociated wavelength, inverse variance, and mask arrays) from 
each science frame for that slitlet and that frame's bspline fit 
in alternating HDUs. 

Step IV: The fourth stage, controlled by SP- 
SLIT_COMBINE and SLITEXAM, combines the separate 
science exposures for each slitlet into one inverse-variance- 
weighted mean, sky-subtracted, cosmic ray-cleaned 2-d 
spectrum (which we will generally refer to below as a 
"Combined" or slit file). First, the B-spline model sky for 
each slitlet is subtracted from its flux. Next, the (generally) 
three 20-minute exposures are combined, weighting them 
according to the inverse variance maps from the previous 
stage. Cosmic ray rejection is also done at this time based 
on the time variability of a given pixel. The result of this 
process is one slit file for each slitlet for each chip, con- 
taining the processed, combined two-dimensional spectrum, 
inverse variance, and wavelength solution, along with various 
diagnostic information (e.g. bad pixel masks). Each row of 
each spectrum has its own wavelength solution; a 2-d map 
of wavelength as a function of position in an image may be 
obtained with the IDE function LAMBDA_EVAL. 

Step V: In the fifth and final stage, 1-d spectra are ex- 
tracted from the Combined 2d spectrum using the routine 
EXTRACT ID (controlled by the procedure DO_EXTRACT). 
Extraction is done using both "optimal" and boxcar (i.e., 
fixed-width sum) extraction techniques (see below), as usual 
processing the blue and red chips separately. As a first step 
in both approaches, each 2d spectrum is rectified in the wave- 
length direction by using whole-pixel shifts (as opposed to in- 
terpolation), making the wavelength scale in each row match 
the central row of the object. 

The optimal extraction (or "home" as it is labeled in the 
FITS headers) loosely follows the algorithm of Home et 
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al. (1986). However, due to the lower signal-to-noise of 
some DEEP2 data, instead of using the light profile evalu- 
ated at a single location along the slit to weight pixels, a 
Gaussian model is used for weighting. First, we measure the 
light profile along each slit by performing an inverse-variance 
weighted average of the flux along each row of the slit, ex- 
cluding regions around bright night sky emission lines. Next, 
an overall offset between the positions of objects along each 
slit predicted from the mask design and their actual observed 
positions is determined by fitting a Gaussian to each peak in 
the individual slits' light profiles and comparing their loca- 
tions to the prediction. We also determine the quantity that 
must be added to or subtracted from (in quadrature) the sizes 
of objects in the CFHT photometry (r^ values) to match the 
observed widths of each Gaussian peak, allowing compensa- 
tion of sizes for differences in seeing. The center of the extrac- 
tion window is determined from either the peak of the spatial 
profile of the object along the slit or (in low-flux cases where 
the peak differs by at least 5 pixels from the prediction) us- 
ing the predicted position with the just-determined shift in the 
spatial direction applied. The extraction center is determined 
independently on the red and the blue sides, ameliorating the 
small shifts between them expected due to differential refrac- 
tion. 

The weighting kernel P used for optimal extraction is then 
a Gaussian with peak at the center of the extraction window 
and a derived from a Gaussian fit to the light profile in the 2-d 
spectrum in high-S/N cases, or the expected seeing-corrected 
size as derived from the photometry in the event of low S/N. In 
our reductions the extracted flux / at a particular wavelength 
is: 

^iPfSi, 

where fi is the flux, P, is the model profile, and 5, is the inverse 
variance at a particular row (at fixed wavelength). Including 
the inverse variance in the weight in this way yielded slightly 
improved results over using the Gaussian profile alone. 

In the boxcar approach, bad pixels are masked (rather than 
interpolated over as in the optimal extraction) and the frac- 
tion of flux missed within these pixels is determined using 
the Gaussian model spatial profile along the slit determined 
for the optimal extraction. The flux from all pixels at a given 
wavelength is then summed, and that sum is divided by the 
fraction of flux which should have come from the included 
rows. Hence, this algorithm differs only in the case of bad 
pixels from a generic boxcar or tophat weighting. For this 
reason, the extraction is called a "boxsprof" extraction in the 
headers of the specld FITS files. The half- width of the ex- 
traction window used is a multiple of the light-profile model 
Gaussian width: 1.75(t (in the optimal case) or 1. 3a (in the 
boxcar case). 

The extractions are saved in a specld FITS file, which con- 
tains the blue and red portions of each object's spectrum for 
each type of extraction, each in its own HDU (i.e. 4 HDUs 
total). Note that in cases where multiple DEEP2 targets are 
placed on the same slitlet, there will be only one set of calib- 
SUt/spSlit/slit files, but separate specld files for each object. 

There is no step in the pipeline data reduction process that 
explicitly flux-calibrates the data using standard star spectra. 
However, DEIMOS is a fairly stable instrument, and a rela- 
tive flux calibration accurate to roughly ±10% over the whole 
spectral range has been determined, based on tests with stan- 
dard stars and the coadded spectra of galaxies on each chip 



(Yan et al., in prep). This error does not include variations 
in absolute throughput due to clouds, seeing, and long-term 
drifts in optics reflectivity (the latter at the 30% level); in ap- 
plications requiring flux calibration, we tie our observed spec- 
tra to the CFHT photometry. Standard star observations are 
taken on a regular basis, but these data are not part of the reg- 
ular reduction process. Throughput plots (for freshly recoated 
optics) are available online.^' 

10.2. Final Data Quality 

Sample raw and reduced spectroscopic images are shown in 
Figures 21 and 22. Figure 21 shows a blowup of a portion of 
a single 20-minute sub-exposure (cosmic rays are clearly vis- 
ible), while Figure 22 shows the fully reduced, CR-cleaned, 
sky-subtracted, and coadded ensemble of three 20-minute ex- 
posures. This particular region represents a typical reduction, 
with sky subtraction that is photon-noise-limited as signaled 
by the random noise pattern under the OH lines (although 
some aliasing is visible in the sky subtraction for slits which 
were not tilted significantly compared to the detector; this 
mask was one of the first observed, and we required a 5 degree 
minimum slit tilt only afterwards). 

Poor sky subtraction happens occasionally, and some ex- 
amples will be shown below. However, quantitative analysis 
indicates that the sky subtraction accuracy is close to photon 
noise-limited even under bright sky lines and under most con- 
ditions. This is illustrated in Figure 23 and Figure 24, which 
show distributions of the quantity 5 /a for sky pixels, where 5 
is the difference between the actual and predicted pixel inten- 
sity and a is the read noise combined with the photon noise 
in that pixel based on the measured sky level at that location 
(inferred from a local median filtering of the inverse-variance 
map). If sky subtraction is noise-limited, this distribution 
should follow a Gaussian with RMS width unity. Significant 
numbers of pixels with bad sky subtraction will depress the 
middle and enhance the tails. 

We consider first Figure 23, which treats pixels between 
the bright sky lines. Two sets of 800 spectra were analyzed, 
one set for which visual inspection confirmed good sky sub- 
traction over the whole spectrum ("normal"), and one set for 
which visual inspection indicated significant systematic sky- 
subtraction errors over at least part of the range (see "bsky" 
comment code in Section 13). The horizontal axis is the ratio 
of the actual flux deviation at that pixel divided by the pre- 
dicted noise due to read noise and photon sky noise at that 
pixel. If the noise is random and is properly predicted, the 
observed curves should be Gaussian with RMS a =\. Both 
sets of spectra do indeed populate the same Gaussian, but with 
RMS width 0.86, not 1.0, indicating that if anything the errors 
have been overestimated. It is reassuring to see that the tails 
of the distribution are not significantly overpopulated. 

Figure 24 repeats this test under more difficult conditions, 
using the same spectra but now with pixels under bright OH 
lines. Small errors in the wavelength scale or in flat-fielding 
will cause the residuals to be large and populate the tails of 
the distribution. The normal, well-subtracted spectra again 
closely follow the same Gaussian curve with a = 0.86; our 
pipeline error model appears to be modestly overconserva- 
tive. The "bsky" spectra look slightly worse but are still re- 
markably similar to the normal spectra here as well as in the 
benign sky regions in Figure 24. The two figures collectively 
demonstrate that our sky subtraction accuracy is very close to 

^ ' http://www.ucolick.org/~ripisc/Go 1 200/tlirough_go 1 200_tilts.gif 
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Figure 21. A section of single raw 2-d 20-min exposure from tlie DEIMOS spectrograph. The [O II] doublet is visible in four objects, but cosmic rays and 
night-sky lines are much brighter. 
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Figure 22. A comparable section of a 2-d image showing the reduced and rectified sum of 3 exposures. The [O II] doublets seen in Fig. 21 stand out strongly, and 
the noise under night sky lines is consistent with random noise except at slit ends (due to their sharper PSF compared to the bulk of the slit). This is an example 
of a typical reduction and subtraction of night sky. Some vertical striping is visible in object regions due to undersampling of the sky spectrum on nearly-vertical 
slits (which were only used in 2002 observations); this correlated noise is absent from the great majority of DEEP2 spectra. 



being photon-limited both on and off the OH Hnes, even in 
spectra that the eye sees as having imperfect sky subtraction. 
Figure 25 shows the major speed gain that has been ob- 
tained by deciding to use high spectral resolution for the 
DEEP2 survey. This figure shows the ratio of exposure times 
needed to obtain the same S/N at each wavelength using a 
spectral resolution of R = 600, which is typical of distant red- 
shift surveys, versus DEEP2's resolution of R ^ 6000. The 
red curve is the raw ratio at each pixel, while the blue curve 
is the red curve smoothed by a 7 A boxcar. The speed gain 
varies rapidly with wavelength depending on whether one is 
on a sky line or not, but the median (unsmoothed) gain from 
7600-8800 A is 1.58. This factor is the savings in exposure 
time. The gain occurs because high resolution confines OH 
photons to a small fraction of the pixels, leaving most of the 



pixels much darker than they otherwise would be. 

Figure 26 illustrates various data reduction problems and 
their associated comment codes. Spectra (a) and (b) show 
regions with poor sky subtraction in the OH lines along the 
slitlets, recognizable by the systematic patches of bright and 
dark sky values under the OH lines. Variable wavelength er- 
rors sometimes cause these residuals to wrap along the slit- 
let, yielding a "barberpole" pattern. Appendix A. 2 derives an 
allowed image-motion criterion of 0.6 px (rms) in the wave- 
length direction in order to keep flat-field fringing-induced er- 
rors in sky subtraction to tolerable levels. DEIMOS's actual 
image motion is two times smaller than this. Thus, the sky- 
subtraction errors seen in spectra (a) and (b) should not be 
due to fringing errors but rather tiny residual sky-line mod- 
eling errors, due either to wrong sky-line profile shapes or to 
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Figure 23. Distribution of tiie quantity A/cr for slcy pixels not under briglit 
sky lines. Here A is the difference between the actual photon count in the 
pixel and the predicted sky count at that location infeired from a local sky 
median, and a is the predicted total noise based on the read noise and pre- 
dicted sky count. Each curve combines the pixels from the sky regions of 800 
2-d galaxy spectra. The blue solid curve is based only on normal spectra; the 
red dashed line is for spectra with poor sky subtraction in at least part of the 
range (designated by the "bsky" comment, see §13). The curves are nearly 
identical, and both closely follow a Gaussian. If the noise were due purely 
to photon statistics plus read noise and were properly estimated, the curves 
should follow a Gaussian distribution with rmd width cr = 1 . The actual dis- 
tributions are nearly Gaussian but with width a = 0.86, indicating that we 
have if anything overestimated the noise in our spectra. Regardless of any 
normalization error, the tails are very nearly Gaussian, indicating that our sky 
subtraction is well behaved. 
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Figure 24. This figure repeats Figure 23 but for sky pixels falling on bright 
OH lines. Wavelength and flat-fielding errors would each produce residuals 
that are larger than the expectation from combined photon and read noise, 
causing an excess at the tails of the distribution. Normal and "bsky" spectra 
are again compared to each other. The normal spectra again closely follow 
a Gaussian curve with a = 0.86. The "bsky" residual distribution is slightly 
worse, but are still remarkably similar to the results for typical spectra on 
sky-lines or for the benign regions depicted in Figure 23. These two figures 
collectively demonstrate that our sky subtraction accuracy is very close to 
being photon-limited both on and oft" the OH lines. 



small wavelength calibration errors. Errors in the latter case 
of only 0.015 px are sufficient to shift the model sky lines 
enough during "de-tilting" to produce noticeable subtraction 
errors. Less than 1% of the total pixel area in the DEEP2 
spectra is noticeably affected by such problems. 

Spectrum (c) in Figure 26 shows an example of burrs or 
dirt along a slitlet, which perturbs the slit function and the 




Figure 25. The speed gain (in exposure time at each wavelength) that is 
obtained by using a spectral resolution 5X/\ = 6000 versus 600. The red 
curve is based on the full-resolution DEEP2 sky spectrum; the blue curve has 
been smoothed with a 7 A boxcar. In regions that are relatively free of sky 
lines (e.g., near 8200 A), the exposure times are equal. In regions near sky 
lines, the high-resolution mode saves a large amount of time by removing 
sky-line photons from those pixels, while in regions directly on top of sky 
lines, the low-resolution mode is faster because the sky brightness is diluted 
there. The blue curve is the best gauge of total savings, as it provides a 
running average of time saved over 7 angstroms. The time savings in regions 
that are dense with sky lines can exceed a factor of 10. 

sky-line PSFs. Spectrum (d) shows scattered light from sky 
lines in the neighboring slitlet that has contaminated the sky 
continuum in this slitlet; such light most often comes from 
the alignment boxes, whose 4"slit widths admit more light. 
The gap between CCDs has eliminated large chunks of the 
spectrum in spectrum (e). Such loss can also be caused by bad 
CCD columns. Spectrum (f) illustrates a common problem in 
which galaxy emission extends into the sky pixels, causing 
the sky spectrum to be too bright, which creates a dim spot 
in the extracted spectrum. Spectrum (g) shows a (rare) bad 
wavelength solution. In spectrum (h), a bad column crosses 
the spectrum, obliterating most of the galaxy, but enough of 
A5007 peeks through to yield a reliable redshift. Spectrum (i) 
illustrates a flux discontinuity across the red and blue CCD 
gap. These cases usually look fine in the 2-d data and are 
caused by a slight mismatch in the extraction window between 
the two CCDs. The final continuum flux is very sensitive to 
the exact width and position of the extraction window. 

10.3. Sample Spectra 

Typical 1-d spectra for eight galaxies are shown in Fig- 
ure 27, arranged by 7? -band magnitude and color. The spec- 
tra have been smoothed with a 15-pixel boxcar. The selected 
galaxies he near z = 0.8 and are sorted into two color bins, 
four from the blue cloud and four from the red sequence. No- 
table features in the 4000 A region are visible, including [O 
ii] A3727, the 4000-A break, Ca H and K, two CN bands near 
3800 A and 4100 A, and the higher-order Balmer lines. The 
lowest four panels illustrate galaxies near the survey magni- 
tude limit (Rab = 24.1). Emission-line redshifts are still quite 
easy at this level for galaxies in the blue cloud (right side), 
but continuum S/N is becoming marginal for red-sequence 
galaxies, as noted in §11. Velocity broadening is detectable in 
most spectra; 66% of all galaxies with redshifts have emission 
linewidths at least 3(j larger than the instrumental signature 
(while 79% of these galaxies have at least one line detectable 
at > 3a). 
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Figure 26. Examples of problematic data and the corresponding comment codes: (a) - (b): Bad modeling of OH lines and resulting systematic residuals. A 
sky spectrum with the wrong profile width (e.g. due to a different effective slit width on sky vs. object regions due to dust or burrs) causes residuals visible at 
both the line center and wings. This is tagged by bsky if problems affect substantial fraction of wavelength range, (c): Dust or burrs along slitlet can cause 
high-order intensity variations in the slit function too strong to be modeled. This is tagged by bsky. (d) Sky continuum errors caused by scattered OH light from 
a neighboring slitlet. The extra-wide alignment-star boxes are the most common sources. If this affects many wavelengths, it gets tagged by bsky. (e) Missing 
columns caused by a set of adjoining bad columns on the CCD. This has caused an inaccurate wavelength solution. Wavelength errors as small as 0.015 px can 
cause positive vs. negative errors across the line, yielding a 'barberpole' effect at OH lines. Slitlets with only partial coverage due to a position overhanging a 
CCD edge often exhibit similar problems. If many wavelengths are affected, the slitlet is useless and the spectrum is assigned quality code Q = -2; if only a small 
region is affected, the spectrum is tagged with bcol and bsky. (f) Emission lines have extended beyond the extraction window into the sky region , so that the 
sky is locally oversubtracted. This case would yield a good redshift but get the iffy comment code to signify that emission-line EWs are incorrect, (g) A slitlet 
with a bad wavelength solution. Gets Q = 2 to indicate need for rereduction. The problem is so rare there is no specific comment code for it. (h) A set of bad 
columns (bcol) cross the spectrum, obliterating most of the galaxy, but a hint of A5007 peeks above and below. If one other feature is clear, this spectrum would 
get 2 = 4 and comments bcol and iffy, (i) Example of a flux discontinuity across the blue-red CCD gap. Signalled by disc. The DEEP2 object numbers for 
these specfi-aare: (a) 11013128, (b) 11033793, (c) 11039367, (d) 11044563, (e) 11043877, (f) 11045853, (g) 13010507, (h) 12029377, (i) 11048040. 



Figure 28 summarizes the mean counts and S/N of the com- 
bined one-hour 1-d spectra as a function of /?-band magnitude. 
Panel (a) shows the median continuum number of photons 
per px for a range of galaxy brightnesses near z ^ 0.8. The 
median value is 14 photons px"' near Rab = 24.1 and scales 
slightly more slowly than the linear relation (shown by the 
dashed line) because brighter galaxies tend to overfill the slit. 
Much of the scatter is due to size variations - larger galax- 
ies at fixed magnitude have lower counts per pixel. Panel 



(b) shows the median inverse variance per px for the same 
spectra, which is nearly flat with brightness, being dominated 
by sky noise (see below). Panel (c) plots the median contin- 
uum S/N per px, which is panel (a) divided by the square root 
of panel (b). The upper, dot-dashed line shows the shape of 
the trend expected if noise is dominated by photon statistics 
from the object, while the lower dashed line shows the shape 
expected if noise is dominated by sky. The actual data are 
bounded by these extremes. 
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Figure 27. Examples of extracted, 1-d spectra from DEEP2, illustrating a selection of red and blue galaxies near z = 0.8. Objects are ordered by brightness 
within the red sequence and blue cloud to illustrate the declining S/N for fainter galaxies. The dip in counts near the centers is due to the gap between the red 
and blue CCDs. The spectra have been smoothed with an inverse variance-weighted boxcar of width 15 pixels. The spectra have not been flux calibrated; the 
vertical bar denotes an intensity of 50 DN/pix. There is a strong bias to preferentially determine redshifts for faint red galaxies that have [O II] emission, as in 
object 41021515 at lower left. 



This figure is useful for predicting the S/N available from 
single and stacked DEEP2 spectra. For example, it says 
that summing together roughly 50 spectra at /? ^ 22.0 will 
yield S/N = 25 per A, which is the recommended minimum 
for absorption-line velocity dispersion and stellar-population 
work (e.g., Schiavon et al. 2007). 

More information about sky counts is given in Table 2. 
Comparing the sky continuum level between the OH lines 
with the read noise tells whether an individual 20-minute ex- 
posure is sky-limited or read-noise limited. The typical con- 
tinuum sky count in the range 7700-8300 A is 13 photons per 
2-d px in 20 minutes while the variance of the read noise is 
only 2.55^ = 6.5e" per px. Thus, read noise is less than sky 
noise but is not negligible. This suggests that some gain in 
continuum S/N could be obtained by using a 600-line grat- 
ing instead of the 1200-line grating. Having explored sev- 
eral tens of thousands of galaxies at very high resolution for 
DEEP2, we have indeed chosen DEIMOS' 600-line grating 
for the DEEP3 survey in EGS, which will capture galaxies 



over a wider redshift range, sample more spectral features, 
and yield somewhat higher continuum S/N, with the tradeoff 
of greater difficulty in measuring linewidths and in identifying 
[O ii] amidst sky lines. 

Figure 29 illustrates various types of information that are 
available from the 2-d spectra. Spectra (a)-(f) illustrate a 
range of galaxy internal motions, from highly ordered and 
resolved rotation curves, as in spectrum (a), to unresolved 
but broadened lines, as in spectrum (f). The quantity 5o.5 = 

\/0.5V}o, + (7^, where Vmt is rotation speed and a is linewidth, 
combines both systematic and random motions into a single 
line-width broadening parameter (Weiner et al. 2006ab) and 
proves to be a remarkably useful and robust internal velocity 
indicator that places the internal speeds of both highly regu- 
lar and highly disturbed galaxies on the same scale (Kassin et 
al. 2008). 

Spectra (g)-(j) show four candidate dual AGNs (Gerke et 
al. 2007b, Comerford et al. 2009). The first two exhibit spa- 
tially separated AGNs, while the third shows a large velocity 
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Figure 28. Mean counts, variance, and S/N of the hour-exposure 1-d spectra 
as a function of S-band magnitude. This figure can be used to estimate the 
S/N for a single spectrum or from stacked spectra as a function of magni- 
tude. Panel (a) shows the median continuum number of photons per pixel for 
a range of galaxy brightnesses near z ~ 0.8. The median value is 14 photons 
pix"' hr"' near Sab = 24.1 and scales slightly more slowly than a linear re- 
lation (shown by the dashed line) because brighter galaxies tend to overfill 
the slit. Panel (b) shows the median variance per pixel for the same spectra, 
predicted from the noise model tested in Figures 23 and 24. The variance 
is nearly flat as a function of magnitude, as it is generally dominated by sky 
noise. Panel (c) plots the median continuum S/N per pixel, which is the 
quantity plotted in panel (a) divided by the square root of the corresponding 
object's value in panel (b). The upper line shows the trend expected if noise 
is dominated by photon statistics from an object, while the lower line shows 
the expected trend if noise is dominated by sky. The actual data are bounded 
by these extremes. 

offset between [O m] and the stellar absorption lines, as would 
be expected for an in-spiralling black hole following a galaxy 
merger Several dozen such offset cases have now been found 
(Comerford et al. 2009). Panels (j)-(l) show examples of close 
pairs that are suitable for measuring satellite galaxy motions. 
Also measurable from the resolved spectra but not shown are 
equivalent widths and line-ratio gradients. 

1 1 . REDSHIFT MEASUREMENTS AND COMPLETENESS 
11.1. The spec Id Redshift Pipeline 

The output of the spec2d pipeline is a set of extracted spec- 
tra, one for each object. Obtaining redshifts and object classi- 
fications from these spectra is the task of a separate software 
pipeline, specld. This IDL code base inherits much from 
the SDSS specBS pipeline produced by David Schlegel, Scott 
Buries, and Douglas Finkbeiner, with modifications and new 
routines by Michael Cooper, Marc Davis, Darren Madgwick, 
Renbin Yan, and Jeffrey Newman. 

The algorithms employed are relatively straightforward; a 



single routine, REDUCEID, guides every step of the pro- 
cess. First, the 1-d extracted spectra are read in and cor- 
rected for telluric absorption (using templates for the A and 
B band absorption provided by C. Steidel scaled as airmass 
to the 0.55 power), for variation in the overall response of 
the DEIMOS spectrograph with wavelength, and for the dif- 
ference between air and vacuum wavelengths (all spectra are 
provided on an air wavelength scale using the lAU standard 
method of Morton (1991). The spectrum is then interpolated 
onto a uniform grid in log[Q(A) with spacing 2 x 10"^ for ef- 
ficiency; in log(wavelength) space, redshifting corresponds 
to a shift in pixel numbers rather than a change in wave- 
length scale. Finally, a version of the spectrum that has been 
median-smoothed using a 2500-pixel window is subtracted 
from it; this renders the redshift fit insensitive to any large- 
scale, coarse features in the spectrum (e.g., due to instrumen- 
tal effects causing variation in throughput or continuum level 
with wavelength). 

Next, the spectrum is separately fit to find the best galaxy, 
QSO, and stellar template matches and redshifts. The figure 
of merit used is x^(z) = ^(fobjeci,i- ftempiatefit,i,zf /erf), where 
E denotes summation over all pixels of the spectrum, f object. i 
is the flux from the object in pixel /, f template f it. i is the flux 
predicted for the best fitting template or linear combination 
of templates at redshift z evaluated at pixel /, and cr, is the 
predicted uncertainty in the flux measurement for pixel / (i.e., 
the inverse variance at that pixel to the -1/2 power). In the 
limit where cr, is a constant, this is equivalent to standard 
cross-correlation techniques (Tonry & Davis 1979),^^ but de- 
weights pixels with poor flux measurements (e.g., due to night 
sky lines). 

For galaxy fits, the best-fitting linear combination of three 
templates is used for each trial redshift: an early-type galaxy 
spectrum (based on the composite luminous red galaxy spec- 
trum of Eisenstein et al. (2003); a model Vega spectrum 
from the Kurucz library^^ (cf. Yoon et al. 2010 and references 
therein) convolved to a velocity dispersion of 84.5 km s"' 
(corresponding to intrinsic velocity dispersion and instrumen- 
tal broadening of 60 km s~' each); and an emission-line-only 
template with line strengths matched to the coadded spec- 
trum of all blue galaxies and with line FWHMs matching 
the predicted width for a DEIMOS observation of an emis- 
sion line having rest-frame velocity dispersion = 60 km s"' 
but observed at the assumed trial redshift. All template spec- 
tra have the continuum subtracted off using the same window 
size and pixelization as for the data. Then we construct the 
function x^(z) with spacing 2 x 10"^ in log[Q(l +z), spanning 
all redshifts from -0.0001 to the largest redshift at which [O ii] 
A3727 will fall within the spectrum (z = 1 .47 for a spectrum 
extending to 9200 A). 

The five deepest minima in x^ are identified, and each one is 
fit with a quadratic over the seven pixels centered at the mini- 
mum. This is used both to obtain an improved estimate of the 
redshift of each minimum and to determine its correspond- 
ing redshift error (as the distance away from the minimum of 
the fit curve where x^ changes by 1). For each of these min- 
ima, an emission-line velocity dispersion is measured by mea- 

^^ Davis and Tonry were able to rapidly do the cross-correlation by shifting 
to Fourier space, but their spectra had an assumed noise that was homoskedas- 
tic, and the strong sky line line at 6300 A was blocked out. In our case at is 
far from constant and we must do the cross-correlation in wavelength space. 

'^ http://kurucz.harvard.edu/stars/vega/vegl000pr25.500000 
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Figure 29. A montage of interesting or unusual 2-d spectra. Tlie liorizontal line corresponds to a restframe velocity of 400 km s"' in all cases. Spectra (a)- 
(f) show rotation curves and/or broadened spectral lines, illustrating the kind of information on internal kinematic motions afforded by our R ^ 6000 spectral 
resolution. Spectra (g) and (h) show galaxies that exhibit two well separated sets of AGN emission lines, possibly indicative of merging black holes. Spectra 
(i) and (j) are examples of significant velocity offsets (^^170 km s"') between the emission velocity and the stellar absorption velocity (marked by the vertical 
white line); these may also be a result of in-spiraling black holes. Spectra (k) and (1) show nearby offset companions (of f ser's), which can be used to measure 
the speeds of satellite motions as in Conroy et al. (2005). The DEEP2 object numbers for these spectra are: (a) 11039144, (b) 11045441, (c), 13033966, (d) 
13058074, (e) 12011735, (f) 12004795, (g) 13025437, (h) 11026433, (i) 11028015, G) 13051207, (k) 11051416, (1) 12025255. 



suring x^ between the spectrum and 40 emission-line tem- 
plates constructed with velocity dispersions a evenly spaced 
from to 360 km s~' (but otherwise utilizing the same as the 
emission-line template used for redshift determination). We 
then determine the best-fit velocity dispersion for a galaxy as 
the minimum of a spline fit to x^i'^) tabulated in 9 km s"' 
bins, and the uncertainty in that dispersion as half the size of 
the range in velocity with x^ no more than 1 larger than the 
minimum. 

For QSO fits, we use a single composite spectrum con- 
structed from the QSO eigenspectra used by specBS, and 
search for the two deepest x^ minima over the redshift range 
0.0033 < z < 5. For stars, we choose the three best (i.e. 
minimum-x^) fits with -0.004 < z < 0.004 using any of the 
SDSS template spectra of stars of type O, AO, F2, GO, Kl, 
MOV, Ml, and LO, as well as the SDSS carbon star template. 
Otherwise, the routines proceed in the same way as for galax- 
ies, but the velocity dispersion is left undetermined. 

At the end of this process, for each DEEP2 spectrum we 
have 10 possible estimates of the redshift plus the best-fitting 
spectral-template match for each one (for galaxies, we specify 
the linear combination of old stellar population, young stellar 
population, and emission lines; for QSOs we only have one 
template; while for stars we specify the spectral type of the 
best match). The results for all objects on a mask are com- 
piled into a single "zresult" file, stored in FITS BINTABLE 
format. We have learned from experience that it is impos- 
sible to choose amongst them automatically, unlike in high- 
signal-to-noise measurements using data with accurate flux 
calibrations (e.g., SDSS; cf. Abazajian et al. 2003). In DEEP2 
data the lowest value of x^ often, but far from always, corre- 
sponds to the correct redshift. Human intervention is required 
to choose amongst these candidate redshifts, or (rarely) to find 



redshifts that the automated pipeline missed. 

11. 2. Visual Redshift Inspection Process: the zspec Tool 

To allow us to check and fit for redshifts interactively, we 
have developed an IDL widget-based program called zspec. 
This program allows the user to examine interactively each of 
the eight candidate redshifts proposed by the specld pipeline, 
or to provide an alternative redshift solution. Using this pro- 
gram, all redshifts from specld have been checked by eye, a 
process typically requiring 1-2 minutes per spectrum. 

The three computer display screens for the zspec IDL 
widget are shown in Figure 30. When a candidate red- 
shift is selected, the zspec control panel (top) displays both 
1 -dimensional and 2-dimensional spectra of small regions 
around any of the six most common absorption or emission 
features seen in DEEP2 spectra. The middle panel shows the 
full 2-d spectrum, which is displayed using the ATV tool de- 
veloped by Aaron Barth (Barth 2001; middle panel). The fi- 
nal window, shown at bottom, displays the full extracted 1-d 
spectrum using the SPLOT tool developed by David Schlegel; 
the user can smooth the 1-d spectrum by varying amounts if 
so desired. The predicted locations of more than 20 common 
absorption and emission features are marked in each of these 
windows. The specld pipeline template that is fitted for a 
given choice of redshift is plotted both in the full view and 
subwindows, and the variance is also plotted in the 1-d spec- 
trum view to aid the user in identifying false features associ- 
ated with sky-subtraction residuals or Poisson noise. Option- 
ally, a smoothed version of the 1 -d spectrum may be plotted 
instead of the raw measurements, and plot ranges may be ad- 
justed arbitrarily to help the user check redshifts. 

In many cases, obvious, resolved multiple emission (e.g., 
the [O ii] doublet or H/3 and [O m] ) and/or absorption (e.g., 
Ca H & K) features will appear in the upper inspection win- 
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Figure 30. Illustration of the zspec widget being used to determine the redshift for a typical star-forming galaxy at z = 0.798. Top: the main control panel, 
showing the four most likely spectral features for the chosen redshift in both 1-d and 2-d. Buttons allow the user to choose a redshift and quality, enter comment 
codes, control the smoothing of the 1-d spectrum, and force fit a redshift if necessary. Middle: the 2-d spectrum for this object (displayed using ATV, Barth 
2001), with tick marks at features for the chosen redshift. An option exists to rectify the slit, removing the apparent tilt. Red tick marks at either side indicate 
the extraction window that was used to produce the 1 -d spectrum. Bottom: the extracted 1 -d spectrum (here smoothed by taking the inverse-variance- weighted 
mean over a rolling 15 pixel window), displayed using SPLOT, with candidate features marked. The aqua spectrum shows the best-fitting linear combination of 
templates, while the red curve shows the variance in the spectrum (rescaled for convenience). 
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dows at the proper positions (compared to the template) for 
some choice of redshift, providing immediate confirmation. 
In other cases, the full spectrum in either the 2-d or 1-d win- 
dows must be examined in order to find rarer emission fea- 
tures or to confirm a redshift using the ensemble of absorption 
lines in the spectrum of an early-type or post-starburst galaxy. 

Confirmation of multiple features (where a resolved [O ii] 
doublet counts as two features) is required to assign a redshift 
to a given object with a reliable redshift quality code (Q = 3 
oiQ = 4; see §11.3). The one-dimensional extracted spectrum 
and the two-dimensional spectrum are both checked to ensure 
that the features are real and not a consequence of perverse 
Poisson noise under a sky line or due to much of the spectrum 
being obscured by a bad column, for instance. Other redshift 
quality codes which the user may assign indicate that an ob- 
ject's spectrum is completely masked by instrumental artifacts 
so that the object was effectively never observed (Q = -2), is 
a star (2 = -1), yields no useful redshift information (Q= 1), 
or may potentially yield redshift information but needs more 
analysis or re-reduction (Q = 2); many, but not all, 2 = 2 cases 
may be resolved with more inspection. Section 13 below for 
a fuller description of the DEEP2 redshift quality codes. 

During the zspec process, the user may also fill in one 
or more comment codes that signify additional information 
about the data; the standard codes are also described in §13. 
These codes are useful, but their application has varied among 
various zspec users, and they are the least homogeneous as- 
pect of the DEEP2 database. Additional description can be 
found in §13. 

In a small minority of cases, none of the 10 x^ -minimum 
redshifts from the specld pipeline are correct for a given 
galaxy, but the true redshift is still measurable with human 
intervention; these constitute less than 1 % of the reliable red- 
shifts in the DEEP2 sample. In those cases, the person in- 
specting the redshifts can use the zspec widget to fit for a red- 
shift using the same algorithms used by specld (described in 
§11.1) but searching only a limited range in redshift around a 
user-specified value and optionally using only a limited set of 
templates and/or spectral features. 

The result of the zspec process for a given mask is a FITS 
BINTABLE "zspec" file, recording the properties of the best 
redshift choice identified by the user (template fit, velocity 
dispersion, etc.), the assigned redshift quality code Q, and the 
user's comment codes, if any. These are compiled into a cata- 
log with one entry for each spectrum that has been checked. 

The results of the zspec process are combiled into a sin- 
gle redshift catalog, labelled the "zcat" . To make this catalog, 
for each mask, we use the zspec results from the most expe- 
rienced zspec user who has examined it. The redshifts (z), 
redshift quality codes (ZQUALITY, or Q for short), and com- 
ments provided by this user, together with all derived param- 
eters determined by the specld pipeline (e.g., the coefficients 
for each template, the measured linewidth, etc.) and selected 
photometric information from the peat files are combined in 
this catalog.^"* There is one entry in the catalog for each ob- 
servation of an object (e.g., if an object was observed on two 
different slitmasks, it is listed in the catalog twice); for conve- 

^'* In a few dozen cases, object numbering changed in the course of the sur- 
vey in regions of the Extended Groth Strip where multiple fields overlap. In 
the zcat file, the integer-format object number, OBJNO, indicates the object 
number for the current best photometry of a given source, whereas the string 
object name, OBJNAME, indicates the object number used in mask design, 
which is used to construct the filename of the corresponding one-dimensional 
spectrum FITS file (i.e. the object's specld file). 



nience, we also distribute a version of the zcat in which only 
a single redshift estimate per object is included^^. 

The redshifts stored in a zspec file are first compared to 
those provided in the corresponding mask's zresult file (i.e., 
the original output from the specld redshift pipeline), which 
contains information on 10 candidate redshifts, selected based 
on reduced y^, for each object. For all objects to which the 
zspec user assigned redshift quality code 2 = 3 or 4 (i.e., ob- 
tained a secure redshift measurement), we compare the red- 
shift identified by the zspecer to the y^ minima compiled in 
the zresult file. If there is a x^ minimum within Aj. = 0.01 
of the user's selected z, then we use the redshift stored in the 
zresult file, rather than the zspec result. This is done to ac- 
commodate objects which went through the zspec inspection 
process before the specld and spec2d code was finalized; in 
such cases, redshifts may change slightly between software 
versions, and we wish to use the best possible estimates of 
object attributes. 

For those objects which were assigned g = 3 or 4 but for 
which A.. > 0.01 for all minima in the zresult file, we attempt 
to re-fit for the redshift, again using the ipecic/redshift-fitting 
code, but restricting the redshift range considered to a limited 
window (half-width = 0.03 in z) centered about the redshift 
chosen in the zspec process. If the re-fit redshift agrees within 
Aj. = 0.01 of the zspec result, then the properties correspond- 
ing to that x^ minimum (its redshift, template coefficients, 
velocity dispersion, etc.) are assigned to that object in the 
zcat. However, if there is no x^ minimum within A, = 0.01, 
the object is assigned redshift quality 2 = 2, as we are unable 
to provide standard pipeline quantities (such as velocity dis- 
persion) in the same way as for a redshift with well-defined 
minimum. There are 138 objects assigned Q = 2 for that rea- 
son. 

Next, we read in a set of ASCII-format tables that provide 
a list of all redshift corrections that have been compiled out- 
side of the usual zspec process (q.v. below); in these cases, 
we override the zspec results. In cases where the the object 
is reassigned a redshift with a quality code of Q = -2, -1, 
1, or 2, we simply update the redshift, redshift quality, and 
comment according to the override information. For those 
assigned a secure redshift (2 = 3 or 4), we again utilize the 
specld redshift-fitting code and search for a x^ minimum in 
a Aj. = ±0.03 redshift window centered about the corrected 
redshift given in the override file. If this fit fails (again a 
Aj. = 0.01 criteria is applied), which occurs for 105 objects, 
we update the redshift and comment, but we set the redshift 
quality to g = 2. When the fit succeeds, we update the zcat 
with the full set of corrected parameters from specld, as well 
as the override-provided redshift quality and comment. 

The measurement of DEEP2 redshifts for difficult cases is 
still a work in progress; in particular, objects with 2 = 2 may 
still be possible to extract redshifts from. A good example 
are the superimposed serendips (supser's; 377 objects); i.e., 
cases where the spectra of two separate galaxies are super- 
imposed in the same spectrum (cf. §13). These have not yet 
been disentangled (and indeed, it is ambiguous which red- 
shift should be assigned in such cases). Another class are the 
spectra that exhibit only one reliable feature (sngls's; 1161 
objects). Many such cases can be resolved using BRI photom- 
etry alone (Kirby et al. 2007), and many more will eventually 

'^ A few objects which were observed on DEEP2 masks (e.g. supernova 
hosts) are not included in the zcat, as they were not selected according to 
DEEP2 target selection criteria. 
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yield to high-accuracy photoz's that are being prepared using 
muhi-band data from m to 8 /im (Huang et al., in prep). These 
revised redshifts are not included in DR4. 

11.2.1. Redshift reinspections 

A few of the objects for which redshift overrides are pro- 
vided resulted from scattered anomalies discovered by the 
DEEP2 team in the course of performing survey science. A 
greater number resulted from two focused efforts: an investi- 
gation of all objects in the Extended Groth Strip with multiple 
zspec inspections or secure redshift codes, conducted before 
DEEP2 Data Release 3; and a reinspection of all 2 = 2 ob- 
jects by experienced zspecers in order to recover additional 
redshifts. 

In preparation for the previous major data release (DR3), 
we undertook a detailed review of the data quality in Field 1 , 
in order to provide the community with the most dependable 
redshift catalog in the EGS region, as well as to ascertain 
the rates of error in our full redshift catalog. This quality 
review was carried out in two stages: 1) A comparison 
and resolution of differences between zspec efforts from dif- 
ferent reviewers; and 2) a review of all g = 3 redshifts in EGS. 

By the summer of 2007 (the time of the review), most 
masks observed in Field 1 had gone through multiple zspec 
checks. The zspec results from all reviewers were combined 
together into a database with 14,509 independent entries, 
where an entry corresponded to a single observation of a 
single target; i.e, multiple observations (on different slitmasks 
or at very different times) of the same target were treated 
as separate entries in the database. Out of all entries in the 
database, 10,763 entries (75%) had gone through the zspec 
process more than once. Several criteria were applied to 
flag cases of possibly erroneous redshifts or redshift quality 
among the multiply-checked targets in the Extended Groth 
Strip. The selected objects (1102 in total) were: 

A) Targets with redshift quality 2 = 3 in all reviews, but 
with differences in assigned redshift greater than 0.001. 
There were 63 such cases, most of which had complex 
spectra that were misidentified by novice reviewers; 

B) Targets assigned 2 = 3 by one or more reviewers and 
2 < 3 by at least one reviewer, the bulk of the sample; or 

C) Targets classified as a 'supser' (cf §13) by any reviewer. 
In certain cases, rare and complex spectra, e.g. from AGNs, 
were assigned this quality class in error. There were roughly 
100 such objects. 

Using zspec, the flagged targets were subsequently re- 
viewed again carefully by expert team members and a final 
judgement was made on their best-fitting redshift and red- 
shift quality. From a total of 1102 flagged targets (10.2% 
of multiply-checked spectra), 544 were verified to have good 
redshifts (2 = 3 or spectroscopically-confirmed stars), 523 
were assigned Q = 2 and 30 were assigned Q = I. The ob- 
jects assigned 2 = 2 were primarily either supser's (which are 
assigned Q = 2 due to their ambiguous redshifts) or objects to 
which zspec users with limited experience had assigned inap- 
propriately high redshift confidence. 

As a result of this process, we developed a procedure to vi- 
sually examine the spectra of all 2 = 3 objects in Field 1 in 
a rapid manner, to ensure that egregious redshift errors would 



be caught. Using custom-purpose software to display multiple 
spectra simultaneously on a common rest-frame wavelength 
scale, we examined approximately 10,000 spectra, 2,600 of 
which had only had a single zspec review at the time. This 
was done in two passes. In the first pass, sections of the 
two-dimensional spectra were examined in regions around 
major emission lines to verify emission line redshifts and 
catch mismatches. Objects which lacked strong emission 
lines were then further reexamined in a second pass using 
their one-dimensional spectra. At the end of both passes, ob- 
jects for which the redshift could not be verified by this quick- 
inspection procedure, 739 in total, were subjected to another 
round of full zspec review by a set of expert team members. 
This guaranteed that every putatively-secure redshift in Field 

1 was reviewed at least twice by a human eye. 

As a result of this process, ^ 250 objects (out of the 10,000 
inspected) were downgraded to 2 = 2, and 8 were downgraded 
to 2 = 1- As a check on the quality of the redshift catalogs 
we distribute, we compared the results of our review with the 
redshift catalog from which DEEP2 Data Release 2 was con- 
structed; it contained the results from 14175 spectra in the 
Extended Groth Strip. We found that out of those objects, 3 1 
galaxies had incorrect redshifts (0.3% of all 2 = 3 redshifts), 
while 249 new, secure redshifts had been obtained for ob- 
jects that previously had Q < 3(largely due to reobservations 
of masks with problematic signal-to-noise), while there were 
493 galaxies with 2 = 3 that we downgraded to 2 < 3 during 
reinspection (i.e., in 3% of cases, the original zspec user was 
judged to be overoptimistic on whether the assigned redshift 
was secure). The results of these redshift reinspections were 
included in the DEEP2 Data Release 3 (DR3) catalog, yield- 
ing a uniform, internally consistent, high quality catalog for 
objects in the Extended Groth Strip. 

As a result of this inspection process, we have now required 
that all DEEP2 masks were inspected by at least one expert 
zspec user, and use their results in the current (Data Release 
4) redshift catalog. Additionally, we have reinspected all ob- 
jects assigned redshift quality Q = 2, searching for reliable 
redshifts that might have been missed in more cursory exam- 
inations. This was done by first having a single expert user 
(Marc Davis) examine all spectra assigned Q = 2 and search 
for an improved redshift estimate. Then, a second expert user 
examined each spectrum assigned a new redshift and judged 
whether they concurred with the result. In cases where both 
experts agreed that a redshift had been recovered, a new red- 
shift, quality code, etc. were assigned and the results incorpo- 
rated into the zcat. As a result of this process, more than 1000 

2 > 3 redshifts were recovered out of more than 5000 spectra 
inspected. 

11.3. Redshift Results 

Redshift results are summarized in Table 7, which gives the 
number of galaxies assigned each quality code in a set of bins 
of apparent magnitude. These are totals over all spectra in 
the zcat redshift catalog; hence, a galaxy observed on mul- 
tiple slitmasks will be counted multiple times, once per ob- 
servation. A total of 52,989 spectra were obtained. 1 .0% are 
so severely compromised by instrumental issues that the ob- 
ject was effectively never observed (Q = -2); 1.8% are stars 
(2 = -1); 16.8% are of such poor quality/signal-to-noise that 
they will likely never yield useful redshifts (2=1; many of 
these are faint blue galaxies beyond our redshift range); 9.8% 
contain information that could be used to determine a redshift 
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Figure 31. Histogram of fractions of DEEP2 spectra in Fields 2, 3, and 4 
yielding redshifts with various quality codes in bins of apparent magnitude. 
Quality code Q = 4 corresponds to redshifts that are >99% secure; Q = 3 
redshifts are assessed to be >95% secure; Q = 2 indicates spectra with a low 
S/N or a known problem as noted in comment codes, but for which a redshift 
may be recoverable; Q = 1 indicates cases with low S/N, that are probably 
not recoverable; Q = indicates objects whose redshifts were never measured 
in zspec; Q = —\ indicates stars; and Q = —2 corresponds to instrumental 
problems so severe that the object was effectively never observed. 

but no definitive determination was made (or no single red- 
shift could be assigned); 11.7% yield a redshift estimated by 
the zspec user to be reliable at the 95% level {Q = 3); and 
58.8% have redshifts estimated to be reliable at the 99.5% 
level {Q = 4). These reliability estimates appear to be accu- 
rate (see below). This information is also shown in Figure 31 
and Figure 32, which show the fraction of spectra with dif- 
ferent quality codes vs. R magnitude in Fields 2, 3, 4 and in 
Field 1 (EGS). The trends are as expected: in Fields 2, 3, 
4, while the fraction with Q = A falls off near the magnitude 
limit, the fraction with (Q=\ rises greatly (owing mainly to 
the onset of the faint blue galaxies), while the fraction with 
2 = 3 is low among the brightest objects but remains fairly 
constant fainter than /?ab = 21. The fraction of potentially re- 
coverable redshifts {Q = 2) is rather constant at about 15 at 
all magnitudes, while the number of stars {Q = -1) and catas- 
trophic instrumental failures {Q = -2) is generally very small. 
The behavior in Field 1 is similar except for the much larger 
fraction of stars at the brightest magnitudes (where there are 
few total targets), which is likely due to differences in target 
selection strategy in that field. 

11.3.1. Reliabilm ofDEEP2 Redshifts 

The accuracy of the redshifts and estimated errors can be as- 
sessed by comparing multiple observations of the same object. 
Duplicate data come from the overlapping areas between the 
top and bottom rows of masks, from the several overlapping 
pointings in Field 4, and from a few other scattered regions 
which have been covered repeatedly (see Figure 1). Due to 
the relatively large fraction of DEEP2 objects observed twice 
(roughly 2.5%), the reliability of our redshift measurements 
can be tested robustly. 

We first consider whether the claimed reliability rates for 
quality codes 3 (95%) and 4 (99.5%) are borne out by real 
data (using only observations outside the EGS, as in that field 
discrepant redshifts were investigated and reconciled). Since 
the claimed error rate for Q = 4 is 0.5%, if we compare two 
independent observations of the same galaxy both of which 



Figure 32. Histogram of fractions of DEEP2 spectra in Field 1 yielding red- 
shifts with various quality codes in bins of apparent magnitude. Aside from 
the larger number of bright stars in EGS due to differences in sample design, 
this figure is similar to that for Fields 2, 3, and 4. 



were assigned Q = 4, they should agree (to within, say, 500 
km s"') 99.0% of the time and disagree 1% of the time. In 
fact, two Q = 4 redshifts disagree only 0.29% of the time (the 
95% upper limit on the mismatch rate is 0.60%), considerably 
better than the claimed quality.^^ Similarly, the mismatch rate 
for g = 3 vs. g = 4 pairs is only 0.68% (upper Umit 1.60%), 
compared to a predicted rate of 5.5%, and the mismatch rate 
for 2 = 3 vs. 2 = 3 pairs is 1 .5% (upper limit 4.4%), compared 
to a predicted rate of 90%. We conclude that the actual catas- 
trophic error rates in zcat redshifts are roughly half as large as 
claimed; i.e., the claimed rates are very conservative. 

To assess the accuracy of our redshift error estimates, we 
consider the set of all objects that have been observed exactly 
twice and for which Q>3 for both observations.^^ Figure 33 
shows histograms of actual velocity differences for pairs of 
objects with g = 3 vs. Q = 3, g = 3 vs. g = 4, and g = 4 
vs. g = 4. The observed RMS difference when comparing 
two redshifts assigned quality 3 is 62 km s~' (with outliers 
removed), compared to a mean predicted RMS difference of 
21.5 km s"' based on the pipeline errors of the same dupli- 
cates, or 21 km s~' based on the mean pipeline g = 3 error 
for the sample as a whole. The observed RMS difference for 
cases where one spectrum was assigned g = 3 while the other 
received g = 4 is 51 km s"', versus 16 km s"' expected based 
on the pipeline errors of the same duplicates, or 17 km s"' 
based on the mean pipeline errors for g = 3 and g = 4 spec- 
tra. Finally, the observed RMS difference for cases where 
both objects received g = 4 is 22 km s"\ compared to a mean 

^^ We note that, when assessing the probability of an event from a Poisson 
process which was observed to occur N times, the typical assumption that the 
uncertainty in the number of events, cr(N) is equal to N is quite inaccurate 
when N ^ S.liv/e assume a flat prior on the true expected number of events, 
fi, at 68% / 95% confidence we can conclude that fi < 1.14/3.00 if we ob- 
serve W = 0; or similarly we obtain upper limits fi < 2.35/4.73 if we observe 
/V = 1, 3.49/6.23 for N = 2, 4.59 / 7.48 for /V = 3, 5.61 / 8.42 for iV = 4, or 
6.52 / 9.032 for W = 5 (as opposed to 68%/95% upper limits of 6.05/8.68 for 
N = 5 for the Gaussian approximation). We present here the best estimate of 
the mismatch rate (=N /Nj^ip, where N is the number of cases of mismatches 
and Niiiip is the number of objects with duplicate observations of the requisite 
qualities), as well as the 95% upper limit on this rate. 

^' Note: all velocity errors and differences in this paper are given in the 
restframe of the object, which means they have been divided by (1 +z). That 
is, Av = c X Az/(1 +2), where Az is the redshift difference. 
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predicted RMS difference of 10 km s"' based on either the 
objects' or the full-sample mean pipeline error. Though small 
number statistics may distort these histograms, it is true in 
general that the observed errors are larger than either set of 
pipeline-predicted errors, and the latter need to be corrected 
upward. From this test, we find that the pipeline Q = 4 RMS 
errors can be made to match the observed Q = 4 RMS errors if 
14 km s"' is added in quadrature to each error estimate, while 
the pipeline Q = 3 RMS errors can be made to match the ob- 
served 2 = 3 RMS errors if 41 km s"' is added in quadrature. 

Figure 34 sheds more light by plotting the individual mea- 
sured velocity differences for each of the above duplicate 
galaxy pairs versus the predicted velocity differences based on 
the pipeline errors of the two galaxies. The predicted pipeline 
error is -\/((o'j,.i)^ + (cr,_2)^)/(l +z), where a^^i is the pipeline 
redshift error of the first observation, CTj. 2 is the pipeline red- 
shift error of the second observation, and z is the mean of the 
two redshift estimates. Points are color-coded by quality code 
(blue for Q = 4vs. Q = 4 pairs and red for 2 = 4 vs. Q = 3). We 
again confirm that the pipeline-estimated errors are too small; 
i.e. the measured differences for galaxies with small pipeline 
errors are larger than would have been predicted. This can be 
cured by adding an extra error in quadrature to the pipeline 
values. The blue lines correspond to the predicted 3- a limits 
if an extra error of 14 km s"' is added in quadrature to the 
pipeline errors of each individual galaxy, while the red lines 
have 43 km s"' total added in quadrature, values based on the 
excess errors determined above. These additions make the 
lines fit the observed 3-0- excursions of the two populations 
fairly well. 

A similar test is shown in Figure 35, which compares er- 
rors for red and blue galaxies rather than dividing up objects 
according to their quality codes (though blue, emission-line 
galaxies predominantly yield Q = 4 redshifts, with a greater 
Q = 3 fraction for red galaxies). By adding an extra 17 km 
s"^ error in quadrature for blue galaxies, or 32 km s"' for 
red galaxies, we can match the observed results reasonably 
well. These tests suggest that an effective method for correct- 
ing DEEP2 pipeline redshift uncertainties to match empirical 

errors is to make the replacement o-y, = A/Cy ,^^,,^^^ + (Tq, where 
(Ty is the predicted empirical error, (Tv^.^^,^^^^ is the error in v pre- 
dicted by the pipeline, and ag = 25 km s"' if 2 = 3, or 11 
km s"' if 2 = 4. We note that the fact that the pipeline un- 
derestimates observed errors is not necessarily a failing of the 
spec2d or specld error model; it could instead reflect miscen- 
tering of objects in slits or the relative rotational velocities of 
the dominant star-forming regions in target galaxies. 

The final way in which we can assess the accuracy of infor- 
mation in the zcat is by checking the consistency with which 
the various quality codes have been applied. Table 8 shows 
agreement among quality codes for duplicate pairs. Star iden- 
tification (2 = -1) seems highly reliable; ^ 90% of objects 
assigned Q = -\ in one observation received the same qual- 
ity code in the other. The assignment of 2 = 4 codes is also 
quite repeatable: galaxies receiving Q = 4m one observation 
received it again 68% of the time and Q = 3 another 19% of 
the time, and hence received a "reliable" rating 87% of the 
time. The Q= I quality code, too, is a fairly consistent indi- 
cator of poor quality, being paired with itself or 2 = 2 78% of 
the time. 

The 2 = 2 codes are more highly scattered, being paired fre- 
quently with other codes either worse or better, as we might 



expect for this marginal category. Nearly half of the Q = 2 
objects observed twice were assigned 2 = 4 on the second 
observation; this could happen because a bad-data problem 
such as bad sky-subtraction or inconsistent continuum levels 
(bsky or bcont, cf. §13) was not present in the 2 = 4 obser- 
vation. 2 = 3 codes also seem to repeat relatively rarely; the 
second observation of a 2 = 3 object is classified with Q = 4 
some 66% of the time. 

A final way to look at these issues is to group all 2 = 1 and 
2 = 2 together to form the "unreliable" class, and all 2 = 3 and 
2 = 4 together as the "reliable" class. Problematic pairs are 
those that have one member in each class. According to Ta- 
ble 8, this happens 14% of the time. Of the discrepant cases, 
roughly half are due to technical problems in one observation 
but not the other: i.e., these are galaxies that would normally 
yield reliable redshifts if the data were good. The remain- 
der appear to be dominated by objects which received 2=1 
or 2 = 2 in one observation due to poor signal-to-noise or a 
single visible emission line. 

11. 4. Target Selection and Redshift Success as a Function of 
Color and Magnitude 

We conclude this section by showing two series of diagrams 
that illustrate the density of objects in the original peat galaxy 
candidate catalog in color or color-magnitude space; then the 
fraction that received slitlets; then the fraction that yielded ei- 
ther 2 = 3 or 2 = 4 redshifts (summed); and finally the fraction 
of 2 = 4 redshifts alone. The first series plots the density of 
galaxies in the apparent R-Ivs. Rab plane. Figure 36 demon- 
strates that relatively blue galaxies dominate a pure R < 24.1 
sample at faint magnitudes. However, only half of them are 
targeted in EGS, as may be seen in Figure 37, due to the fact 
that we de-weight faint nearby galaxies (most of which are 
very blue) in target selection. The same plot for Fields 2,3,4 
(Figure 38) is much more uniform, reflecting the mild 25% 
roll-off of weights over the faintest 0.3 mag in these fields. 

Figure 39 shows the fraction of targeted galaxies (combin- 
ing data from all fields) that yield either 2 = 3 or 2 = 4 red- 
shifts as a function of color and magnitude (i.e. over the same 
plane as Figs. 36 - 38. The steep loss for faint blue galaxies 
is due to the fact that many are beyond our z ^ 1 .4 redshift 
limit. There is also a slight loss in the last 0.3 mag for red 
galaxies, reflecting the difficulty of measuring z's for objects 
with weak emission. Finally, Figure 40 repeats this figure but 
including Q = 4 redshifts only. There is relatively little change 
in the faint blue corner but a large change for faint red galax- 
ies, showing that the latter tend to have 2 = 3 redshifts, in 
keeping with the difficulty of measuring their redshifts. 

The second series of plots is similar to the first, but now 
plots objects in apparent B-R vs. R-I color-color space. 
Figure 41 plots the density of the full target sample in this 
diagram, while Figure 42 shows the target sampling density 
in Field 1 (EGS). Several factors are evident here. First, the 
de-weighting of nearby faint galaxies reduces the sampling 
rate for low-redshift blue galaxies, which may be found in the 
lower-left corner of the z < 0.7 region (above the color cut 
line). Second, the de-weighting of all nearby galaxies in EGS 
reduces the sampling rate of objects to the upper left of the 
color pre-selection boundary compared to those below and to 
the right. Third, the high sampling rate at the outer edges of 
the distribution reflects the high priority given to objects with 
peculiar colors, which are outliers in this diagram. 

Figure ?? repeats this diagram for Fields 2, 3, 4. The main 
features evident are the strong color pre-selection, the slow 
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Figure 33. Histograms of heliocentric recession velocity differences for multiple observations of the same galaxy. Comparisons are separated according to the 
by quality class of each redshift. There are 70 pairs in the left histogram (Q = %/Q = 3), 306 pairs in the middle histogram (2 = 3/2 = 4), and 1 195 pairs in the 
right histogram (2 = 4/Q = 4). The RMS velocity difference for 2 = 3 vs. Q = 3 is 62 km s"' (with outliers removed). The RMS for 2 = 3 vs. 2 = 4 is 51 km s"' 
while for the most common 2 = 4 vs. 2 = 4 case the RMS is 22 km s"' . 
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Figure 34. Velocity differences between observations of doubly observed 
galaxies in which blue points represent cases where both observations yielded 
redshift quality 4, while red points represent cases where one observation 
yielded redshift quality (2) of 3 and the other redshift quality 4. The hori- 
zontal axis is the pipeline redshift error estimate for the velocity difference 
(combining the enors estimated for each observation in quadrature); the verti- 
cal axis is the observed redshift difference. The blue and red solid lines are the 
expected 3-cr errors if the estimated pipeline eiTors for the two measurements 
are combined in quadrature with 20 km s"' or 41 km s"', respectively; the 
conesponding excess for the smaller number of cases which yielded 2=3/3 is 
58 km s"'. The quantities added in quadrature were determined by calculat- 
ing the RMS excess error over the pipeline prediction for 2=4/4 and 2=3/4 
cases, respectively; adding this to the pipeline predictions in quadrature cor- 
responds reasonably well to the 3-(T range of the data points in this plot. The 
pipeline-predicted redshift errors appear to be modestly overoptimistic; it is 
likely the excess eiTor is due to variation between observations (e.g., slightly 
different slit placement on rotating galaxies will lead to variation in the mea- 
sured redshift at this level), rather than due to photon statistics. 



roll-off in the selection function near this boundary due to 
"pre-whitening," and, again, the strong preference given to 
objects with peculiar colors in the outer parts of the distri- 
bution. Redshift success rates are plotted in Figure 43 and 



Figure 35. As in Figure 34 but with samples now split by color rather than 
redshift quality. Blue points represent blue (restframe (U -B) < 1) star- 
forming galaxies whose redshifts are based predominantly on emission lines. 
Red points represent red-sequence galaxies, whose redshifts are based mainly 
on absorption lines. The blue and red solid lines are the 3-cr errors that are 
produced by adding the estimated pipeline error for each measurement in 
quadrature with an extra 17 km s"' or 32 km s"', respectively, to roughly 
match the observed 3-a points of the observed distributions. These quantities 
conespond to the median excess error over the pipeline prediction for blue 
and red galaxies respectively. Pipeline-predicted redshift errors again appear 
to be modestly overoptimistic, by an amount similar to found in Figure 34. 

Figure 44. When considering quality codes Q > 3, the red- 
shift success rate is a relatively flat function for galaxies of 
all types save in the bluest corner of color-color space, where 
z > 1 .4 galaxies may be found. It is worth noting that even in 
the most favorable regions of color space the redshift success 
rate is approximately 90%, not 100%, however Since Q = 4- 
redshifts are predominantly assigned to blue galaxies, it is lit- 
tle surprise that the objects assigned Q = 4 are more localized 
in color space; regions dominated by intrinsically red galaxies 
will generally yield a lower rate of Q = 4 redshifts. 
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Figure 36. Map of the density of DEEP2 galaxy candidate objects in Fields 
1,2, 3, and 4 in apparent R vs. R — I color-magnitude space. Likely stars 
(i.e. objects with p^„i < 0.2) are not included. The color at each position is 
proportional to the square root of the number of objects in the catalog with 
that color and magnitude, on an arbitrary scale from to 100. 



Figure 38. The fraction of the candidate targets shown in Fig. 36 which were 
placed on slitmasks in Fields 2, 3, and 4. 




Figure 37. The fraction of the candidate targets shown in Fig. 36 which were 
placed on slitmasks in Field 1 (EGS). 

The net result of all of these selection and success func- 
tions multiplied together is summarized in Figure 19, which 
illustrates the fraction of galaxies that meet the DEEP2 selec- 
tion criteria that receive slitlets and the total fraction that yield 
Q>3 redshifts, as a function of distance to the third-nearest 
neighbor Overall, 61% of qualifying targets in Fields 2, 3, 
4 are placed on slitlets (59% in Field 1, EGS); and 43% out 
of all qualifying targets are in fact observed and yield Q>3 
redshifts (43% in Field 1, EGS). 




Figure 39. Fraction of the above objects placed on slitmasks that yielded 
either Q = 3 (> 95% confidence) or 2 = 4 (> 99% confidence) redshifts. 
Highly-problematic data (Q = —2) and stars are not counted as either suc- 
cesses or failures, so this is (24 + Q3)/(Ql +Q2 + Q3 + Q4), where Q\ is the 
number of objects assigned quality Q= I, Q2 is the number given Q = 2, etc. 
Hence, this is the fraction of galaxies which were targeted for a useful spec- 
trum for which we obtained a secure redshift. The lack of redshifts for faint, 
very blue galaxies reflects the fact that most of them are beyond our redshift 
limit of z ~ 1 .4. 



12. MULTIPLE GALAXIES 

Section 8 listed blending in the CFHT ground-based BRI 
imaging as a potential source of biases in the DEEP2 target 
sample. Such blends have two possible bad effects. First, 
galaxies that have ostensibly reUably-determined redshifts 
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Figure 40. As Fig. 39 but including only 2 = 4 (> 99% confidence redshifts) 
in the numerator. The lack of high quality redshifts for very faint/blue galax- 
ies continues; in addition there is a deficit of very faint red galaxies in the 
highly-secure category, due to their poor continuum S/N. 



Figure 42. The fraction of the candidate targets shown in Fig. 36 which 
were placed on slitmasks in Field 1 (EGS). 




Figure 41. Map of the density of DEEP2 galaxy candidate objects in Fields 
1, 2, 3, and 4 in apparent B-R vs. R—I color-color space. Objects rejected as 
stars (Pgai < 0.2) are not included. The intensity at each position is propor- 
tional to the square root of the number of objects in the catalog with that color 
and magnitude on an arbitrary scale from to 100. The white line indicates 
the DEEP2 color cut used to select high redshift objects; it corresponds to the 
dot-dashed line in Fig. 9 

may in fact consist of two superposed but separate galaxies, 
in which case the photometry, as well as derived properties 
dependent on continuum strength such as equivalent width, 
reflect the properties of the blend rather than the particular 
galaxy providing a redshift. Second, it may happen that the 
combination of two galaxies, each of which is just below the 




Figure 43. Fraction of the objects placed on slitmasks that yielded either 
2 = 3 or Q = 4 redshifts, over the same plane as Fig. 41. Instmmentally- 
compromised data and stars are not counted as either successes or failures, so 
the plotted quantity is (Q4 + Q3)/(Ql + Q2 + Q3 + Q4). 

Rab = 24.1 magnitude limit, is bright enough that their com- 
bined light brings them above that limit. These pairs may or 
may not contaminate the "reliable" {Q > 2) redshift sample, 
as they may be identified as supsers in the zspec process; 
however, even if they do not, their presence still falsely en- 
larges the target sample and could distort its statistical prop- 
erties. However, the presence of multiple subclumps within a 
single galaxy is not a significant problem from either of these 
standpoints. Similarly, in cases where the light from a pair 
of galaxies with sufficiently high brightness ratio (e.g., 5:1) is 
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Figure 44. As Fig. 43, but counting only redsliift quality 2 = 4 as successful. 



combined, one galaxy dominates and few conclusions would 
be altered. 

The number of superpositions (hereafter called "multi- 
ple galaxies") has been estimated by comparing the low- 
resolution ground-based DEEP2 peat catalog to a high- 
resolution Sextractor catalog based on the HST ACS mosaic 
in EGS by Lotz et al. (2008). This latter catalog consists of 
all objects detected in either the V or / ACS images. Four 
percent of DEEP2 candidate objects with Rab < 24.1 prove 
to have multiple Sextractor matches within a .75 radius, and 
10% have multiple matches within a 1 .0 radius. Inspection 
by eye confirmed that objects within either of these separa- 
tions frequently look single in the CFHT images, depending 
on seeing. 

To understand these objects in more detail, we start with 
the quarter that were targeted for spectroscopy and yielded 
reliable redshifts [Q > 3). Note that this fraction is only half 
as large as the 50% targeting + redshift rate obtained for peat 
objects overall in Figure 19; we will return to that point later. 
The fraction of multiples as a percent of all galaxies with red- 
shifts is plotted as a function of z in Figure 45. 

The number of objects found in the larger aperture vs. the 
smaller remains constant at about 2.5 at all redshifts. This 
number is larger than the ratio of the areas of the two apertures 
(1.8). One possible explanation is that pairs with separations 
below ^ 0.5" are treated as single objects by Sextractor, so 
the ratio of the two numbers is really the ratio of the areas 
of two annuli, the inner one relatively narrow. Alternatively, 
to the degree to which close pairs are likely to be associated 
with each other, it could represent the impact of physical pro- 
cesses: l"corresponds to ll/i~' kpc comoving, or 5.6 /;"' kpc 
physical, at z = 1 . Close pairs could have merged, preventing 
us from seeing them; alternatively, the background galaxy in 
a pair might be hidden from view by dust in the foreground 
object. 

Figure 45 shows a clear increase in the rate of multiple 
galaxies near z ^ 1, apparent with either aperture size. To 
explore this, we have inspected all HST images of the mul- 



Figure 45. The fraction of multiple galaxies vs. redshift for galaxies with 
secure z measurements (Q > 3). Multiples are found by matching DEEP2 
peat galaxies to objects detected in V and / images in the EGS HST mosaic. 
A multiple galaxy is defined as having two or more Sextractor matches within 

a radius of .75 or 1 .0 of the peat position. Below z ~ 1, only about 40% 
of the multiples found here are serious - the remainder are due either to data, 
subclumps within single objects, or very high-brightness-ratio pairs. If these 
are set aside, the remaining fraction of multiples (1-2%) is consistent with the 
expected rate of random superpositions of physically unassociated galaxies 
along the line of sight. The fraction of multiples rises significantly beyond 
z ~ 1 . Though some of the increase is due to subclumps (as seen in the 
HST images), much of the rise must be due to an increase in the number of 
separate but physically associated companion galaxies, which appear to be 
more prevalent at high redshifts. 

tiple galaxies with redshifts by eye to see if Sextractor is re- 
liably identifying the same kind of galaxies at all redshifts. 
It appears that this is not the case. At low redshifts (where 
the multiplicity fractions are low), we find that only 40% of 
the multiple galaxies consist of well-separated, distinct galax- 
ies of comparable brightness - the remainder are either sub- 
clumps within the same galaxy, pairs with high brightness ra- 
tios, or regions with bad ACS data. 

Since these latter cases are not a problem for DEEP2 data, 
we can multiply the fractions of multiples below z ~ 1 in Fig- 
ure 45 by 0.4 to find the true numbers of problematic galaxies. 
We find that ^ 1 % of DEEP2 galaxies at z < 1 have a contam- 
inating galaxy within 0.75", while ~2% have a potential con- 
taminant within I". These numbers are quite close to those 
expected if the companion objects consist of R <^ 25ab mag 
galaxies (i.e., roughly those down to the brightness limit of the 
Sextractor catalog) distributed randomly on the sky. We there- 
fore estimate that at redshifts z < 1, 1-2% of DEEP2 galaxies 
are seriously contaminated by unrelated superimposed fore- 
ground or background galaxies; we further conclude that most 
of these superpositions are not caused by neighboring galax- 
ies, but rather by physically unassociated objects at widely 
separated distances along the line of sight. 

The situation changes above z ^ 1, where the fraction of 
multiples in Figure 45 increases. We would not expect the 
number of random superpositions within a fixed angular sep- 
aration to depend on redshift. As a consequence, the analysis 
we have done for z < 1 galaxies demonstrates that randomly 
superimposed galaxies can make up at most a small fraction of 
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Figure 46. Smoothed contoure .showing the frequency of multiple objects: 

i.e. DEEP2 /jca? objects having multiple companions within 0. .75 radius in 
the Sextractor catalog based on the HST mosaic in EGS (see Figure 45). Here 
we consider all objects in the peat photometric catalog with errors on both 
(B-R) and (R-I) lower than 0.2, not only those for which spectroscopy was 
obtained. The highest peaks lie off the main body of points and are dominated 
by just a few galaxies - their colors are likely compromised by bad photom- 
etry due to multiplicity. The remaining multiples lie mostly in the extreme 
bottom and lower-left of the disitribution, which is dominated by distant blue 
galaxies, most at z > 1 .4. Inspection of the HST images reveals many sub- 
clumps within single objects but an even larger fraction of apparently distinct 
companion galaxies (see §12). Based on their location in color-color space, 
many of these multiples must lie near z ~ 2. 

of the multiple companions at z > 1 . Based on our inspection 
of ACS postage stamp images, the fraction of multiples which 
are in fact single objects with several subclumps doubles near 
z ^ 1 ; this is reasonable, since more distant galaxies tend to 
look clumpier (both because we are observing them at shorter 
rest-frame wavelengths and due to their intrinsic properties; 
cf. Papovich et al. (2003). However, this increase cannot ac- 
count for the full rise, and we are forced to conclude that the 
number of physically associated but distinct companions also 
rises beyond z = 1 . 

These trends can be followed to even greater z by using 
the location of multiple galaxies in the color-color diagram 
as a crude redshift guide. Figure 46 plots the distribution of 
the fraction of peat objects which have multiple components 
within .75 in the HST imaging over color-color space. The 
highest peaks in multiplicity rate lie away from the bulk of 
the galaxy color distribution; they tend to be regions contain- 
ing only one or two galaxies. Inspection suggests that the 
apertures used to measure CFHT 12K photometry for these 
objects were incorrect due to the multiple components, so 
that magnitudes and colors were compromised, yielding pe- 
culiar colors. Setting aside these outliers, the main body of 
multiples cluster strongly in the lower-left part of the dia- 
gram, where z ;SL 1 .5 blue galaxies may be found (cf. Fig- 
ure 13). Evidently most multiples are quite distant, beyond 
the z ~ 1 .4 limit where [O ii] passes beyond our spectral cov- 
erage. This explains why relatively few multiples were both 
targeted for spectroscopy and yielded reliable redshifts (as 
mentioned above). 

The increase in multiplicity rate at higher redshifts inferred 



from Figure 46 agrees with the trend found previously in Fig- 
ure 45, which showed a large increase setting in around z ~ 1. 
Inspection of the HST images of the bluest and presumably 
most distant multiples in Figure 46 shows a similar mix as 
at z ^ 1: a small subset of pairs with large brightness ra- 
tios, a moderate number of galaxies consisting of multiple 
subclumps, and a large fraction of apparently well separated, 
distinct companions which are blended in ground-based pho- 
tometry. 

To sum up, the fraction of DEEP2 galaxies down to Rab = 
24. 1 whose ground-based CFHT photometry is blended due 
to purely random superposition of physically unrelated ob- 
jects along the line of sight is of order 1-2%. These objects 
dominate the population of significantly blended galaxies be- 
low z ^ 1. At higher redshifts, the rate of single galaxies 
which are broken up into multiple clumps in the HST data 
by Sextractor goes up; at the same time, the occurrence of 
physically-associated pairs of galaxies that are well-separated 
in HST images also increases. This population appears to ex- 
tend to higher redshift (z ^ 2. It is intriguing to speculate 
that a different phase of galaxy formation has been detected 
at these higher redshifts, in which a larger fraction of the mass 
accretion occurs in major mergers rather than in minor merg- 
ers or smooth flows. A further inference is that ground-based 
photometry of blue galaxies beyond z ~ 1 is significantly con- 
taminated by companions as much as 5-10% of the time. 

This investigation into multiplicity is qualitative and has 
barely touched the surface of this subject. It does, however, 
give some rough indications for the number of contaminating 
blends expected in ground-based data and their likely effects. 
These issues are likely to be important for deep weak lens- 
ing studies, as close pairs may preferentially align along fil- 
amentary large-scale structure and their photometric redshift 
estimates will be contaminated; photometric redshift calibra- 
tion requirements for next-generation projects such as LSST 
are extremely tight (Ma, Hu & Huterer 2006). Future work 
will benefit from comprehensive spectroscopic redshift sur- 
veys extending beyond z = 1 .4, better photometric redshifts 
for distant galaxies, and near-IR HST images, which will 
establish whether the separate-appearing blue clumps in the 
//ST optical images are in fact distinct galaxies or merely blue 
star- forming clumps within larger potential wells. The new 
CANCELS survey, which includes the Extended Groth Strip 
amongst its fields, will provide an ideal testbed for this study 
(Grogin et al. 2011; Koekemoer et al. 201 1). 

13. DATA TABLES AND DR4 REFERENCE GUIDE 

The two tables in this section contain the main data from 
DEEP2 Data Release 4 (DR4), which is the first release of 
the complete set of spectra from the DEEP2 survey. Future 
data releases are planned as redshifts are further improved 
and additional spectral quantities (e.g., EWs, velocity widths) 
are added. Only the first few lines of each table are pre- 
sented here; the complete tables are available electronically at 
http://deep.berkeley.edu/DR4/. In addition to the tables sum- 
marized herein, DR4 includes a variety of other survey data 
products, such as the design parameters for each DEIMOS 
slitmask, sky-subtracted 1-d and 2-d DEIMOS spectra, and 
the CFHT 12K BRI photometric catalogs from which the tar- 
gets were selected. Below, we provide a brief list of the 
primary DEEP2 DR4 data products along with the relevant 
URLs: 

DEEP2 Redshift Catalog (Table 10): 
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http://deep.berkeley.edu/DR4/zcatalog.html 

Summary of Slitmask Observations (Table 9): 

http://deep.berkeley.edu/DR4/masktable.html 

Slitmask Design Parameters: 

http://deep.berkeley.edu/DR4/maskdesign.html 

1-d and 2-d DEIMOS Spectra: 

http://deep.berkeley.edu/DR4/spectra.html 



2-d Completeness Maps (Figs. 1 and 2): 

http://deep.berkeley.edu/DR4/completeness.html 



CFHT 12K BRI Photometric Catalogs: 

http ://deep .berkeley.edu/DR4/photo .html 



Table 10 lists data on the individual slitmasks, one line per 
observation. In some cases, multiple masks with the same 
mask number but covering different objects were observed 
due to the mask centers being adjusted early in the survey; 
each mask/date combination in the table corresponds to a 
unique slitmask. For each mask, this table records the four- 
digit mask ID number; the observation date; a nominal RA 
and Dec near the center of the mask; the position angle of 
the long axis of this mask; the number of objects from this 
mask in the DR4 zcat data release (and hence the number in 
Table 10 below); the estimated signal-to-noise ratio per pixel 
(S/N) for an R = 23.5 galaxy spectrum on the mask; the esti- 
mated seeing; the redshift completeness obtained on this mask 
for all galaxies targeted; and the redshift completeness ob- 
tained for red galaxies only (R-I > 0.5). 

The first two digits of the mask ID number indicate the 
CFHT pointing, using the same numbering scheme as in Fig- 
ures la and lb, while the last two digits are a position code 
within each pointing.^** Completenesses are based on reli- 
able redshifts only, which are defined as having quality code 
2 = 3 or 4, and seeing is the FWHM measured from stars 
in the slitmask alignment boxes. The S/N estimate used is 
the continuum signal-to-noise per 1-d pixel in the extracted 
spectrum, computed from skyline-free regions near 6900 A 
with 3(7 outliers rejected. The value reported is based on a 
median-median line fit to the S/N vs. magnitude for objects 
with 23.1 <Rab< 24.1, evaluated at/? = 23.5. 

Both redshift completeness fractions are useful indicators 
of mask quality and are plotted along with seeing and S/N 
in Figure 47. We report both as, for galaxies bluer than 
(R-I) = 0.5, the redshift completeness obtained does not cor- 
relate with S/N or seeing; instead, it likely reflects only cos- 

^^ In Fields 2,3,4, the DEEP2 survey as designed consisted of tliree CFHT 
12K pointings per field (see Figure 1. The position code in each pointing 
has values 0-39 in the lower-declination row and 40-79 in the upper row, 
both in order of increasing RA. Vertical masks in the "fishtail" at low RA are 
numbered 80 and 8 1 . This numbering scheme permits a total of 243 masks 
per field, not 120, but not all numbers are filled. Field 1 (EGS) is covered by 
four overlapping pointings divided into eight blocks (see Figure 2. In each 
pointing, there are 16 masks that run perpendicular to the strip, numbered 
0-39 and ordered by increasing RA. There are 14 masks, seven in each block, 
running parallel to the strip, numbered 40-79, and again ordered by increasing 
RA. This allows for a total of 320 masks in EGS, but again not all numbers 
are filled. 



mic variance in the number of bright z > 1 .4 galaxies in a 
given mask region. The trends are as expected, with the frac- 
tions of reliable redshifts increasing strongly with both S/N 
and seeing, which are closely correlated. Probably the sin- 
gle most useful mask-quality indicator is the seeing FWHM, 

which ranges between .4 and 1 .5 with a median value of 

0'.'76. 

Table 10 is the master data table for DEEP2, presenting 
redshifts and associated data for each individual candidate 
galaxy targeted for spectroscopy. Each line corresponds to 
a single observed slitlet spectrum. Multiply-observed objects 
have multiple entries, though a Zhest for each object is provided 
based on information from all spectra (see below). Data are 
given for targeted objects only; many serendipitous objects 
are visible on the 2-d spectra, but they have been uniformly 
analyzed. Their presence is signaled by the comment code 
offser ("offset serendip") in the comment field of the targeted 
galaxy (see below). 

Table 10 is similar to the DEEP2 zcat table at 
http://deep.berkeley.edu/DR4/, but certain confusing or 
rarely-used quantities have been omitted (we describe those 
quantities in §13.0.1), while some other quantities derived 
from mask design files have been added (the dropped quan- 
tities continue to be available in the zcat available online.) 
Table 10 is thus a melding of data in the DEEP2 DR4 spec- 
tral database, the previously-released peat photometry cata- 
logs (Coil et al. 2004b)^'', and data from the slitmask design 
files, collected together for easy reference. 

Objects in Table 10 are ordered by their object number in 
the DEEP2 peat catalogs. The variable names in caps are de- 
fined here or, if pre-existing, agree with published catalogs 
and webpages. Sources of all entries are given in footnotes to 
the table. The columns are as follows: 

Col. (1) - OBJNO: A unique 8-digit DEEP2 object num- 
ber drawn from the peat photometric catalogs. The object 
number format is XXyyyyyy, where XX indicates the CFHT 
field/pointing used to derive the photometry (e.g., "23" for 
DEEP2 Field 2/Pointing 3) and yyyyyy is a unique identifier 
within the pointing. The fields and pointings are shown in 
Figure 1 and Figure 2. 

Col. (2) - RA: Right ascension of the object, in decimal 
degrees. 

Col. (3) - DEC: Declination of the object, in decimal de- 
grees. 

Col. (4) - MAGB: B-band apparent AB magnitude niB from 
the peat CFHT photometry. This is defined as mB = mR + {B- 
R)i", where (B-R)i'iis the apparent color measured through 
a 1 "-radius aperture and niR is the apparent R magnitude (next 
column). Using a small fixed aperture for the color yielded 
the most robust colors and total magnitude corrections from R 
to B and / (Coil et al. 2004b). 

Col. (5) - MAGR: /?-band apparent AB magnitude m« from 
the /5Cflf photometry. This is measured through a circular aper- 
ture whose radius is three times the Gaussian radius, r^, de- 
fined in column (10). If 3r^ < 1", the magnitude is measured 
through a 1 "-radius aperture (Coil et al. 2004b). All magni- 
tudes provided have been corrected for Galactic extinction. 

Col. (6) - MAGI: /-band apparent AB magnitude m/, de- 
fined as nij = mR-(/?-/)i/', where (R-I)irns measured through 
a l"-radius aperture. See column (4). 

Cols. (7), (8), (9) - MAGBERR / MAGRERR / 

^' See also http://deep.berkeley.edu/DR4/photoprimer.html 
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Figure 47. Four panels comparing various indicators of tiie quality of observations in eacli slitmaslc. In all cases the y-axis is the percentage of reliable redshifts 
with > 3. The left two panels evaluate redshift success using all galaxy targets on a mask, while the ones on the right include only "red" galaxies with 
R — I> 0.5. The measure of signal-to-noise (S/N) used here is the average continuum signal-to-noise value per 1-d pixel (throwing out 3(t outliers) computed 
in a window around 6900 A in the extracted 1-d spectra, based upon a fit to S/N vs. magnitude for a given mask evaluated at i? = 23.5. "Seeing" indicates the 
average FWHM of the spectra of alignment stars. The fractional redshift success for the full set of galaxies is lower than when we restrict to "red" galaxies 
because the former sample contains large numbers of blue objects which lie beyond z = 1 .4. The excluded blue galaxies influence the all-galaxy redshift success 
due to variations in their number from mask to mask from clustering; they fail to yield redshifts even in good conditions because of their lack of features in the 
DEEP2 spectral window. As would be expected, for red galaxies our redshift success rate correlates with both signal-to-noise and seeing (which itself strongly 
influences S/N). 



RG: Gaussian radius rg of a circular 2-d Gaus- 



MAGIERR: RMS B-band//?-band//-band magnitude error ( 
sigmaB I <JR I cTi). These error estimates include sky and pho- 
ton noise only. Based on the analyses presented in Coil et 
al. (2004b), there may be additional systematic zeropoint er- 
rors internal to each pointing of order 0.04 mag in B and 0.02 
mag in R and /, but mean zeropoints should differ from point- 
ing to pointing by < 0.01 mag. 

Col. (10)- 
sian fit to the /?-band image, expressed in units of CFHT pix- 
els (one pixel is .207). 

Cols. (11), (12) - EL / PA: Ellipticity ez (EL) and position 
angle (PA) of the object's CFHT /?-band image, derived from 
an analysis of image moments. The ellipticity parameter 62 
is defined as e2 = (1 -b/d) = 2e/(l -I- e), where e is the con- 
ventional eccentricity, and PA is relative to north. Note that 
position angles determined from ground-based images do not 
necessarily match those from high-resolution Hubble images; 
e.g., see Wirth et al. (2004). 

Col. (13) - PGAL: Probability of being a galaxy pgai, based 
on the R-band image. A value pgai = 2 indicates a source with 



unusual BRI colors; pgai = 3 indicates an extended source. All 
sources with pgai = 2 or 3 are treated as galaxies in target se- 
lection. A value of Pgai in the range 0.0-1.0 is the Bayesian 
probability of a compact object's being a galaxy (not a star). 
Objects with pgai < 0.2 are treated as stars and excluded from 
selection for spectroscopy; objects with pgai between 0.2 and 
1 .0 are included with a selection weight proportional to their 
value of Pgai; and those with Pgai>^ are given the same weight 
as an object with pgai = 1. For further explanation see §6.2 and 
Coil et al. (2004b). 

Col. (14) - SFD_EBV: Galactic reddening E(B-V) from 
Schlegel, Finkbeiner & Davis (1998). This value was used to 
correct the photometry for a given object for Galactic extinc- 
tion. 

Col. (15) and (16) - M_B / UB_0: Absolute CFHT 12K 
B-band magnitude Mg and restframe U -B color corrected 
for Galactic reddening, U-Bq , both computed from the 
BRI photometry as described in Willmer et al. (2006). We 
use the subscript here to indicate that the calculations are 
done with z = passbands; no correction for dust internal to 
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DEEP2 galaxies has been applied. Mb is calculated assum- 
ing an LCDM cosmology with h= I; i.e., the listed value is 
MB-5logmh. 

Col. (17) - OBJNO: repeated from Col. (1). 

Col. (18) - MASK: Slitmask number on which the object 
was observed, from Table 10. 

Col. (19) - SLIT: Slitlet number on which the object was 
placed. Slitlets are numbered in order of position along the 
long axis of the mask, starting with Slitlet lying to the far 
left on Figure 16. 

Col. (20) - DATE: UT date of observation in YYYY-MM- 
DD format. 

Col. (21) - MJD: Modified Julian date of observation. 

Col. (22) - SLITRA: RA of sUtlet center 

Col. (23) - SLITDEC: Declination of slitlet center. 

Col. (24) - SLITLEN: Slitlet length in arc seconds. 

Col. (25) - SLITPA: Slitlet PA relative to north. 

Col. (26) - Z: Observed best-fitting redshift z from this 
spectrum, in the geocentric reference frame. This is the value 
obtained from the zspec process. 

Col. (27) - ZBEST: Best heliocentric -reference-frame red- 
shift combining information from all spectra obtained for an 
object, z'lf^,. Uses the highest-Q value or, when there are 
multiple choices of equivalent quality, uses the first value ob- 
served. A value of z'^l, is given only if there is at least one 
redshift with Q>3. 

Col. (28) - ZERR: Redshift error <t, from this observation, 
as explained in §11. 

Col. (29) - ZQUALITY: Redshift quality code Q for this 
observation. Codes are as follows: Q = -2 (data so poor that 
object was effectively never observed; tracked for statistical 
purposes); Q = -l (star); 2=1 (probable galaxy but very low 
S/N; data are not likely ever to yield a redshift); 2 = 2 (objects 
with low S/N and/or for which data are in some way compro- 
mised, but for which a redshift may be obtainable with extra 
effort. The reason for assigning 2 = 2 is generally listed in the 
COMMENT column (below) - all supser's and sngls's auto- 
matically get Q = 2), Q=3 (reliable redshift with probability 
of accuracy > 95%; 2 = 4 (reliable redshift with probability 
of accuracy > 99%). Further details are given in §11, while 
the relevant comment codes are explained below. 

Col. (30) - COMMENT: A comment provided in the course 
of the zspec and catalog compilation process. The comment 
can be used to identify a variety of problems or issues with the 
data (see also http://deep.berkeley.edu/DR3/comments.html). 
As noted in § 1 1 , their application has varied, and they should 
be used with caution; many are a matter of judgment and may 
have been applied differently by different users. See further 
remarks below. 

There are two separate categories of problems that are given 
comments, those that have a high probability of impairing red- 
shift quality and those that may or may not do so. All Q = 2 
redshifts should have at least one comment from the first cat- 
egory; where none is assigned, the "marg" keyword (marginal 
S/N) may be assumed. The first category consists of: 

bsky: Bad sky subtraction. This comment highlights 
badly subtracted sky continuum or sky lines or both. An in- 
correct sky spectrum shape results in incorrectly subtracted 
sky, and consequent errors in the extracted object spectrum. 
Examples are shown in Figure 26. 

bcol : Bad column(s). One or more bad CCD columns 
spoil a significant part of the extracted spectrum and/or the 
neighboring sky region. Note that "columns" here run parallel 



to the dispersion direction; i.e., the name of the comment code 
refers to columns on the detector, rather than columns within 
the extracted 2d spectrum. This is the opposite convention 
from Figure 16, where they run parallel to the slitlet direction. 
The pipeline interpolates across bad columns and manages to 
repair a single bad column fairly well, but three or more ad- 
joining columns seriously degrade data quality, bcont (and 
bs) may also be assigned when a piece of spectrum that is 
curved due to optical distortion falls onto a vignetted region 
columns or onto a gap between CCDs. An example is shown 
in Figure 26. 

bcont : Bad continuum shape. This is a catchall phrase 
expressing the fact that the continuum shape does not look 
right. It is often associated with bsky or bcol. 

bext: Bad extraction window. The extraction window 
used by the pipeline for the spectrum is displaced from the 
actual object, or is too wide or too narrow, or is contaminated 
by a companion. 

edge : Edge. The object is too near the end of the slitlet 
and consequently has poor sky subtraction and/or extraction. 

In the case of these problems (bcol, bcont, bext, 
bs and edge), whether an object is assigned Q = 2 (red- 
shift possibly measurable), Q=\ (redshift not measurable due 
to the object properties) or 2 = -2 (object was effectively not 
observed) will sometimes involve judgement calls, which may 
differ from one zspec user to another 

disc : Discontinuity. There is an unphysical jump in the 
continuum level between the blue and red-side spectra. This 
is typically caused by mismatched extraction windows on the 
red and blue sides or by bad sky subtraction on one side. An 
example is shown in Figure 26. 

marg : Marginal : The S/N is low and the spectrum barely 
meets standards for its assigned quality code. 

sngl : Single. Only one feature is visible (we count a re- 
solved [O ii] doublet as two features). All sngVs, automatically 
get 2 = 2. 

supser: Superimposed serendipitous object 

("serendip"). There are two redshifts measurable from 
the spectrum, and hence it is ambiguous what redshift should 
be assigned to the target. If one or both spectra have only one 
line, supser sngl is used. Either z can be entered into the 
zspec file, and the second z is not necessarily recorded. All 
supser's get 2 = 2. We note that supser's are among the most 
difficult phenomena to spot, and many have probably been 
missed. 

The second category of comments represents conditions 
that can often be present without impairing redshift quality: 

fix: Fix. Some element of the reductions went wrong and 
should be redone. The redshift is not necessarily affected. 

if f y : Iffy. The extracted 1-d spectrum may be adequate 
for determining a redshift (as indicated by the Q code), but 
there is something wrong that makes the observation unreli- 
able as an integrated spectrum of the galaxy, and therefore 
unsuitable for other analyses such as measurement of equiv- 
alent widths, rotation curves, linewidths, or SEDs. Examples 
include: (1) a bad extraction window that gets only part of a 
galaxy or is contaminated by a companion, (2) emission ex- 
tending beyond the extraction window and therefore under- 
summed in the 1-d spectrum (this may also result in oversub- 
tracted sky, as shown in Figure 26.), (3) a severe flux discon- 
tinuity across the blue/red CCD boundary disc. With rare 
exceptions, iffy was used only for objects assigned 2 = 3 
or 2 = 4. It may be further qualified by adding bs, bcont, 
bcol, disc, bext, or edge to explain the nature of the 
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problem. 

offser: Offset serendip. A second spectrum is present 
in the slitlet but is well enough separated or dim enough that 
the extraction window and sky-subtraction of the target object 
are not seriously affected. If only one line is visible, offser 
sngl is used. The redshift of the serendip is not necessarily 
recorded. Serendip spectra may be obtained by download- 
ing the 2-d slitlet spectra at http://deep.berkeley.edu/DR4/. 
Offser's are easier to spot than supser's but, again, some have 
probably been missed. 

fill-gap failed: Identifies those instances where 
a redshift was not able to be identified because the spec2d 
data reduction pipeline failed to construct a robust one- 
dimensional spectrum for the object. In the vast majority of 
cases, this resulted from an exceptionally poor wavelength so- 
lution, by which the blue and red halves of the object spectrum 
(see §7 and §10.1) overlapped in wavelength space rather than 
being separated by a ^ 5 Acorresponding to the gap between 
the two DEIMOS CCDs. 

ZREVISED : Identifies those objects for which the redshift 
was "fixed" based on information obtained outside the first 
zspec analysis of the mask. This includes objects whose red- 
shifts were assigned following the Q = 2 rechecking process; 
for more details see §11. 

ZREVISED: zfix failed: Identifies those cases 
where a redshift revision was attempted (see ZREVI SED), but 
failed. This failure generally resulted from the inability of the 
specld redshift pipeline to find a suitable fit that matched the 
revised redshift value (e.g., within Az = 0.01). 

ZMATCH-FAIL: Identifies those cases where a redshift 
that was previously identified as 2 = 3,4 was refit, but the 
specld redshift pipeline failed to yield a good fit that matched 
the previous redshift value (e.g., within Az = 0.01). One mo- 
tivation for this refitting process was to account for changes 
in the spectral reductions after a given slitmask had been vi- 
sually inspected using the zspec tool. 

13.0.1. Additional Parameters in Online DR4 Catalogs 

In addition to the quantities described above, the zcat red- 
shift catalog distributed at the DR4 website includes several 
additional parameters which are of more limited utility or as 
yet only partially vetted. These are: 

OB JNAME : The DEEP2 object number in string format, 
used to define filenames. This can differ from OBJNO in 
overlap regions between CFHT pointings, where OBJNAME 
was defined from the object number used for object selection, 
whereas OBJNO indicates the instance of an object with the 
best photometry. 

CLASS : Indicates the class of template which yielded the 
best redshift fit in the zspec process; possible values are 
'STAR', 'GALAXY', or [broad-line] 'AGN' (see §11.1). It 
should be noted that this only indicates the best-fitting tem- 
plate; many AGN (both broad- and narrow-line) in the DEEP2 
sample will have 'GALAXY' as their class. 

SUBCLASS: For objects fit with stellar templates, indi- 
cates the spectral type of the best-fitting template. Note that 
only a sparse set of spectral types was used in the fits, so this 
does not necessarily match the actual spectral type of the tar- 
get. 

RCHI2 : Reduced chi-squared for the best fit template 
compared to the observed spectrum. Based on results of the 
zspec process, this value is not necessarily indicative that a 
given redshift is good or not. 



DOF : The number of degrees of freedom used to determine 
RCHI2. 

VDISP: The estimated velocity dispersion of the given 
galaxy (in km s"'), determined as described in §11.1. VDISP 
values appear to be a good estimate of velocity dispersion, but 
may be subject to small systematics, especially at low VDISP. 

VDISPERR: The estimated uncertainty in the VDISP 
value for a given galaxy. 

14. TWO-DIMENSIONAL SELECTION FUNCTIONS 

As part of DEEP2 Data Release 4, we are releasing a set 
of files describing the two-dimensional selection function of 
the DEEP2 Galaxy Redshift Survey: i.e., the probability (as 
a function of position on the sky) that an object meeting the 
DEEP2 target criteria (magnitude limit and color cut, as ap- 
plicable) is selected for observation and successfully yielded 
a redshift. The calculation of these selection function (or, 
as they are sometimes labelled, "window function") maps has 
been described in Coil et al. (2008) and references therein; we 
review here. 

The two-dimensional selection function maps take into ac- 
count the actual placement and geometry of the overlap- 
ping slit masks used for DEEP2, as well as vignetting in 
the DEIMOS camera and the locations of gaps between the 
DEIMOS CCDs. In producing these maps, we account for the 
multipass nature of DEEP2 targeting and the fact that some 
overlapping DEEP2 slitmasks were designed simultaneously 
and some were not (as the mask design code will not place an 
object on multiple masks designed at the same time and in the 
same row of masks; cf. §7) in computing selection probabili- 
ties. The redshift completeness for red (R-I > 0.5) galaxies 
is used to determine the probability a targeted object yields 
a redshift; the redshift success rate for bluer objects is both 
low and not correlated with observing statistics such as see- 
ing (rather it primarily reflects the density of z ^ 2 galaxies at 
a given location on the sky), so the red galaxy completeness 
provides the best estimate of the probability that a targeted 
galaxy at a redshift where DEEP2 can obtain a z will actually 
yield a reliable redshift measurement. 

In these maps, we also mask out all regions where the pho- 
tometric data are affected by either saturated stars or CCD de- 
fects and hence no galaxies were targeted. The region where 
these maps are nonzero provides a geometrical outline of the 
regions where DEEP2 targeted on the sky, while the actual 
value at a given position represents the probability that we 
both targeted a galaxy (which depends on how many masks 
overlap a given area, whether the object is in the first or 
second-pass region of those masks, and whether they were 
designed simultaneously) and then successfully measured its 
redshift. Note that the provided files do not include pointing 
14 in the northernmost region of the EGS. Due to poorer pho- 
tometry, the mask design differed in that pointing from the 
rest of the survey, such that the spatial selection is not uni- 
form between it and the rest of DEEP2. No DEEP2 clustering 
measurements use data from that pointing. In an upcoming 
paper (Matthews et al. 2012) we will provide improved pho- 
tometry for this region based on imaging obtained as part of 
the CFHT Legacy Survey (Gwyn 2012). 

The selection function maps are distributed at the DEEP2 
DR4 website in the form of FITS-format images, with World 
Coordinate System headers describing the mapping from right 
ascension and declination to pixel, and the value at a given 
pixel (ranging from to 1) being the combined selection and 
redshift success probability for a DEEP2 target galaxy at that 
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position. Outside of the EGS, where the mask-making was 
done independently for each CFHT pointing, we provide one 
file per pointing. For the EGS, we provide one file describing 
the selection function over the entire field. In order to pro- 
duce random catalogs for calculating correlation functions, 
one should use the value of the selection function at a given 
position as the probability of keeping an object (between 
and 1) at that location within the catalog. We note that for 
some measurements, it may also be necessary to correct for 
the dependence of selection probability on source density at 
small scales (cf. §8); in DEEP2 science papers, this has been 
done with mock galaxy catalogs (e.g. Yan, White, & Coil 
2004). 

15. TRENDS WITH REDSHIFT 

We conclude this overview of the DEEP2 survey by illus- 
trating certain major trends in the data as a function of red- 
shift, summarizing properties of the DEEP2 dataset. 

Figure 48 shows redshift histograms for each of the four 
fields, using only reliable iQ = 3 and 2 = 4) redshifts. These 
figures illustrate graphically that, fractionally, the contribution 
of 2 = 3 (> 95% reliable) redshifts is rather small; the great 
majority (83%) of secure redshifts have Q = 4- (> 99% reli- 
able). Field 1 (EGS) is well sampled at all redshifts, whereas 
our BRI color pre-selection has efficiently eliminated fore- 
ground galaxies in Fields 2, 3, and 4. As expected, strong 
peaks due to large-scale structure are evident, which differ in 
detail from field to field, demonstrating the need for statisti- 
cally independent samples to beat down cosmic variance. 

Figure 48 also sheds light on whether a significant number 
of of redshifts are lost when important spectral features fall 
on atmospheric absorption bands (e.g. the A band of O2 at 
7620 A). The [O 11] A3727 emission line should fall on the A 
band for objects at z = 1.045. There is no depression in the 
diagrams at that redshift, indicating that any loss is small. 

Figure 49 plots apparent magnitudes and colors versus red- 
shift. The top row is for Field 1 (EGS); the bottom row 
combines all objects in Fields 2-4. The most striking fea- 
ture of these diagrams is the clear bimodality visible in the 
color plots. The lower, densely populated sequence consists 
of blue-cloud galaxies, while the more thinly populated upper 
trace is the red sequence. The latter fades out near z ^ 1, due 
partly to a real reduction in numbers there (Bell et al. 2004, 
Willmer et al. 2006, Faber et al. 2007) and partly because the 
survey begins to lose red galaxies rapidly at that redshift. This 
is a combined effect of such objects having difficulty meeting 
the DEEP2 magnitude limit (as the CFHT R band is shortward 
of the 4000 A break past z = 1) and poorer redshift success for 
faint red objects without emission lines (cf. the discussion of 
color-magnitude diagrams below). 

Also notable are the two small islands of very blue, z ^ 0. 15 
galaxies visible in Fields 2, 3, and 4. These reside in the 
t)lue corner of the color-color diagram where the colors of 
low-redshift, high star-formation rate galaxies overlap those 
of more distant blue galaxies (cf. Figure 13); as a result, a 
BRI color pre-selection cannot reject them cleanly without 
also rejecting real high-redshift galaxies. Some of these are 
exceedingly low in luminosity, with Mb as low as -14 (see be- 
low), and likely present an interesting population of objects 
for follow-up studies. 

The vertical striping visible in these figures (especially for 
Field 1) is a consequence of cosmic variance: redshifts where 
the density in our fields are higher than average will have more 
objects of all colors and magnitudes, and those with low den- 



sity will have few. The stripes correspond directly to the fea- 
tures visible in the redshift histograms in Fig. 48. In Fields 2- 
4, the impact of large-scale structure is much reduced, as we 
are effectively averaging the properties of three statistically 
independent regions of the Universe. A significant dearth 
of galaxies near z = 1.15 is nevertheless still visible, despite 
combining roughly 22,000 galaxies in three widely separated 
fields, illustrating how difficult it is to reduce cosmic variance 
to a small error at every redshift (at z = 1.15, [O 11] A3727 
overlaps with one of the cleanest regions in the night sky spec- 
trum, eliminating sky subtraction errors as a possible culprit). 
The impact of cosmic variance increases the smaller the red- 
shift bins (fe) considered. As a result, we find no detectable 
correlation between the DEEP2 redshift histogram and the 
DEIMOS night sky spectrum (evaluated either at the central 
wavelength of [O 11] or at the wavelength of either doublet 
component). Variations in the DEEP2 redshift success rate 
on the scale of the width of a night skyline are completely 
swamped by the large cosmic variance on those scales, if they 
exist at all. 

Figure 50 plots absolute B magnitude and restframe (U - 
B)o color versus redshift. These are calculated from the 
CFHT 12K BRI photometry using the methods of Willmer 
et al. (2006); we use the subscript to indicate that rest- 
frame passbands are used (no correction for dust internal to 
DEEP2 galaxies has been applied). The color bimodality and 
the low-redshift blue galaxies are again evident. The most in- 
teresting aspect of this figure is the curved locus in (U -B)q 
vs. z traced out by blue-cloud galaxies. The solid line shows 
the median color of blue cloud galaxies (defined as having 
{{U-B)o < 1.0) as a function of z, while the dotted and dashed 
lines shows quintiles of the distribution. Figure 5 1 suggests 
that this curvature is not due to errors in the K-correction pro- 
cedure for restframe (U-B\) but is instead caused by the inter- 
action of the fixed /?ab = 24. 1 magnitude limit with the color- 
magnitude distribution of galaxies in the blue cloud (see be- 
low), combined with the impact of any color evolution (Blan- 
ton 2006). 

Figure 5 1 plots restframe color-magnitude diagrams binned 
by redshift. The well-known color trends within both the red 
sequence and blue cloud are clearly evident, with brighter 
galaxies being redder in both cases. Straight lines approxi- 
mate the boundaries imposed by the fixed ^ab = 24. 1 survey 
limit at different redshifts. In each panel, the heavy solid line 
shows this limit for the far side of the indicated redshift bin 
while the dashed line is the corresponding limit for the near 
side of that bin. Light dotted lines repeat the lines from all 
other bins to guide the eye. 

The aggregate trends of (L'^ - B)q with z seen in Figure 50 
can be explained by what is seen here. At the highest red- 
shifts, the survey tends to include bluer blue-cloud galaxies 
but lose redder ones (at fixed Mb), as the limit lines become 
highly slanted in restframe color-absolute magnitude space. 
The mean color of the most distant detected blue-cloud galax- 
ies is therefore by necessity bluer than for lower-distance sam- 
ples. At lower z, the R = 24.1 limit line is more vertical, per- 
mitting redder blue-cloud galaxies to enter the survey across a 
wider range in Mb- At even lower redshifts, the survey is able 
to pick up abundant, very faint blue galaxies, which causes the 
median color of blue-cloud galaxies selected to fall again. We 
have not modeled these effects quantitatively, but it is clear 
that the most obvious trends in color plots like Figure 50 need 
careful interpretation. 
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Figure 48. The histogram of reliable redshifts in each field. The histogram including only Q - 
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16. CONCLUSIONS AND ACKNOWLEDGEMENTS 

This paper has presented DEEP2 Data Release 4, the first 
spectroscopic database resulting from a complete DEEP2 
Galaxy Redshift Survey sample, containing 52,989 spectra 
and 38,348 reliable redshift measurements for objects as faint 
as /?AB = 24.1. The information in this paper is intended as 
a handbook for users of the survey, as well as for those us- 
ing the DEEP2 data reduction pipelines to reduce DEIMOS 
data. Both raw and reduced spectra, photometry, and all red- 
shifts are available for download from the Data Release 4 
website (http://deep.berkeley.edu/DR4/) as well as other web- 
sites listed in §10. We do not expect this to be the final data 
release from the DEEP2 survey. Future releases may include 
improved photometry, revised redshift estimates (e.g. from 
incorporating photometric redshift information for 2 = 2 ob- 
jects), improved data reductions and extractions, and/or addi- 
tional derived parameter measurements not included here. 

The discoveries from DEEP2 have led to the development 
of the DEEP3 Galaxy Redshift Survey (Cooper et al. 2011, 
2012a). This project has obtained more than 8,000 spectra, 
primarily of objects in the Extended Groth Strip. The primary 
goals of DEEP3 were to obtain redshifts for three different 
samples: 



• Objects of interest identified through AEGIS multi- 
wavelength observations; 

• ^AB < 24.1 galaxies that were not previously targeted 
by DEEP2 to avoid slit overlaps, enhancing the sample 
of objects which can be used to trace galaxy environ- 
ments in the EGS, and providing an enlarged DEEP2- 
like sample of galaxies all of which have deep multi- 
wavelength imaging; and 

• Faint, highly star-forming galaxies down to /?ab = 25.5, 
allowing us to test whether they continue to follow the 
same scaling relations as brighter galaxies. 

In order to maximize complementarity to DEEP2 data, 
DEEP3 observations were conducted with a 600-line/mm 
(R ^ 3000) grating on DEIMOS, with spectra typically cov- 
ering the wavelength range from 4600 A to 9800 A. Hence, 
DEEP3 provides worse kinematic information than DEEP2, 
but covers more spectral lines in a given galaxy. We expect 
the first data release from DEEP3 to occur in the next year, 
and that the legacy of DEEP2 will continue into the future. 
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Figure 49. Apparent magnitudes and colors of galaxies versus redshift for the completed DEEP2 survey. Only galaxies with reliable (i.e. Q = 3 and Q = 4) 
redshifts are shown. The upper row shows objects from the EGS, while the lower row consists of objects from Fields 2, 3, and 4. Color bimodality is clearly 
visible even in apparent colors. 
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Figure SO. Absolute magnitudes and colors of galaxies versus redshift for the completed DEEP2 survey. Galaxies are the same as in Figure 49. The upper row 
represents EGS, the lower row Fields 2, 3, and 4. Lines in the right panels represent the 1- and 2-a ranges for blue-cloud galaxies only. We argue that most of the 
complicated trend evident in the plot of restframe I] —B color {(U —B)q) vs. redshift is created by the apparent S-band magnitude selection limit beating against 
the intrinsic distribution of blue-cloud galaxies in the color-magnitude diagram at each redshift. See discussion in § 1 5 and Figure 5 1 . 
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APPENDIX 
A. APPENDIX 

A. 1 . Sky Subtraction 

This appendix describes the steps we have taken to allow photon-limited subtraction of night sky emission, even in the limit of 
bright OH lines. There are two major requirements necessary to achieve this goal. 

First, we must ensure that the photometric response of all pixels remains constant to high accuracy between afternoon flat- 
fielding and evening observations. If this response varies, the flux calibration between object and sky pixels will differ and 
sky subtraction will be imperfect. As explained below, the response in each pixel varies periodically with wavelength at red 
wavelengths, a phenomenon commonly known as "fringing." Hence, in order to ensure accurate sky subtraction we must keep 
the wavelength of light falling on a given pixel constant to high accuracy. We can attain this by keeping the spectral image of 
a given mask accurately centered on exactly the same pixels at all temperatures and spectrograph position angles. Since natural 
flexure in large spectrographs is larger than this tolerance, we have incorporated a "flexure compensation system" in DEIMOS 
to keep the images accurately centered. This system has essentially solved the problem of fringing by ensuring that the fringe 
pattern is identical in flat-field and science frames. 

The second key requirement is that we be able to employ brightness information in the sky spectrum offset from a target object 
to infer sky brightness at the same wavelength on the object spectrum accurately. This requires a highly accurate determination 
of the wavelength of light falling on each pixel compared to its neighbors. The challenge is very high because of the extremely 
steep gradient in brightness on the shoulders of the OH sky lines. For the 1200-line grating used for DEEP2, the slit image has a 
FWHM of 3.9 px, but the native FWHM of the spectrograph is only 2.3 px (Table 2). It is this latter quantity that sets the gradient 
in the measured sky line flux. Since the native PSF is quite sharp, the brightness falls off extremely fast, and shifts of only 0.015 
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Figure 51. Color-magnitude diagrams for all fields, binned by redshift. The straight lines approximate how the apparent Sab = 24. 1 limit of the survey translates 
into restframe color and absolute magnitude limits at each redshift. The solid line in each panel is the completeness limit for the far side of (i.e. highest redshift 
within) that bin and represents the absolute completeness limit for that bin; the dashed line indicates the limit at the near side of a bin. Dotted lines repeat the 
limits from other redshift bins, using the same colors in all panels. For blue-cloud objects, the magnitude limit favors rather blue galaxies at the highest redshifts, 
redder galaxies at intermediate redshifts, and finally bluer galaxies again at low redshifts. These trends are a major reason why the median (U —B)q line is curved 
for blue-cloud galaxies in Figure 50. 

px between the sky spectrum and the object spectrum produce detectable systematic residuals. We therefore need a wavelength 
solution that is better than this over all pixels of the 2-d spectrum, plus a model sky spectrum that tracks the steep gradients with 
wavelength in the sky spectrum to high accuracy. 

The next section discusses the basic physics of fringing and derives an image-motion tolerance needed to keep flat-field errors 
due to fringing negligible. We then describe the flexure compensation system and its performance in light of this tolerance. 
Finally, we describe the methods used to determine wavelength solutions, the b-spline modeling of the sky spectrum, and the 
sky-subtraction process. 

A.2. Fringing Physics and Image Motion Tolerance 
The peak OH sky-line brightness in a one-hour exposure is roughly 3000 photons per pixel with the 1200-line grating. Hence, 
photon-limited sky subtraction requires a flatfielding accuracy much better than l/-\/3000, or 1.8%. Lower-dispersion gratings 
have higher peak counts and need even higher accuracy, as do multiple stacked exposures. In practice we chose the target of 
±0.25% for DEIMOS. A major contributor to flat-fielding errors is CCD "fringing," as explained above. The CCD silicon layer 
is a resonant cavity that supports constructive and destructive interference through internal light reflections (Newton's rings), and 
the pixel sensitivity varies according to whether the interference is constructive or destructive. Tests on the Kast spectrograph 
detector at Lick and the LRIS spectrograph detector at Keck show that the fractional sensitivity variation, F, is sinusoidal as a 
function of the wavelength shift in A, s. Hence, it can be described by the function F{s) = A sin =^, where A is the fractional 
fringe amplitude on a pixel and P is the phase wrap period in A. The latter is given hy P= ^, where A^ is the interference order 
N = ^, n is the refractive index of silicon (3.6 at these wavelengths), A is the wavelength of light falling on a given pixel, and 

t is the detector thickness, which for DEIMOS is 40/im. Thus, A^ ^ 300 at 9000 A, so P ~' 30 A. The interference is extremely 
high-order because the effective thickness of the detector is much larger than the wavelength of light. 

The greatest sensitivity to wavelength shifts occurs where the derivative of F is at its maximum value, ^ = ^^. The fringe 
amplitude A in the thick CCDs used for DEIMOS is only ±2% at 9000 A, which minimizes fringing. Keeping changes in F 
below 0.25% therefore requires keeping s ^ 0.6 A, corresponding to an image shift of roughly L8 pixels with the 1200-line 
grating. However, the tolerance is tighter for lower-dispersion gratings because their dispersion is lower and because the OH sky 
brightness per pixel is higher Coadding many spectra for the same object also requires better flat-fielding for flat-field errors 
to remain negligible compared to uncertainties in sky subtraction from photon statistics. Allowing for the use of other gratings 
and/or longer exposure times motivated a tighter design tolerance of 0.6 pixels RMS flexure in the wavelength direction. Errors 
in the direction along the slit are less problematic; their design tolerance was set at 1.0 pixels RMS. 
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A.3. The DEIMOS Flexure Compensation System 

In order to ensure such a high degree of stabiHty, a closed-loop flexure compensation system (FCS) was built into the basic 
design of the DEIMOS spectrograph. The system contains two movable elements: the grating tilt mechanism, which is used to 
steer the image in the dispersion direction, and a stage in the dewar, which moves the CCD perpendicular to the dispersion. Two 
pairs of optical fibers pipe light into the spectrograph at opposite ends of the slitmasks. The light in spectroscopic mode comes 
from a ThAr arc lamp (an LED is used for direct imaging mode), which produces arc -line images (spots) on the two separate 
flexure-compensation CCD detectors, which are shown in Figure 16). 

Despite the small size of these detectors (600 x 1200 15-/i pixels), the ThAr spectrum is so rich that spots are available with 
nearly all possible gratings and tilts. The FCS spectra are determined from ^ 1-3 second exposures on the FCS CCDs, and then 
the locations of the ThAr spots are processed. Any shifts needed to recenter the spots in their proper locations are then fed back 
into the grating tilt and CCD actuator mechanisms. The loop takes roughly 10 sec to complete. A sensing system turns FCS 
tracking off and on automatically whenever the shutter is opened or closed. An automatic computer script is used to ensure that 
the afternoon calibrations and evening observations are taken at the same spot locations. 

Several issues had to be taken into account when designing the DEIMOS flexure compensation system. The capture range of 
the actuators had to be large enough to accommodate the maximum flexure expected. The lens element mounts in the camera 
turned out to be softer than planned and other mechanical flexure also cropped up, so the initial design capture range, ±15 pixels, 
proved to be barely adequate. Image rotation is also not corrected by the system, necessitating rigid mounts for all fold mirrors 
and gratings, which can rotate the image. 

The most important limitation is that all parts of the image do not move uniformly when certain optical elements are moved, 
owing to higher-order optical distortions. Some of these distortions were anticipated, such as the distortion due to slit curvature, 
but a worse effect is that the central parts of the image move slightly differentially from the edges due to differential sags among 
the heavily-curved camera elements. Since the EC spots are at the edges of the focal plane, their motions do not precisely 
represent image motions near the center of the FOV. This last factor limits the ultimate accuracy of the system, causing image 
motions of 0.30 pixels RMS along the dispersion direction and 0.50 pixels RMS perpendicular to the dispersion. However, both 
motions more than meet their specs; the 0.30 pixel motions along the dispersion direction are 6 times smaller than the tolerance 
of 1.8 pixels that must be achieved to flat-field the 1200-line grating used for DEEP2. 

A.4. Pipeline Sky Subtraction 

Our method for modeling the sky spectrum and subtracting it depends on the fact that the flexure compensation system will 
reliably place a given wavelength on a given pixel, and thus that the afternoon wavelength calibration closely matches the wave- 
length calibration of the evening observation. We use the technique of b-spline fitting to produce a very precise model of the 
night sky spectrum. The code used is based on similar code implemented by Scott Buries and David Schlegel for a data reduction 
pipeline designed for the Sloan Digital Sky Survey fiber spectrograph. *"" Related techniques are discussed in detail in Bolton & 
Buries (2007). 

B-splines (short for "basis splines") are smooth curves that can accurately represent an arbitrary, continuous and continuously 
differentiable function given a sufficiently dense set of control points. Essentially, b-splines provide the interpolation power of 
cubic splines, but unlike cubic splines, they need not pass through the control points used to determine the interpolating curve. 
The technique has found wide application in smooth modeling of irregular high-dimensional surfaces and can filter out noise and 
reject outliers, which is very important for sky modeling. The cubic (n = 4) b-splines used in the DEEP2 DEIMOS reduction 
pipelines are comprised of localized cubic basis functions that are least-squares fitted to the measured intensity in each 2D sky 
pixel of a spectrum as a function of the wavelength of that pixel (for a general introduction, see de Boor 2001); these basis 
functions are continuous at the "breakpoints" between successive segments through the second derivative. 

For example, suppose that we have a slitlet with five arc seconds of slit-length free of any objects, and hence suitable for 
estimation of the sky spectrum. The DEIMOS plate scale corresponds to eight pixels per arc second, so there are 40 separate 
sky spectra, each consisting of 4096 pixels per chip (i.e., 4096 on the blue CCD and the red CCD separately). We tilt each 
slitlet by at least 5 degrees relative to the detector rows so that the 40 different measures of the sky brightness each sample 
the sky spectrum at slightly different wavelength. This tilting is critical, as it results in a strongly oversampled sky spectrum 
free of significant gaps in wavelength coverage. By means of KrArNeXe arc spectra taken during the afternoon, we fit a 2-d 
wavelength solution to each slitlet with a fifth-order polynomial describing the variation of wavelength with pixel number in the 
dispersion direction (i.e. along a row of the extracted spectrum), a second-order polynomial describing the average derivative of 
wavelength with pixel number along the slit direction as a function of pixel number in the dispersion direction, and a constant 
offset for each row compared to the value predicted from the polynomial terms (e.g. to correct for any shifts in the average 
wavelength at some position along the slit due to dust grains partially blocking the slit). In order to achieve a robust model for the 
local wavelength solution on each pixel, the program takes the second and higher Legendre-polynomial coefficients from the 2-d 
wavelength calibration done each afternoon, but supplements them with a linear term and zeropoint shift determined from fits to 
the cross-correlation shift between the observed sky spectrum and a template sky spectrum from the Keck HIRES spectrograph 
(Osterbrock, Fulbright, & Bida 1997) in 100 A windows, measured as a function of wavelength. 

We use 6144 breakpoints evenly spaced in wavelength for the 1-d b-spline, which is simultaneously fit with outlier rejection to 
data from all sky rows on the slitlet (after rectification to remove the spectral curvature). The b-spline minimizes the deviations 
of a locally cubic polynomial from all measurements of sky flux as a function of wavelength (40*4096 points in our example). 

See http://spectro.princeton.edu/idlutils_doc.html and http://spectro.princeton.edu/idlspec2d_doc.html for 
details. 
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The FWHM of the DEIMOS PSF of ^2.3 pixels means that the image of the sHtlet is critically sampled or better (it is the width 
of the native PSF that is important in determining gradients in the measured sky brightness, not the slitwidth itself, which is 3.9 
pixels FWHM). The number of breakpoints used (1.5 times the number of pixels along a row of the extracted, rectified spectrum) 
is therefore more than adequate to ensure full sampling.^' After an initial fit, the code then discards points that are outliers by 
more than 20-(t, and the process runs again for a maximum of three iterations. 

We use individual b-splines to estimate the 1-d sky spectrum of each frame of spectroscopic data for each DEIMOS slitlet. 
Almost all of the time, the b-spline technique produces a sky spectrum that is extremely clean and allows us to reach the Poisson 
noise limit, as shown in §10.2. As a result, it was not necessary for us to dither objects up and down along slits to ensure good 
sky subtraction, and doing so would have been a waste of detector real estate. The method breaks down when major pieces of the 
2-d spectrum are missing due to CCD gaps or bad columns, or in cases where large amounts of skyline light is scattered in from 
neighboring alignment-star boxes, but overall the system works very well, as demonstrated by Figures 24 and 23. 
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Table 1 

Selected Publications Describing the DEEP2 Survey & Science Highlights 



Paper 



Content 



Survey design, infrastructure: 
Davis et al. (2003) 
Faber et al. (2003) 
Coil et al. (2004b) 
Yan et al. (2004) 
Davis et al. (2005) 
Davis et al. (2007) 

Bimodality, red sequence: 
Willmer et al. (2006) 
Faber et al. (2007) 
Bundy et al. (2006) 
Gerke et al. (2007) 
Schiavon et al. (2006) 
Barker et al. (2006) 

Groups and clusters: 
Newman et al. (2002) 
Gerke et al. (2005) 
Gerke et al. (2012) 

Interactions, morphologies: 
Lin et al. (2004) 
Lin et al. (2008) 
Lotz et al. (2008) 

Environments: 
Cooper et al. (2006) 
Cooper et al. (2007) 
Cooper et al. (2008) 
Lin etal. (2010) 
Cooper etal. (2010) 
Cooper etal. (2012) 

Clustering, satellites: 
Coil et al. (2004a) 
Coil et al. (2006a) 
Coil et al. (2007) 
Conroy et al. (2007) 
Coil et al. (2008) 
Coil et al. (2009) 

Star formation, feedback: 
Yan et al. (2006) 
Lin et al. (2007) 
Noeske et al. (2007ab) 
Weiner et al. (2007) 
Marker (2008) 
Konidaris (2008) 
Yan et al. (2009) 
Salim et al. (2009) 
Weiner et al. (2009) 
Sato et al. (2009) 

Structural properties: 

Kassin et al. (2007) 
Dutton etal. (2011) 

AGNs and black holes: 
Nandra et al. (2007) 
Pierce et al. (2007) 
Gerke et al. (2007) 
Georgakakis et al. (2008) 
Bundy et al. (2008) 
Montero-Dorta et al. (2009) 
Comerford et al. (2009) 
Yan etal. (2011) 



Science objectives and early results 

DEIMOS integration and testing 

CFHT BRI photometry; peat catalog 

DEEP2 mock catalogs 

Survey description; early group and cluster counts 

AEGIS data sets 

Luminosity functions; growth of red sequence since z ~ 1 

Lum. funcs for DEEP2+COMBO-17; quenching of blue galaxies feeds RS 

Stellar masses and mass functions; down-sizing of quenching mass 

Quenching of red galaxies stronger in groups; onset at z ~ 1.3 

Red sequence contains recently quenched populations from z ~ 1 to now 

U — V color of RS evolves slowly due to continuing new arrivals 

Measuring w with group and cluster counts 
Group and cluster catalog, group properties 
Final group and cluster catalog 

Pair counts; galaxy merger rate is not as fast in past as thought 

Pair counts by type; separate merger rates for "wet," "dry," and "mixed" 

Gini/M2o morphologies; galaxy merger rate is not as fast in past as thought 

Color-environment relation at z '^ 1 ; bright blue galaxies in dense regions 

Red galaxies in dense environments increase continuously after z ~ 1 .3 

Sloan-DEEP2 comparison; disappearance of massive blue galaxies at z ~ 

Dry and mixed merger rates increase rapidly with density, wet merger rates do not 

Color-density relation existed but was weaker at z ~ 1, even for constant-mass samples 

The impact of environment on the size-stellar mass relation at z ~ 1 

First ^(r) at z ^ 1 for large sample; clustering trends resemble local ones 
First ^(r) for groups and clusters at z ~ 1 ; galaxies do not follow NEW profiles 
Galaxy-QSO clustering; QSOs cluster like blue galaxies, similar halo masses 
Satellite motions; halo masses, mass-light ratios, stellar mass growth 
Definitive DEEP2 ^(r) vs. color and luminosity 
Galaxy-AGN clustering; AGN cluster similarly to their (predominantly red) hosts 

[O II] in SDSS and DEEP2 red sequence galaxies is due to LINERS, not SFR 

SFR is enhanced by X 2 in close pairs 

Star-forming "main sequence", staged SFR model, downsizing in SFR 

Extinction, SFR, and AGN tracers based on emission-line indices 

Optical spectra at z ~ 0.75 consistent with the "main sequence" SFR model 

[O II] in DEEP2 red sequence galaxies is due to LINERs, not SFR 

Post-starburst galaxies; "absolute" environments constant with redshift 

Most mid-IR (24 fi) dust emission comes from older stars, not current SFR 

Ubiquitous outflows from star-forming galaxies at z ~ 1.3 

Outflows at low z more often in IR-luminous or post-starbursts 

5o.5 linewidths; TF relation already in place at z ~ 1 for normals and mergers 
Observed Tully-Fisher evolution agrees with LCDM-based models 

AGNs populate RS, green valley, upper blue cloud; persist after quenching 

AGNs are found mainly in E/SO/Sa galaxies, rare in mergers 

Dual merging supermassive black hole in early-type galaxy 

Population of deeply buried, persistent, dust-obscured AGNs on the RS 

AGN "trigger rate" matches red sequence quenching rate after z '^ 1 

Numbers and environments of red sequence Seyferts and LINERS evolve after z '^ 1 

Thirty-two candidate in-spiraling supermassive black holes 

Fractions of obscured and Compton-thick AGN at z ~ 0.6 were at least as large as today 
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Table 2 

Instrament and Exposure Parameters 



Parameter 



Value 



Slitmask parameters: 

Total slit length 

Usable slit length 

Maximum mask width 

Typical number of slitlets per mask 

Slit width 
Slidet tilts 

FWHM along slit for point source 
Spectral parameters: 

Grating 

Typical spectral range 

Dispersion 

FWHM of sky line in spectral direction 

Spectrograph image quality, point source 

Spectral resolution" 

Equivalent velocity u 

Order blocking filter 

Peak throughput, telescope+spectrograph 

RMS image stability* 
Detector parameters: 

Detector mosaic 

Pixel size 

Total array size 

Blue-red gap in spectral direction 

Gap between CCDs, spatial direction 

Detector scale 
Gain 

Read noise 
Exposure parameters: 
Nominal exposure 
Typical sky level between OH lines 
Brightest OH sky line 



16.7 

16'3 

5 '3 
140 

I'.'O 

Up to 30° 

0.5-1 .2 (best/worst seeing) 

1200 1 mm"', gold-coated, 7500 A blaze 

6500-9300 A 

0.33 A px-' 

1.29-1.39 A; center best, worse at detector edges 

2.1 px (field center), 2.7 px (corners) [FWHM] 

R = X/A\ =5900 at 7800 A 

24 km s"' 

GG495 

29% with fresh coatings 

0.3 px along spectrum, 0.5 px perpendicular 

Eight 2Kx4K CCDs in 2 X 4 array 

15^ 

8192 px X 8192 px; 12.6 cm X 12.6 cm 

7 px 

Approx. 8" 

0'.'ll85px-' 
1.25 e" per DN 
2.55 e- 

3 X 20 min. 

13 photon px"' in 20 min. near 8000 A 

'~2000 photons px"' in 20 min. 



Notes: 

° AA is the FWHM of an OH night-sky line. 

* Daily, spanning evening observations and matching afternoon calibrations. 















Table 3 
























DEEP2 Field Information 












Fid 


RA 


Dec 


I 


b 


Red." 


Area 


Area 


Pre- 


Masks 


Masks 


Targ. 


!2>3 


No 












nom.* 


done"^ 


sel?'' 


nom."" 


done^ 


done* 


z's* 


1 


14:19 


52:50 


96.52 


-1-59.57 


08-11 


0.60 


0.60 


N 


120 


104 


17,745 


12,051 


2 


16:52 


34:55 


57.38 


-t-38.63 


15-24 


0.93 


0.62 


Y 


120 


85 


10,201 


6,703 


3 


23:30 


00:00 


83.79 


-56.55 


34-49 


0.93 


0.90 


Y 


120 


103 


12,472 


8,126 


4 


02:30 


00:00 


168.10 


-53.99 


22-40 


0.93 


0.66 


Y 


120 


103 


12,494 


7,943 



Notes: 

" E(B — V) range within nominal field area in units of 0.001 mag, from Schlegel, Finkbeiner & Davis (1998). 

* Nominal area of field as originally planned, in square degrees. 
^ Area actually covered, in square degrees. 

'' Indicates whether targets are pre-selected using BRI color cuts. 
'^ Nominal number of masks in field as originally planned. 
^ Number of successful masks actually obtained. 

* Number of targeted candidate galaxies on slitmasks. These are the objects listed in Table 10; duplicate observations are included. Counts objects eventually 
found to be stars but not stars in alignment boxes. 

* Number of reliable galaxy redshifts in this field with 2 = 3 or 4, not including stars or duplicate redshifts. 
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Table 4 

Existing and Planned Data in DEEP2 Field 1 (Extended Groth Strip)" 



Data 


Wavelength 


Depth* 


Area^ 


Contact 


Keck DEIMOS (DEEP2) 


6500-9100 A 


R = 24.I 


16'xl25' 


M. Davis 


z = 0-1.4 










Keck DEIMOS (DEEPS f 


4800-9600 A 


R = 25.5 


I6'x60' 


S. Faber 


z = 0-1.4 










MMT Hectospec 


4500-9000 A 


R = 22.5 


17' X 120' 


C. WiUmer, A. Coil 


CTC-EMIR GOYA 


1-2.5^ 


K = 24.5 


~0.3 deg^ 


R. Guzman 


Chandra ACIS 


I-IOkeV 


200 ksec 


17' X 120' 


P Nandra 


Chandra ACIS (AEGIS-X) 


I-IOkeV 


800 ksec 


I7'x40' 


P Nandra 


XMM EPIC 


0.1-15 keV 


70-82 ksec 


30' diam 





GALEX UDeep imaging 


FUV.NUV 


154, 270 ksec 


1.25° diam 


S. Salim 


GALEX gri.sm 


FUV.NUV 


291, 291 ksec 


1.25° diam 


C. Martin 


Hubble ACS 


F606H',F8I4W 


28.7,28.1 


10'x67' 


A. Koekemoer, J. Lotz 


Hubble ACS (CANCELS) 


F606H',F8I4W 


28.9,28.8 


7'x26' 


A. Koekemoer 


Hubble NICMOS 


JH 


25.0,24.8 


0.013 deg^' 


S. Kassin 


Hubble WFC3 (CANDELS) 


FI25W,F160W 


27.2,27.3 


7'x26' 


A. Koekemoer 


CFHT8KxI2K(/jcaO 


BRI 


24.5,24.5,23.5-^ 


4x28'x42' 


J. Newman, A. Coil 


CFHT Megacam (CFHTLS)' 


ugriz 


ugri ~27, z ~25.5 


Ideg^ 


S. Gwyn 


MMT Megacam 


u'g'i'z' 


26.2,27.2,26.0,26.0 


Ideg^ 


M. Ashby 


LBT 


UY 


~ 25,23.9 


17' X 110' 


B. Weiner 


Subaru 


R 


27.0 


Ideg^ 


M. Ashby 


KPNO 4m NDWFS* 


BwRI 


27.1,26.2,25.8 


1.4 deg^ 


A. Dey 


Palomar 5m (POWIR)' 


.IK, 


23.9,21.7-22.5 


0.2,0.7 deg^ 


K. Bundy 


CFHTWIRCAM 


YJHK 


~25 


20'x20' 


R. Pello 


CFHT WIRCAM (WIRDS)' 


JHK 


24.8,24.6,24.5 


3x20'x20' 


C. WiUott 


Subaru 


K 


24.5 


7'x40' 


T Yamada 


KPN04 m NEWFIRM 


JK 


24.4,23.9 


1.4 deg^ 


M. Dickinson 


KPNO 4m NEWFIRM (NMBSf 


JH(meA.y 


^: = 23.4 


28'x28' 


P. van Dokkum 


IRAC GTO 


3.6,4.5,5.8,8.0 


0.9,0.9,6.3,5.8MJy 


10' X 120' 


P Barmby 


IRAC GTO ("Handle") 


3.6,4.5,5.8,8.0 


1.0,1.5,9.3, 12.0pJy 


60'x25' 


I. Huang 


IRAC (SEDSf 


3.6,4.5 


25.7,25.7 fiSy 


I2'x75' 


G. Fazio 


AKARI" 


I5At 


115-150 All' 


I0'x70' 


M. Im 


MIPS GTO 


24,70^ 


77AiJ,10.3mJy 


10' X 120' 


J. Huang, R. Hickox 


MIPS FIDEL" 


24,70,160^ 


30/iJ,3mJ,20mJy 


I0'x90' 


M. Dickinson 


Herschd PACS (HERMES)" 


110,170^ 


5.2,7.4mly 


I0'x67' 


D. Lutz 


Herschel SPIRE (HERMESf 


250,350,450^ 


lI.I,I5.2,12.9mJy 


I0'x67' 


S. Oliver 


Scuba! Legacy Deep 


850 p 


3.5mJy 


Ideg' 


R. Ivison 


VLA 


6cm 


0.6mJ" 


30'x80' 


S. Willner 


VLA 


20cm 


100^ J 


30'x80' 


R. Ivison 


GMRT 


50cm 


75mJ 


I0'x90' 


A. Biggs 



Notes: 

" Planned or in progress are in italics. For a general overview, see http://aegis.ucolick.org/astronomers.html. 

All magnitudes are AB mags. Limiting magnitudes are 5 a unless otherwise stated. 
' Areas are approximate. 

DEEP3 will acquire '-^8,000 new spectra and double the redshift sampling density in the Hubble ACS mosaic region. 
^ NICMOS: 63 pointings. 
^ CFHT peat: 8 cr. 

^ CFHTLS = CFHT Legacy Survey (http://www.cfht.hawaii.edu/Science/CFHTLS). 
^ NDWFS = NOAO Wide-Deep Field survey; EGS is an extension of the main NDWFS. 
' POWIR = Palomar Observatory Wide Infrared Survey (Conselice et al. 2008). 
^ WIRDS = WIRCAM Infrared Deep Survey (http://terapix.iap.fr/rubrique.php?id_rubrique=256) 

NMBS = Newfirm Medium Band Survey; five medium-band filters from / through H (van Dokkum et al. 2009). 

SEDS = Spitzer Extragalctic Deep Survey (G. Fazio, priv. comm.). 
'" AKARI: 50% of the field is at the two quoted depths (5 cr; M. Im, priv. comm.). 

" FIDEL = Far-Infrared Deep Extragalactic Legacy Survey (http://irsa.ipac.caltech.edu/DATA/SPITZER/FIDEL) 
" HERMES = Herschel Multi-Tiered Extragalactic Survey (http://astronomy.sussex.ac.uk/ sjo/Hermes). 
'' VLA, 6 cm: 10-cr (S. Willner, priv. comm.). 
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Table 5 

Existing Data in DEEP2 Fields 2,3,4 



Data 


Wavelength 


Depth" 


Area'' 


Contact 


Field 2 (16:52, 34:55): 










Keck DEIM0S/DEEP2 


6500-9100 A 


R = 24.1 


28'x84' 


M. Davis 


z = 0.75-1.4 










Chandra ACIS' 


1-10 keV 


9ksec 


30' X 100' 


S. Murray 


CFHT Megacam 


/,z 


24.6,23.2 


Ideg^ 


L.Lin 


Palomar/KPN04m (POWIR)'' 


JK 


22.4,21.5 (variable) 


0.2 deg^ 


K. Bundy 


IRAC^ 


3.6,4.5,5.8,8.0 


1.9,2.9,17,21 


45' X 100' 


J. Huang 


MIPS* 


24,70,160m 


0.227,40,200 mJy 


0.7 deg^ 


B. Weiner 


Field 3 (23:30, 00:00): 










Keck, DE1MOS/DEEP2 


6500-9100 A 


R = 24.1 


28' X 126' 


M. Davis 


z = 0.75-1.4 










Chandra ACIS' 


1-10 keV 


9ksec 


30' X 100' 


S. Murray 


CFHT Megacam 


/,z 


24.6,23.7 


2deg^ 


L.Lin 


Sloan Equatorial Stripe 


ugriz. 


~21 AB 


All 





Palomar/KPN04m (POWIR)'' 


JK 


22.4,21.5 (variable) 


0.3 deg^ 


K . Bundy 


Field 4 (02:30, 00:00): 










Keck, DE1MOS/DEEP2 


6500-9100 A 


R = 24.1 


28'x90' 


M. Davis 


z = 0.75-1.4 










Chandra ACIS' 


1-10 keV 


7ksec 


30' X 100' 


S. Murray 


CFHT Megacam 


/,z 


24.5,23.7 


2deg^ 


L.Lin 


Sloan Equatorial Stripe 


ugriz 


~21 AB 


All 





Palomar/KPN04m (POWIR)'' 


JK 


22.4,21.5 (variable) 


~0.25 deg^ 


K. Bundy 


CFHT WIRCAM 


J 


24.0 AB 


0.45 deg^ 


L.Lin 


IRAC 


3.6,4.5,5.8,8.0 


4,6,25,25 ^liy 


0.95 deg^ 


R. Hickox 



Notes: 

" All magnitudes are AB mags. Limiting magnitudes are 5 a unless otherwise stated. 

Areas are approximate. 
' Chandra ACIS program #9900045 (Murray et al. 2008). 
'' POWIR = Palomar Observatory Wide Infrared Survey (Conselice et al. 2008). 
^ Spitzer program #50660 (Jones et al. 2008). 
'' Spitzer program #40689. 
^ Spitzer program #40455. 



56 



Table 6 

DEEP2 Compared to Other z r 



1 Surveys 





DEEP2 


DEEP2- 


TKRS" 


VVDS- 


VVDS- 


zCOSMOS- 


PRIMUS*^ 






EGS 




deep* 


wide^ 


bright'' 




Fid name(s) 


Fields 2-4 


Field 1 


GOODS-N 


F02,CDFS 


F02, 10, 14,22 


COSMOS 


10 fields 


Mag lim^ 


24.1 (R) 


24.1 (R) 


24.4 {R) 


24.0 (/) 


22.5 (/) 


22.5 (/) 


23(/) 


Spectral res, R 


5900 


5900 


2000 


230 


230 


600 


7-100 


Nom targs* 


45,000 


17,775 


2,018 


35,000 


~100K 


20,000 


270,000 


Compl targs* 


35,214 


17,775 


1,987 


10,949 


31,196 


10,643 


270,300 


Reliab z's' 


24,785 


12,617 


1,536 


4,506 


8,961 


900 


79,419" 


7. err, km s"' 


30 


30 


;S 100 


275 


275 


55 


1800-3300 


Nom area^ 


2.80 


0.60 


0.046 


2.0 


16.0 


1.7 


10.0 


Compl area* 


2.18 


0.60 


0.046 


0.61 


6.1 


— 


9.1 


z'sdeg^' 


11,400 


21,000 


21,300 


9,400 


2,300 


11,800 


14.300 


zeffid-, 


0.711 


0.717 


0.725 


0.399 


0.280 


0.565 


0.45 


zeffid{ 


0.734 


0.725 


0.762 


0.432 


0.418 


0.592 


0.48" 


Environ merit" 


503 


— 


5.1 


33.3 


13.2 


115 


— 



Notes: 

" Wirth et al. (2004). * Le Fevre et al. (2005). " Garilli et al. (2008). '' Lilly et al. (2007,2009). 

" A. Coil (priv. comm.). 

f Magnitude limit in AB mags. 

^ Number of target slitlets in survey as originally planned. 

* Number of target slitlets completed to date. Includes stars but not serendips or (for DEEP2) objects without usable spectra (Q = -2). Roughly 1 . 1 % of DEEP2 
targets have Q = -2. 

' Reliable galaxy redshifts to date with probability of correctness being >95%, stars not included. For all surveys but PRIMUS, this means quality codes 3 and 4 
or equivalent. VVDS-wide number includes bright galaxies that are also part of VVDS-deep. 
' Total field area as originally planned, in square degrees. 

* Total field area covered to date (single-pass only in parts of VVDS-deep and VVDS-wide). 

' Reliable galaxy redshifts deg^. Includes only galaxies with quality codes 3 and 4 (no serendips or stars) and attempts to be representative for non-uniformly 
covered surveys (VVDS-wide and VVDS-deep). DEEP2 numbers use galaxies in this paper and actual field areas minus regions lost to stars, CCD gaps, and 
incomplete mask coverage. TKRS is similar. VVDS-deep estimate uses existing redshifts in Field F02 with an assumed field size of 0.4 1 deg^ ; VVDS-wide 
estimate uses numbers of published galaxies in F02, FIO, F14, and F22, which have mostly single pass; does not include the smaller area in F02, which is mostly 
double-pass. zCOSMOS-bright estimate uses the target density and star rate from Lilly et al. (2007) together with the reliable z-success rate measured from the 
zCOSMOS DR2 sample (Lilly et al. 2009) to predict a final density for 03+24 galaxies of 10,400 galaxies deg^. 

"" Overall redshift efficiency, defined as the fraction of all targeted slitlets that yield reliable galaxy and/or QSO redshifts (> 95%). Efficiency for VVDS-deep is 
calculated from the released set of 8981 redshifts, of which 35% have / < 22.5 and may also be in the VVDS-wide sample. 

" Galaxy redshift efficiency, defined as the fraction of targeted galaxies and/or QSOs that yield reliable redshifts (> 95%). Total failures are allocated amongst 
stars and galaxies proportional to their numbers amongst the objects with successful redshifts; if all failures are in fact galaxies, the redshift efficiency would be 
lower than listed here. 

" Figure of merit indicating a survey's weight for environmental and clustering measures, given by N~^ jA, where N is the number of reliable redshifts (in 
thousands) and A is field area in degrees. Values used are the existing numbers for all surveys except zCOSMOS-bright, for which the design numbers are used. 
Combined numbers for DEEP2 and EGS are used. No entry is given for PRIMUS since its redshifts are not accurate enough to localize individual objects within 
the web of large-scale structure (which requires redshift errors significantly smaller than the correlation length, <~ 500 km/s; cf. Cooper et al. 2005). However, 
PRIMUS should provide accurate measurements of larger-scale overdensity around individual objects, and can provide information on the relationship between 
galaxies and dark matter on smaller scales via projected cross-correlation functions. 
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Table 7 

Redshift Quality Codes Assigned to all Targets 



Rab mag range 


Color" 


e=-2* 


(2=-l- 


6=1'' 


2=2'' 


Q=3f 


g=4« 


Total 


<2().0 


Blue 





11 


2 


13 


6 


168 


200 




Red 


2 


29 





10 


8 


215 


264 




Total 


2 


40 


2 


23 


14 


383 


464 


20.0-21.0 


Blue 


6 


16 


5 


63 


23 


401 


514 




Red 


3 


48 


5 


45 


34 


416 


551 




Total 


9 


64 


10 


108 


57 


817 


1065 


21.0-22.0 


Blue 


9 


9 


29 


104 


64 


1020 


1235 




Red 


14 


89 


30 


122 


213 


1658 


2126 




Total 


23 


98 


59 


226 


277 


2678 


3361 


22.0-23.0 


Blue 


29 


12 


346 


396 


217 


2203 


3203 




Red 


70 


211 


306 


571 


954 


6158 


8270 




Total 


99 


223 


652 


967 


1171 


8361 


11473 


>23.0 


Blue 184 


17 


5871 


2050 


1071 


6033 


15226 






Red 


217 


504 


2327 


1812 


3634 


12906 


21400 




Total 401 


521 


8198 


3862 


4705 


18939 


36626 




Full sample 


Blue 228 


65 


6253 


2626 


1381 


9825 


20378 






Red 306 


881 


2668 


2560 


4843 


21353 


32611 






Total 


534 


946 


8921 


5186 


6224 


31178 


52989 



Notes: 

" Color: Blue = R-l< 0.5, Red = R-I> 0.5. 

* Q = -2: data so poor that object was effectively never observed. 

^ Q = -l: star. 

''2=1: probable galaxy but very low S/N; data not likely to yield redshift. Many of these are distant galaxies beyond the nominal redshift limit of z - 

' Q — 1: low S/N or data are somehow compromised for reliable redshift; reason given in comments. 

^2 = 3: reliable redshift with probability of accuracy >95%. 

^2 = 4: reliable redshift with probability of accuracy >99%. 



1.4. 



Table 8 

Quality Codes for Duplicate Observations 



2 


-1 


1 


2 


3 


4 


-1 


43 










1 


4 


264 








2 


2 


186 


52 






3 





48 


54 


68 




4 





82 


114 


295 


1038 



Note: Entries represent the number of duplicate pairs with the indicated combination of quality codes. 
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Table 9 


















Mask Data (Sample) 










Mask No. 


Obs. Date 


RA 


Dec 


PA 


Nobj 


<5/A'> 


Seeing 


%(tot) 


%(red) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1100 


2003-05-03 


213.74732 


52.07412 


-48.86 


152 


0.718 


0.677 


66.23 


73.85 


1101 


2003-05-06 


213.79093 


52.09807 


-48 


86 


155 


0.348 


1.017 


53.64 


56.79 


1102 


2003-06-30 


213.82309 


52.12490 


-48 


86 


158 


0.897 


0.737 


76.47 


81.54 


1103 


2003-05-03 


213.86179 


52.14789 


-48 


86 


150 


0.746 


0.705 


71.62 


81.82 


1104 


2003-05-30 


213.89094 


52.17260 


-48 


86 


162 


0.355 


0.767 


66.24 


69.51 


1105 


2003-07-01 


213.92597 


52.19302 


-48 


86 


147 


0.555 


1.177 


73.47 


78.82 


1106 


2003-05-05 


213.95854 


52.22014 


-48 


86 


151 


0.449 


0.999 


65.56 


72.62 


1107 


2004-05-22 


213.99695 


52.24292 


-48 


86 


156 


0.340 


1.028 


65.16 


67.90 


1108 


2004-05-21 


214.02546 


52.26719 


-48 


86 


153 


0.515 


0.928 


71.71 


78.12 


1109 


2003-05-06 


214.06047 


52.28759 


-48 


86 


154 


0.187 


1.080 


52.67 


57.97 


1110 


2003-05-04 


214.08909 


52.31194 


-48 


86 


159 


0.439 


0.798 


62.99 


66.23 


1111 


2003-05-28 


214.12444 


52.33257 


-48 


86 


155 


0.734 


0.823 


75.48 


84.81 


1112 


2003-05-04 


214.15082 


52.35534 


-48 


86 


157 


0.413 


0.979 


64.47 


63.22 


1113 


2003-05-30 


214.18373 


52.37425 


-48 


86 


155 


0.477 


0.957 


67.74 


72.29 


1114 


2003-06-01 


214.20993 


52.39691 


-48 


86 


156 


0.688 


0.740 


71.61 


71.43 


1115 


2003-05-05 


214.24393 


52.41659 


-48 


86 


155 


0.451 


0.996 


55.26 


56.99 


1140 


2003-05-05 


213.76523 


52.12802 


41.14 


150 


0.360 


0.864 


53.10 


48.48 


1141 


2003-05-06 


213.80431 


52.10401 


41.14 


142 


0.450 


0.893 


53.62 


55.07 


1142 


2004-04-21 


213.84704 


52.08417 


41.14 


150 


0.535 


1.047 


70.07 


72.46 


1143 


2003-05-05 


213.88612 


52.06016 


41.14 


143 


0.400 


1.008 


57.86 


74.29 


1144 


2003-05-30 


213.92886 


52.04032 


41.14 


143 


0.552 


0.670 


74.29 


84.62 


1145 


2003-05-30 


213.96794 


52.01631 


41.14 


156 


0.435 


0.720 


64.00 


67.11 


1146 


2003-05-06 


214.01067 


51.99647 


41.14 


151 


0.503 


0.788 


61.38 


56.60 


1150 


2003-05-04 


214.04399 


52.32691 


41.14 


147 


0.554 


0.764 


70.83 


72.41 


1151 


2003-05-28 


214.08325 


52.30289 


41.14 


146 


0.575 


0.674 


74.83 


74.19 


1152 


2003-05-04 


214.12617 


52.28306 


41.14 


145 


0.469 


0.784 


65.49 


62.32 


1153 


2003-06-01 


214.16543 


52.25904 


41.14 


141 


0.873 


0.743 


73.91 


73.91 


1154 


2004-05-22 


214.20835 


52.23921 


41.14 


153 


0.610 


0.863 


75.00 


80.00 


1155 


2003-06-30 


214.24761 


52.21519 


41.14 


152 


0.848 


0.616 


82.55 


87.06 


1156 


2004-04-19 


214.29054 


52.19536 


41.14 


150 


0.410 


1.159 


63.45 


66.25 



Note. — (1) Mask number: the first two digits represent the CFHT pointing and the last two digits are a position code within the 
pointing (see Section 13); (2) UT observing date in YYYY-MM-DD for the night contributing the majority of a mask's data; (3)-(4) 
RA/Dec (2000.00) of mask center; (5) PA of long axis of mask; (6) number of target objects on mask not counting ahgnment stars; 
(7) median FWHM of alignment star spectra in arcsec (seeing measurement); (8) continuum SjN per pixel near 6900 A(based on 
median flux and measured variance, omitting the atmospheric B-band region); (9) percentage yield of reliable redshifts for overall 
target sample on mask {Q > 3); (10) percentage yield of reliable redshifts for redder galaxies (R — I> 0.5) only. 



Table 10 

Galaxy Data, Part 1 (Sample) 
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ID 


RA 


Dec 


B 


R 


/ 


o-B 


o-R 


o-/ 


rg 


62 


PA 


Pgal 


E(B-V) Mb 


U-B 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


11001673 


213.86870 


51.95644 


23.49 


23.14 


22.58 


0.026 


0.022 


0.029 


1.370 


0.046 


-46.95 


1.000 


0.011 


-999 


-999 


11001699 


213.81047 


51.94232 


22.07 


20.03 


19.55 


0.009 


0.001 


0.002 


1.957 


0.026 


-0.40 


3.000 


0.011 


-19.23 


1.130 


11001770 


213.84843 


51.94888 


24.14 


24.10 


24.02 


0.068 


0.077 


0.152 


2.184 


0.332 


-9.43 


3.000 


0.011 


-999 


-999 


11001800 


213.83176 


51.95255 


25.34 


23.51 


23.08 


0.169 


0.033 


0.050 


1.796 


0.069 


-2.61 


0.510 


0.011 


-999 


-999 


11001860 


213.83255 


51.95417 


24.38 


23.40 


22.57 


0.079 


0.034 


0.035 


2.023 


0.048 


-12.92 


3.000 


0.011 


-20.08 


0.752 


11001861 


213.81732 


51.95325 


23.25 


22.82 


22.21 


0.035 


0.028 


0.035 


2.823 


0.114 


-38.11 


3.000 


0.011 


-20.22 


0.532 


11001878 


213.81459 


51.93851 


23.97 


23.42 


23.13 


0.056 


0.032 


0.053 


1.844 


0.077 


62.49 


0.663 


0.011 


-17.67 


0.395 


11001898 


213.81851 


51.96013 


24.18 


23.14 


22.60 


0.090 


0.040 


0.054 


3.049 


0.359 


-81.78 


3.000 


0.011 


-18.75 


0.646 


11001909 


213.82000 


51.94717 


23.10 


21.18 


20.52 


0.025 


0.004 


0.005 


2.059 


0.037 


-24.73 


3.000 


0.011 


-21.02 


0.813 


11001922 


213.82240 


51.93072 


22.01 


21.26 


21.06 


0.014 


0.008 


0.016 


2.415 


0.089 


21.47 


3.000 


0.011 


-14.92 


0.563 


11001927 


213.82467 


51.96038 


25.29 


22.05 


20.96 


0.208 


0.012 


0.010 


2.547 


0.145 


-28.47 


3.000 


0.011 


-20.43 


1.253 


11001934 


213.82196 


51.95469 


22.89 


22.14 


21.83 


0.022 


0.013 


0.021 


2.430 


0.108 


43.00 


3.000 


0.011 


-19.66 


0.417 


11001950 


213.83030 


51.94160 


23.60 


22.41 


22.12 


0.052 


0.017 


0.028 


2.445 


0.295 


-75.92 


3.000 


0.011 


-999 


-999 


11001974 


213.80832 


51.94306 


23.05 


22.47 


21.92 


0.019 


0.012 


0.016 


0.799 


0.053 


-2.44 


0.999 


0.011 


-20.59 


0.464 


11001978 


213.78782 


51.95398 


24.46 


23.85 


23.55 


0.064 


0.042 


0.070 


1.339 


0.096 


36.55 


0.769 


0.011 


-999 


-999 


11002016 


213.79909 


51.95631 


24.10 


23.91 


23.44 


0.046 


0.044 


0.064 


1.664 


0.042 


78.09 


1.000 


0.011 


-19.49 


0.460 


11002019 


213.78484 


51.94730 


23.99 


22.75 


22.26 


0.052 


0.018 


0.025 


1.909 


0.197 


-43.56 


0.429 


0.011 


-19.15 


0.639 


1 1002024 


213.75619 


51.93770 


24.12 


23.60 


23.07 


0.069 


0.040 


0.053 


1.979 


0.395 


-19.89 


3.000 


0.012 


-18.69 


0.516 


11002039 


213.78086 


51.95369 


25.60 


24.09 


23.36 


0.183 


0.058 


0.059 


1.644 


0.042 


72.13 


0.980 


0.011 


-999 


-999 


11002051 


213.77227 


51.94480 


24.61 


23.20 


22.40 


0.110 


0.032 


0.033 


2.262 


0.192 


87.64 


3.000 


0.011 


-999 


-999 


1 1002064 


213.73228 


51.93605 


24.63 


23.52 


23.22 


0.131 


0.048 


0.077 


2.336 


0.437 


-23.33 


3.000 


0.012 


-999 


-999 


11002085 


213.76673 


51.95897 


23.55 


23.33 


23.05 


0.032 


0.030 


0.051 


1.898 


0.062 


61.91 


0.941 


0.011 


-999 


-999 



Note. — See Section 13 for details. Columns are: (1) Unique DEEP2 object ID frompcat; (2)-(3) object RA and Dec (2000.00); (4)-(9) total BRI magnitudes and errors from 
peat; (10) Gaussian radius (11)-(12) ellipticity and PA from proprietary team peat; (13) probability of being a galaxy determined from photometry; (14) Galactic E(B — V) from 
Schlegel, Finkbeiner & Davis (1998); (15)-(16) absolute B magnitude (Mb -5 ■ logjQ/i) and rest-frame (U—B) from Willmer et al. (2006) K-corrections (sometimes labeled (U-B)o 
in DEEP2 papers). 
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